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ERRATA 
line 4 from bottom; I should be 1. 


„ line 11 from bottom of text; before 


"Or see," insert footnote reference,2. 
bottom; add footnote B:"See 
Supplement PF." | 

line 9; Integral should have lower 
limit, ty. . 

prob. 25; xp should be xs. 

line 11 from bottom; delete "the." 
line 5 from bottom; delete "for 
each arc," 

line 4 from bottom; delete "differ- 
ence for each arc, and hence the 
correct range." 

line 20; after "theoretically," add 
fooinote reference, 3. 


„ bottom; add footnote 5: "Зее p. 63." 
» footnote, last eq.; Up 


Should read: (уе + Е 

line 15; for "& function of," read 
"dependent on." 

line 17; for “equations 47," read 
"equation 47." 

line 7 from bottom; for "book," 
read "chapter." 

line 6 from bottom; for "some in 
Chapter XI," read "that for site." 
ae 5 from bottom; for "Ag" read 
" * 


line 19; instead of "whence," the 
line should read: "whence, by 
theorem 9, page 213" 

line 9 from bottom; add footnote 
reference, 2. 

line 4 from bottom; add footnote 
reference, 3. 

bot tom; add footnotes: 2 See p. 44," 
апа "3 See Problem 55, pe 69, post." 
line 2; first word should be "in." 


Р. 67, 
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Р. 88, 
Р. 88, 
Р. 93, 
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P.113, 
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P. 117, 
P. 126, 


last line; to read: "radius of the 
earth, and was the angle of fall." 
lines 1, 2; change "x - X" io "X - 
x" in each line. 

lines 12, 13; change "As" to " Ўв" 
in each line. - 

fig. 6; the "J" which marks the angle 
should be "j." 

first eq.; "(E 4 J)" should be 

"(E + j)." 

eq. 79; delete the comma in the 1019. 
ед. 94; (x'Q - ути )" should be 
Je Mx )." 

eq. 95; change "cos & to "віп ос" in 
first eq. only. 

85; column to right of the T column 
should have had all horizontal lines 
broken, to indicate that several 
columns have been omitted. 

No. 56; Integral should have upper 
limit, "T." 

first eq.; "A" should be "Mtg." 
values of A, B, С; lower limits on 
integrals are б. 
value of А; "(х'р = y'v )" should be 
"(Vy -px')." 

line 7 from bottom; integral should 
have upper limit, "T" 

fig. 13; reference should be to fig. 
12, instead of to fig. l. 

footnote; the word "second" should 
not be italicized. 

line 3 from bottom of text; delete 
the footnote reference. 

line 1; line should read: "Substitut- 
ing from equations 117 and 118 in 
equations 116;" 

last 5 lines; change "dx" to "dX," 
four times. 

footnote; the first A in the equation 
should have for subscript, "t - 1." 
col. 2; "tA" should be "tA " 
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EXTERIOR BALLISTICS. 
INTRODUCTION. 


The work. of the ballistic computer is divided into three parts: 
(1) the computation of the elements of standard trajectories; (2) the 
computation of differential corrections, whereby the elements of a 
standard trajectory may be corrected for nonstandard conditions; 
and (3) the utilization of the foregoing to construct range tables 
from firing records. 

As is shown in Chapter VI, the World War ushered in & new era 
in the handling of ballistic problems. "This may be called the period 
of numerical integration. 

The approximations, which had been used to modify the simple 
Newtonian equations of motion into such form that they could be 
formally integrated, gave place to precise numerical integration of 
these equations in their original form. 

Practically parallel progress was made in all of the allied countries 
during the war. 

In America the great step in the computation of TN was 
the introduction of numerical integration. Numerical integration 
had long been used in astronomical calculations, and so it was natural 
that an astronomer and mathematician, Prof. F. R. Moulton, of the 
University of Chicago, while serving as a major in the Ordnance Depart- 
ment of the United States Army, should have applied this method to 
ballistics. The practical work of this method was materially reduced 
by formulas for ‘‘integrating ahead“ later introduced. 

The work can be still further reduced by a variant of Maj. Moulton’s 
method, known as the tangent-reciprocal method. But this method 
has the pedagogical drawback of obscuring the physical meaning of 
the steps involved, and hence will not be given first place in this 
book.! 

Another improvement is the change in the analytic interpretation 
of the equations of motion. Formerly the z-axis was conceived of 
as tangent to the earth at the gun, the system being Cartesian. In 
the modern conception, the x of a point is measured along the curved 
surface of the earth, and the y is measured vertically from this surface. 

The progress in the computation of the differential corrections has 
involved more steps. Differential equations for the corrections were 
devised at the same time that numerical integration was introduced. 


1 Bee Supplement A. 
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But these had the difficulty of requiring an independent computa- 
tion for each correction. 

This difficulty was removed by the discovery at Aberdeen of a 
method of solution by means of a set of adjoint equations involving 
several auxiliary variables. All the differential corrections can be 
expressed in terms of these auxiliary variables. A physical deriva- 
tion of these variables, and of the corrections based upon them, was 
_ later found. Matters were further simplified by reducing all of these 
variables to expressions in terms of one variable and its derivatives. 
A physical derivation of this step was at once forthcoming. 

One further step should be noted, namely the development of the 
weighting-factor curves for zero elevation, which have been of great 
value in interpolation. 

The present method of constructing range tables out of firing 
records is a logical result of substituting for the Ingalls tables the 
new methods of computation. Tables, to take the place of the com- 
putations now necessary, are now being constructed by the tech- 
nical staff at Washington. 

The credit for the above-described development is largely due to 
Mr. J. J. Arnaud, Master Computer, Ordnance Department; Prof. 
A. A. Bennett, of the University of Texas, then Captain, Ordnance 
Department, U. S. Army; Prof. G. A. Bliss, of the University of 
Chicago, Technical Expert, Aberdeen; Mr. Philip Franklin, Com- 
puter, Aberdeen; Dr. T. H. Gronwall, Mathematics and Dynamics 
Expert, Technical Staff; Prof. Н. H. Mitchell, of the University of 
Pennsylvania, Master Computer, who organized the range table 
computation work at Aberdeen; Dr. J. F. Ritt, of Columbia Uni- 
versity, Master Computer, Technical Staff, and their associates, in 
addition to those mentioned elsewhere herein. 

Most of the written material on the ballistic progress made during 
the war consists of scattered blue prints, some printed at Aberdeen and 
some at Washington. These pamphlets overlap in spots, contain 
some hiatuses, and do not agree їп symbology and nomenclature. 
Through a three-cornered correspondence between the Technical Staff, 
Ordnance Office, (War Department), at Washington, D. C., and the 
Ordnance School of Application and the Ballistic Section, Aberdeen 
Proving Ground, Md., a uniform symbology and nomenclature have 
been established as standard. 

The first course of instruction in these new ballistic methods ever 
given in this country was given at the Ordnance School of Applica- 
tion in the winter of 1919-20 by Capt. Roger Sherman Hoar, Coast 
Artillery, then in charge of the Ballistic Section of the Proof Depart- 
ment at Aberdeen. This present book is based upon the papers 
used in that course, and uses the standard symbology and nomen- 
clature established as above. 
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It is assumed that the student is thoroughly grounded in algebra 
and plane trigonometry, and knows enough calculus to appreciate 
the meaning of a derivative, a differential, and a definite integral. 
On that basis, this book gives, in Chapters I to IV, the irreducible 
minimum of higher mathematics necessary to understand all points 
involved in the later chapters. 

‘The book then takes up in succession: An introduction to modern 
ballistic methods (Chaps. V and VI); the computation of trajecto- 
ries (Chaps. VII and VIII); the computation of differential correc- 
tions (Chaps. IX to XV); and the construction of range tables (Chap. 
XVI). Alternative methods, elaborations of certain points, and a 
brief mention of the more involved mathematical processes necessary 
to the computation of antiaircraft range tables are reserved for 
supplements. 

Each chapter is followed by a series of questions, designed to bring 
out the salient features of the chapter. The answers to most of 
these questions will be found categorically stated in the text, but 
some under each chapter will require a small degree of original 
thought on the part of the student. 

Throughout the book the attempt is made to explain as much as 
possible from the viewpoint of physics rather than from the view- 
point of abstract mathematics. 


CHAPTER I. 
PARTIAL DIFFERENTIATION. 


Before defining the partial differentiation" of a function, let us 
define the word function.“ u is called a function“ of z, y, z, etc. 
if, when z, y, 2, etc., are given, the value of u is determined. Note 
the broadness of this definition. Thus u=zy can be regarded as a 
function of z, y, and 2, although it is evident that w is not in the 
least dependent on z. 

The functional relation is expressed in the general form: 


* (x, Y, 2, ). 
One should be sure to notice the fundamental fact that this ex- 
pression, as it stands, does not take into consideration any relation- 
ship which may exist between any of the variables in the parenthesis; 
in other words, these variables may or may not be independent in 
that expression. 

Usually the functional equation can be altered so as to make z 
explicitly a function of u, y, 2, etc.; y explicitly a function of v, т, 2, 
etc. In case such a conversion is either impossible or even merely 
inconvenient, it is better to regard the equation, in its original form, 
as defining u as a function of т, y, 2, etc., or z as a function of u, y, 2, 
etc.; and to differentiate it as it stands, using the differential method. 
Such an equation is that in problem 3, to follow. "Thus, in an equa- 
tion containing n variables, any one of these variables can generally 
be regarded as dependent, and the remaining л—1 as independent. 

We are now in a position to define the term “partial derivative." 
If u=f (т, y, 2 . . . ), then the partial derivative of u with 
respect to х is obtained by treating as constants all the other variables 
in the parenthesis, and differentiating with respect to z, in the ordi- 
nary manner. 


But the partial derivative is written ES instead of B so as to in- 


indicate that z is one of several independent variables, instead of being 
the sole independent variable. Thus оч means the rate of change 


of u with respect to a change in z alone out of several independent 
variables. 

An important point to note in this connection is that the symbol 
2v is highly ambiguous, i. e., it has entirely different meanings, and its 
partial derivatives have entirely different values, according to what 
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are regarded as the independent variables. To avoid this ambiguity, 


the subscript notation should be used. Thus Qus means that v 


should be expressed as a function of z, y, and 2, and then differentiated 
with respect to z: Sar means that v should be expressed as a func- 
tion of X and Y before differentiating, etc. 

In practice, the subscripts may be omitted whenever it is self- 
evidént, from the problem, just what are to be treated as the inde- 
pendent variables in each differentiation. 


PROBLEMS. 
Oz дг 


Oy 
Oz! oy ^^ 


Find = d 55 in each of the following: 


(1) 2=2 log y. 
(2) z=ar -+ 2bry — сд. 


e 
G) eap 
(4) 22+ +22 = а. 


All of the theorems of partial differentiation can be derived as 
special cases of the following general theorem: 


1. Fu is a function of X, Y, Z. - , each of these in turn 
being a function of z, y, z - - +, then (with certain assumptions as to 
continuity): 

OUxys... uz. ; Хув... ‚у О%хп... ; OV ays. э. 
Ox OX d oY or 


and similarly for 5 ete. 


This theorem can be derived by the use of undetermined coeffi- 
cients, or see Osgood, p. 296. 

This theorem is called general,“ not in the sense that it is the 
fundamental basis of the other theorems, for many of them have a 
simpler derivation, and some, in fact, may serve as'steps in the deri- 
vation of theorem 1. Nor із it advisable to use theorem 1 when some 
simpler formula is available. But, as theorem 1 has a form easy to 
remember, and as any other formula can be derived as & special case 
thereof, it serves as & good memory peg on which to hang the whole 
subject of partial differentiation. 

The student should particularly note at this juncture that, although 
any partial derivative may be obtained by the differential method; 


1 For further exposition of the subscript notation, see Osgood, Differential and Integral Calculus, p. 306. 
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yet, once formed, its numerator and denominator are inseparable and 
neither can be canceled. Witness the absurd results which would 
follow from performing all possible cancellations in theorem 1. 

The following general principles may be used in deriving special 
formulas from theorem 1: | 

(a) Whenever in partial differentiation any given variable is re- 
garded as dependent on one independent variable alone, then, in the 
expression for the derivative of the former with respect to the latter, 
the operator 0 should be changed to the operator d. | 

(b) The derivative of a variable, with respect to another variable 
of which the first is regarded as independent, is zero. 

(c) The derivative of a variable with respect to itself is unity. 

(d) If, for all values of its variables, a given function is explicitly 
a constant, then the derivative of that function with respect to any 
of such variables is zero. 

The following special formulas may be derived from theorem 1: 


ди ди ди 
2. du yv X t3 ydY +5742. 
This is called the expression for the total differential" of u. 
9. Fu be a particular function of X and У, namely XY, then: 


ди ОХ oY 
Эх“ op TAG! 


which is the analog of the following formula of total differentiation: 
du- YdX + Ха F. 


4. Ju equals it then: 


— — — — — — — — 
— 


oz Z Or 7 Or Лох X Ox Y ox ZOx 


Similar special formulas can be derived for other special relations: 
which may exist between u, X, Y, Z, etc. 


PROBLEMS. 


Nore —The following problems should be treated exactly as though the symbols 
y’, x’, E, etc., were the simpler-looking symbols of the preceding problems and. 
explanation. For the purposes of any given differentiation, х, 2, £”, y, у, y^, 
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and E have no relation to each other except that given in the hypotheses. But the 
work on these problems should be carefully saved for later reference. In Chapter V 
@ meaning will be given to each of these symbols, and in Chapter IX it will be seen 
that these problems, taken in order, constitute almost the entire proof of some im- 
portant ballistic formulas. 


PROBLEMS. 


(5) Given that z” and y” are each functions of the independent 
variables 2˙, y’, and y, evaluate dz” and dy" in terms of dx’, dy’, 
and dy and some partial derivatives. 

For problems 6 to 13, the following is given: 


Ex 
y" = — Ey'—g 


where g is a constant; and E is a function“ of т”, у’, y and z, although 
not dependent on z. 


Oz! oy". Oz'' ду’! 
б 52 oS (8) 5— 09) SE 
Ox“ òy” Òr” oy" 


(10) G (11) or (12) oy ` (13) oy 


The foregoing problems can be done either by one of the special 
theorems or by theorem 1. It is advisable not to use theorem 1; 
but if it be used, care should be taken to observe that x’ and y’ 
each enter into the expressions for z” and / in two capacities: 
1. e., as а variable of the X sort and as a variable of the z sort. To 
ilustrate this, perform the following problem: 


(14) 2''= Е2 
Е =f (z, Y’, y) 
= р! 


Evaluate 92. 
т 


Also by theorem 5. 


-QUESTIONS ON CHAPTER I. 
1. What is a partial derivative? 
2. How, in general, is partial differentiation performed? 


Oz ду Oz 
3. Does dy Or equal oz ? If so, why? If not, why not? 


4. If w=f (p, q, т), and oe is evaluated in terms of p, g, and r, and we then learn 
that a certain special relationship—extraneous to the equation w=/ (p, q, r)—exists 


between p, q, and т, will this fact change the value of ои» If the special relation- 


ship is such that w equals a constant, will that fact make or If so, why? If 


not, why not? 
5. Explain the meaning of, and the need for, the subscript notation. 


CHAPTER II. 
SUCCESSIVE APPROXIMATIONS. 


An equation which determines the. numerical value of a quantity 
may generally be expressed in a variety of alternative ways. Thus 
the fact that z is the square-root of w may be expressed: 


(а) х= + Jw 


(b) z— 7-0 


L9 
(c) 1 


y" 


Here (a) defines z ezplicitly; (b) defines z implicitly; and (c) consists 
of simultaneous equations in z and y, each defined as a function of 
the other. 

Forms analagous to (a) are not always forthcoming. Consider, for 
example, the determination of а number z, such that = is equal to 500, 
plus 1,000 times its own trigonometric sine, x being measured in 
minutes of arc. There is no known explicit form for z. We have, 


however: | 
(b) х=500-++1,000 sin z. 


S LI 
V |y 1,000 sin x 


х= 500 y 
= Sin 2 


An implicit equation may frequently be replaced by a system of 
equations (similar in general to equations of the sort c) so chosen as 
to be convenient for solution by a computational procedure known as 
“successive approximations." This, as its name suggests, is а method 
of starting with an apt number, largely arbitrary, and by successive 
substitutions securing a sequence of approximate evaluations of z, 
approaching the precise value to within any desired degree of precision 

14 
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For example, extract the square root of w by the following equa- 
tions: 


Expressed as a formula: 
Уа 
1 
Тат = 2 (Yn + Tn) 
Extracting the square root оѓ 2 by this pair of equations, first 


assuming 2 to be approximately its own root, gives the following suc- 
cessive values: i 


n 2 у 
1 2 1 

2 з Ы 

3 it 17 
4 205 


and so on. Expressed in ordinary language, we have the following 
rule for the extraction of square root: Divide the number by an ap- 
proximation to the square root desired; the arithmetic mean of the 
divisor and quotient is a new approximation. 

It can be shown that if any approximate square root checks with 
the preceding approximate root to n figures, then the new approxi- 
mate root is correct to at least 2n-1 figures. 

Had a negative value been taken initially for z,, the negative 
square root of 2 would have been approached. The separating value, 
zero, causes the method to fail, as one would expect. 

Àn attempt to extract the square root of 2 by the equations: 


_1 
у=: 
r= wy 

produces the following results: 


Nn 


m сә М н 
— d e 8 8 
2 — te ч 
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and soon. A similar repetition would occur for any initial value for 
z, other than zero or infinity. This shows that not all sets of simul- 
taneous equations are adapted to solution by successive approxima- 
tions. 

The method of successive approximations has the advantage that 
the accuracy of each step is independent of the accuracy of the pre- 
ceding step. A mistake in any single step, therefore, while it may 
prolong the work, will not vitiate the final result. 


PROBLEMS. 


(15) Extract the square root of 100, taking 25 as the approximate 
square root and carrying each division to three decimal places. Con- 
tinue until two successive answers check to three decimal places. 

(16) Assume that 1.41459 is an approximation to the square root 
of 2. Perform one step of getting the root with greater precision, 
and give the result to only the number of places certain to be correct. 

(17) Replace z= 500 +1,000 sin т, by the system: 


y = 2,000 sin z 
2=5 (1,000 +y) 
1 22 
and use as formulas for successive approximation: 
| Уһ = 2,000 sin 23 
лн = (1,000 + Yn) 


z is expressed in minutes of arc. Find z and y correct to the nearest 
unit. | 

(18) y is a tabular function of the time, t. The following is а tabu- 
lation showing the value of y corresponding to each of certain values 
of t. 


t y a 
60 1, 569. 4 
62 1, 052. 4 —517.0 
64 512.5 —539. 9 
66 —46. 6 —559. 1 


а, b, and c are, respectively, the first, second, and third differences of 
y; i. e., the a of any line is obtained by subtracting, from the y of 
that line, the y of the line before; the b of any line is obtained by 
subtracting, from the a of that line, the a of the line before, etc. 
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To find the value of t corresponding to some given, non-tabular 
value of y: 
Let у, — the value of y at time, t; 
y, = the value of y at the time, to; 
and 4 = the tabular interval in t. 


Representing by At the value of — which is to be considered with 
its algebraic sign, the formula to be used is: 


= У: — Vo 
6. А о 5 +40 FAD TAJ | 


or ot 3| 


where a, b, and c are the values taken from the same line on which t, 
and y, occur. The At so obtained must be multiplied by 7, to obtain 
—%. 
The formula is solved by successive approximation, the first 
approximation being 
At- 2t 9» 
a 


Find the vidus of t corresponding to y, = O, using ,=64. Check by 
using t, = 66. 

This formula, obtainable from the result of problem 30, page 31, 
may be used to interpolate і in either direction from any tabulated 
values, but it requires the use of “receding differences;" i. e., differ- 
ences that, in the method of writing used above, occur on the same 
line with to and y,.! 

At the close of Chapter IV problems will be given involving а com- 
bination of numerical integration and successive approximations. 


QUESTIONS ON CHAPTER II. 


1. Define ‘‘successive approximations.”’ 

2. Give the rule for extraction of square root. 

3. Is any set of simultaneous equations solvable by successive approximations? 
4. What is the test of solvability? 

Б. What are the two chief advantages of this method? 


1 When it is required to interpolate forward from the first item in a table, е. g., downward from the t= 60 in 
the table above, the formula for “advancing differences" must be used with the corresponding values of 
у, a, b, and c; e. g., in this case, y= 1,569.4; a= —517.0; b= —22.9; andc=+3.7. The “advancing difference 
formula may be obtained from 6 by changing the sign of every At throughout the equation, and of the 
odd-ordered differences, а and c. See Supplement Н. 
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CHAPTER III. 
EFFECT OF DIFFERENTIAL VARIATIONS. 


This chapter deals with the mathematical determination of the 
effect of a disturbance, on the subsequent motion of a particle which 
moves (except for the disturbance) according to some definite differ- 
ential equations of motion. | 

Let us consider a particle moving in time, along a plane trajectory. 
‘Trajectory’ is here used in the general sense of the path of any moving 
particle, rather than in the specialized ballistics sense of the path of 
a projectile, or the still further specialized sense which will be the 
meaning employed in later chapters, namely, the path of a projectile 
moving under certain so-called standard conditions as to atmosphere, 
wind, gravity, etc. 

Reverting, then, to the motion of a particle along its trajectory, 
it is evident that at any given instant of time (¢) there will be a 

^ F uniquely corresponding value 

of z, y, z', у, £”, Y”, etc., 
where z and y are the coordi- 
nates of the particle, z’ and 
с у’ the two components of 
velocity, x“ and y” the two 
components of acceleration, 

A etc. Primes are thus seen 
FIG. | to represent time derivatives, 

and will be used in that sense throughout this book. In this chapter 
the general symbol и may be used in place of z, y, x’, у’, z”, and /“, 
in theorems true as to any of them. Hence v, also, is a function of t. 

In the illustrative examples of this chapter, t,x coordinates and 

t,y coordinates will frequently be employed, in order that the student 
may become accustomed to considering the elements of the tra- 
jectory as separately plotted against time, and to using time deriva- 
tives. 
. In the case of motion of the sort which will be considered in this 
book, and in fact in the case of most motions, the values of z”, , 
and higher derivatives are, in the absence of disturbing capses, de- 
termined, for any instant t, by the values of z, y, x’ and , or some 
of them; and hence x”, y”, etc., need not be discussed for the 
present. 

The path of the particle will, of course, be a single curve, which may 
be graphed by plotting y against z. But а chronological record of its 

18 
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motion may be represented more completely by four curves, obtained 
by plotting z, y, x’, and y’, respectively, against t. 

Consider now any one of these four curves, represented by ABC 
in figure 1. | 

Suppose & disturbance from B to E, so that the curve takes the 
shape BE during the interval of disturbance, and suppose that 
thereafter there is no further disturbance and the form of the curve is 
EF. This is the general case. The amount of the total disturbance 
up to any instant is the difference in ordinates between the original 
undisturbed curve ABC and the disturbed curve BEF, such as HI 
or GE. 

For convenience, disturbances may be treated as of three sorts: 

(a) Those disturbances which produce a finite effect in a single 
instant; as, for instance, if the curve took the shape ABGEF. 

(6) Those disturbances which vary during the total time of dis- 
turbance, so that if the total time be divided up into an infinite 
number of equal infinitesimal time intervals (dta), an infinitesimal 
рагі (dôu) of the total disturbance (ôu) will occur during each dia, 
and the total disturbance may be represented as: 


u= { dôu, or as 


T dóu 
ou = Я di, dh 


where t is the time the disturbance starts and T is the time it enda. 
This is the most general case. 

(c) Those disturbances of which a proportional part ooours during 
any part of the total time. These may be regarded as a special 
сазе of b. | 

ади Ci dta 


óu Ci ( т— AP 


where c, is some constant.! 

It is essential that one element of the trajectory be considered as 
remaining unvaried, so as to furnish a basis for measuring the varia- 
tions of the other elements. Accordingly time will be selected for 
this purpose. Both the time (t) of which the other elements (z, y, 
x’, /, z”, yY”, eto.) are functions, and the time (ta) at which a dis- 
turbance occurs, will be considered as unaffected by the disturbance. 
Accordingly we can say that: 

ôt =0 

dt, =0 
г As here given, c is a special case of b. But it is possible to regard both b and а as special cases of c. 
Thus, if the eanstant с, becomes the variable Fi we have case ò; whereas, if the time interval T—4 


becomes infinitesimal and ди remains finite, we have case a. 
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and, differentiating: 


4. 


Let us first consider disturbances of the first sort. A particle is 
moving through space according to some definite law of motion, 
expressed by differential equations. At a given instant of time 
(ta to) an instantaneous disturbance takes place, which changes 
the value of z, y, z’, and y’, or some of them. Thereafter, the par- 
ticle proceeds according to its original law of motion, of course not 
on а continuation of its original regular curve, but on another such 
regular curve, also satisfying the original differential equations. 

| € he amount of the instantaneous changes 

p. in z, y, z', and y’ will be designated, re- 

spectively, by ôx, dy, dz’, and dy’. Any one 

of these expressions can represent either 

FIG.2 a positive or a negative change. The 

changes which will be considered in any practical application of the 

principles of this chapter, and the effects resulting therefrom, are so 

minute in comparison with the elements affected by these changes, 

that the following restrictive definition can be given: 62, dy, eto., 

are small finite increments of z, y, eto., and are so minute that 

second and higher order terms (such as ёбх ог ôzôy, for instance) 

are of no consequence, as compared with z, y, etc., and hence any 

terms containing them may be disregarded and dropped from a sum 
or series in which first order terms ocour. 

The operator 6 is the operator employed in that branch of mathe- 
matics known as the calculus of variations, but not always with the 
above restriction. No further understanding of the principles and 
methods of this branch than here given is necessary to the purposes 
of this book. 

A certain resemblance between the operator ô and the familiar 
differential operator d will be noted. The distinction between the 
two should also be noted. Consider a curve in an z,y coordinate sys- 
tem. At any point on this curve, dy represents a continuous infini- 
tesimal change in y along the curve, corresponding to an infinitesimal 
change in z. dy is the slope of the curve at the point in question. 
ôy and ò represent a very small finite break in the curve in question. 
dy and dz may be considered as taking place in an infinitesimal period 
of time (dt); whereas бу and dz may be considered as taking place 
instantaneously, or as cumulating during a finite time interval. 
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Let us now derive some of the basio theorems relative to the 
operator 6. 
8. Independent changes in u and u' may be made at any time ta. 
Proof: Consider the curve AB in 2, t coordinates, and W any 
given point thereon. (See Fig. 2.) 
The curve can be moved up or down, thus Балое the х of W, 
without changing its ¢ or its slope. Or the curve can be rotated 
about W, thus changing its slope, without changing its t or its 2. 


Thus ôz and ô 105 ; are independent. Q. Е. D. 


Similarly, by plotting y against t, and y against z, it can Фе shown 
that ôr, 5x’, бу, and ôy’ are all independent. 
9. Terms containing more than one 8 may be dropped. 

This was one of the fundamental hypotheses of the definition of 
the operator 6 on page 20. 

It will now be demonstrated that the four general formulas of 
differentiation still hold true when the operator 6 replaces the opera- 
tor d; in other words, that the small finite increments of this chapter 
obey certain laws, already familiar in form for the case of differ- 
entials, although, of course, these increments are quite different 
from differentials. 

These four theorems are as follows: 


10. : ^. 6 (cu) c bu. 

11. ô (u +v) = du 4- 0v. 

12. $ (uv) =v bu+u dv. 

13. 5 (50 nes ôv 
v v 


Only one of these (namely, 13) will here be proved, the proof of 
the other three being reserved for problems. 
Proof of theorem 13: 


Take the expression —, and give u and v each an increment. Then: 


p! 


E — — —— XXL — фә 


I8 1 би u v du Ai bv 


Expand this fraction by dividing the numerator by the denom- 
inator, as follows: 


v би -u ёр v би — u oV ov ôu u ôV Óv 
v? +v dv v? v? v 


ЖО 


by theorem 9. Therefore: 
uN v du—u bv 
(0 ä Q. E. D. 


Next let us derive the expression for ‘‘total increment," analagous 
to the expression for ‘‘total differential" (formula 2 of Chap. I). 
14. Vu is a function X and Y, then: 


Ou Ou 


bu 6X T x 


. “Oe 
Proof: Let u=f (X, Y). Give X and Y the increments &X and 
8 F, respectively. Then: 
au f (X 4-0.X, Y +8Y)-f (X, Y). 
Subtract and add the quantity f (X, FTF). Then: 
zu = (X +8X, Y -6Y) -f (X, Y +Y) +f (X, Y 49 Y) -f (X, Y). 


Applying the law of the mean (see Osgood, p. 230) to each of these 
two differences gives: 


ри Oaths XY D, y 42 (X, rth lV y, 


Now, if these two partial derivatives are continuous, each would 
approach the corresponding partial derivative of f (X, Y) if ö A and 
6Y both were to approach zero, and hence will differ but slightly 
when óz and Sy are very minute. 

Consequently we may express 


dX +X, Y +Y) „ NL, Y) | , YX, Y), 
Uem o — шоу: — nec 


and 
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From the right member all terms, except the first, may be dropped, 
of( X, Ү+6,8Ү) a s AX Р) Y), 5 (X, Y) Y) 

oY ӘУ 


Substituting these values in the expression for ёи, and substituting 
4 for f (A, Y), gives: 


ш= Du X +24 bY && sou 7 57 


from which the last two terms can be dropped by theorem 9. 
An extension of this derivation gives: 
i= Dy * ЗҮ 295.5 Ze Q. E. D. 


15. (bu) = 20и) 
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Proof: Consider a particle moving along a ourve from 4 to Р. 
When it reaches P, at time і, let 2 and y instantaneously receive an 
increment (ôx and ôy, respectively) which will place the particle at 
р, and thereafter let the particle move undisturbed along the curve 
-pD (the curves AP and pD being defined by the same differential equa- 
tions of motion). Let Q be a point on the curve AP, such that, if 
no disturbance had taken place at P, the particle would have reached 
after a small finite time 
interval At from the time it 
left P. Let q be the point on 
the curve pD reached by the 
particle, At seconds after 
leaving p. Then a pair of 
changes ór--Aóz and бу+ 
Aéy, occurring at time é, + 
At, would produce the same 
-effeot as the pair of changes 
dr and ôy occurring at timet, ^ FIG. 3 

Let óz and dy, although arbitrary at time , be thereafter con- 
sidered as restricted by tne condition that at any instant thereafter 
their value must be such as to produce at that instant the same 
situation as would have existed at that instant, had the changes been 
made, with their initial values, at time 4. From this point of view, 
the average rate of change of ôr and dy, during the interval At, is 
respectively A and Ay 
time of disturbance (ta). 


Ate =R M—PG—-GM-—PR =pr-PR; 
Ady =Lq — Mr =rq— ML =rq— RQ. 


Let At approach zero. The rates of change at time f, thus become 
dôr d dày | 
dt, ^75 аһ 


D 


- At is here regarded as a change in the 


Now, the unaffected z-component of velocity (z’) at time 11 is the 
limit of A as At approaches zero, At being here regarded simply as 
a change in the time (t) of which the elements of the trajeotory are 
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functions. The value of z’ at time ti as affected by the variations, is 


A Therefore: " 
, „WL 
827) (27-58 "At J^ di 
Similarly: 


uu RQ\ d 
ку) =lim( 74- A )- -% QED. 


16. The operator 5 does not always signify an independent incre- 
ment. 

Examples: In the expression for the total increment given under 
theorem 14, if any three of the four increments are conceived of as 
independent of each other, then the fourth must of necessity be 
dependent upon the other three. 

In equation 38 in Chapter IX, 5X may be taken as a constant, and 
or as dependent on бу, dz’, and dy’. 

Throughout this book, ôt and 5t4 will be taken as zero, by theorem 7. 

So we may say that changes represented with the operator 6 are 
essentially independent, unless rendered otherwise by some express 
condition of the problem confronting us. 

Let us now consider the following problem. If a particle, moving 
according to some definite law up to time їл, suffers a disturbance 
and then moves on according to the original law, without further 
disturbance, what is the relation between the value of its z at some 
later time T and the value which z would have had at time Tif there 
had been no disturbance? Let X represent the undisturbed value of 
z at time T, and let X +ôX represent the disturbed value. The 
problem is to express 6X in terms of the 6z, dy, dz’, and dy’ occurring 
at time ta. 

If y, x’, and у’ are not changed at time ta, then 6X is of the form 
L ёх, plus terms containing more than one $, which can accordingly 
be dropped by theorem 9. J will have a value which will depend on 
the trajectory in question and on the point on that trajectory at 
which the change àz occurred. Thus for any given trajectory L is a. 
function of t4, but is not dependent on óz. 

Similarly, if z, z’, and у’ are not changed at time ta, the resulting 
6.X can be expressed as Méy, etc. 

If all four of z, y, z', and y’ are changed at time ta, the resulting ô X 
wil equal Lóz-- Móy-- Nóz'-- P dy’, plus terms containing combi- 
nations of ёт, dy, etc., which may therefore be dropped by theorem 9. 

Thus we have, as the fundamental equation for the X-effect at 
time T, due to & set of small arbitrary changes 1 in z, y, z’, and у’ at 
lime ta: 


17. 6X Lor T Мӧу+ Nóz' + Póy' 
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Formal expressions for L, M, N, and P may be obtained from theorem 
14, but have not yet been shown to be of any practical value. 

The fact the ôr, бу, ôr’, and dy’ can each be given any arbitrary 
value at any time of change (ta) enables us, if we wish, to assign any 
arbitrary value we please to 6X and any three of these, and then 
satisfy the equation by a proper choice of value for the fourth. 


18. 5X 48 not a function of ta. 


From the viewpoint of theorem 15, the values of ö, dy, eto., 
although initially arbitrary, are regarded as restricted so as to change 
value at such a rate, during the interval 
di,, as not to alter the resulting value of 
6X. . E. D. | 

Thus far, we have been considering 
merely instantaneous changes. Let us 
now consider the second sort of changes 
listed at the beginning of this chapter, 
namely changes which vary throughout a 
finite time, starting as zero at the begin- 4 
ning of the interval. A velocity change of this latter sort is readily 
seen to be made up of an infinite number of infinitesimal changes 
of velocity throughout the interval. Similarly a coordinate change 
of this sort is seen to be made up of an infinite number of infinitesi- 
mal changes in position throughout the interval. 

Consider the finite time interval, from time f, to time Т, as divided 
up into an infinite number of infinitesimal time intervals, dta. An 
infinitesimal part of the total ёт, dy, dz’, and dy’ occurs during each 
of these infinitesimal] intervals, and thus causes an infinitesimal part 
of the total ö X which will occur at time T. Thus the result, at time 
T, of the disturbance during any interval аёл, is 


19. dà X = Ldoz4- Mdóy + N dêr’ +P diy’. 


It should be noted that déz, etc., are here used in quite a different 
sense from that of theorems 15 and 16. Неге dôr means the change 
which óz undergoes along the actual disturbed path of the particle, 
during the interval dta. But here also, d and ô may be shown to 
be commutative, as in theorem 15. 

Consider z plotted against t. Let ABC, Fig. 4, represent the t, 
curve of the undisturbed trajectory. Let a disturbance start at B 
changing the shape of this curve to BED. Let EF (=GH) represent 
one of the infinitesimal] time intervals, dta. 
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The normal undisturbed rate of change of z is HI. The actual 
disturbed rate of change of z is FJ. Thus: 


HI-da 
ЕЈ -—dz--5dx 


Let EK be parallel to GT. 
The total бт at the beginning of the dt, interval is GE. The total 
ôx at the end of the interval is JJ. Thus: 


dór—-l1J—GE-IJ—IK-KJ 
$02 = FJ — HI = FJ— FK= KJ 


аёх dz 
28022) Q. E. D. 


The right-hand member of formula 19 can now be transformed by 
transposing d and б, and by both dividing and multiplying each term 
by dta, and by applying this principle of commutativeness. 


20. dô X —-Lóz'dt,-4- Мбу'ай + N oz''dt^ +P y dt a 


But: 


Therefore: 


Integrating this from & to Т, we get 


T T T T 
21. 5X - | L óz'di4 + [ M by’ dt, + [ Мзш'йһ+ | Р ày''dt4 
lo to to fo 


This is the formula for the effect, at time Т, of changes which vary 
throughout the finite time interval from f, to T. 

For solution of any of the integrals, the dz’, dy’, ö“ or ) therein 
contained must be replaceable by something which is constant either 
in value or in algebraic form throughout the interval. 

Let us now consider the third sort of changes listed at the begin- 
ning of this chapter, namely, changes which are proportional to 
time. Consider a ôr of this sort. Then, from formula 19: 


dó X —Ldóz. 


Substituting dt, for dor, and integrating within the limits of the 
disturbance, we get: 


T 
* | L dta 
b 


If L is a constant, this becomes 


5X -L(T—1). 
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PROBLEMS. 


(19) Prove theorem 10. 

(20) Prove theorem 11. 

(21) Prove theorem 12. 

(22) Prove that ô sin u=cos u ôu and that ô cos u= sin u du. 
Suggestion: Turn to some book on the calculus, and follow the 
analogy of the similar differential formulas. 


QUESTIONS ON CHAPTER III. 


. What is the object of this chapter? 
. What is the meaning of the operator 5 in ballistics? 
. Do formulas involving the operator d hold true with respect to the operator 5? 


. When would zy represent the slope of a curve? 


. In this chapter what is meant by the word “‘trajectory "? 
. Distinguish between the two meanings of dôu and ódv, 
. What variations of time are considered in this chapter? 
. Is óu an infinitesimal? 
. Distinguish between óu and du, 
0. Interpret M of formula 17 by means of formula 14, 


bt O очо m © М ны 


CHAPTER IV. 


FINITE DIFFERENCES. 


Integration by finite differences is based upon the principle that 
inasmuch as a derivative is a rate of change, the value of an integral 
of any smooth function can be computed step by step, if successive 
values of its derivative are known at sufficiently close intervals. 

Suppose that u is a function of t, such that u can not be integrated 
with respect to t by means of any of the expressions tabulated in any 
table of integrals. In other words, formal integration is impossible. 

Let us now imagine a plotted curve with u for its ordinates and t 


ten 
for its abcissas. Then udt is the area between the curve, the 


t axis, and the ordinates at f= t and f= tn. Such as area exists for 

every continuous curve, and hence every continuous curve has an. 
integral, even though that integral 
is not expressible by means of the 
ordinar, elementary functions. 

In such а case, various approxi- 
mations are possible, some crude 
and some so refined that they can 

©-AXIS produce results to any desired de- 
gree of precision. 
FIG. 5 | All of the methods here described 
require that и be first tabulated for successive finite values of t. 

Let us consider the portion of the curve |, ing between -i and tk, 
assuming the curve to be concave downward. 

The entire required area can be divided up into sections such as 
this. If we can evaluate each section, the sum of these evaluations 
will be the value of the whole. 

A first approximation would be the area under the chord. 

Thus: | 


му, HU 
dee (Кы — -I) 


If we space ¢ with unit intervals, then: 
22. A, = T Un 


This approximation is too small. 
A second approximation would be the area under the tangent at the 
point (ux, fx). Thus: 


A. ur ua А-4) 
28 
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If we space ¢ with unit intervals, then: 
23. A 97 2м 


‘This approximation is too large, but the error is about one-half the 
error of A,. А, the true area, lies between A, and A,. Thus: 


A,>A>A, 


If the curve were concave upward, the above inequality would be 
reversed. In either case it is evident by inspection that: 


„2А, T A, 


4 


Therefore a third, and very close, approximation is: 


24. ＋ A, 
3 


This is known as Simpson’s rule. 


1 
24. А, 7$ (u,-, + Au, + ик). 


PROBLEMS. 


Evaluate the following to four decimal places, by formal integra- 
tion, and by formulas 22, 23, and 24, and determine the percentage 
of error in each case: 


(23) [= ах 


e) | 


In using Simpson’s rule, when the tabular interval is other than 
unity, the formula is: 


h 
25. А, = 3 (D + 4u, + Мул) 


1 Formula 24 can be proved as follows: Any ordinary function of one variable can be expressed as a 
powereeries of that variable. Thus the general curve which we have been considering so far in this 
chapter can be expressed: 

t= Dot Pit p p pd. 


This equation represents a straight line, another straight line, a quadratic, a cubic, а quartic, ete., 
according as we drop all but the first, all but the first two, etc., terms of the right member of the 
equation. For all approximations up to and including a cubic, formula 24 is precise. Proof: By this 
formula, | RE 


ENTE. 
f tidie (tio ＋ 4 + tan) 
0 
By formal integration. 
зь 
й pf ри 
- ＋ 2 аена] 
Í Pot 2 3 4 ]t—2h 


These are identical. With A sufficiently small, similar pairs of expressions, involving quartics or 
higher, are nearly mutually identical, being indeed coincident in as far as the first four terms are con- 
cerned, the discrepancy in the higher terms being very minute. Q. E. D. 

For another proof, see Osgood, '' Differential and Integral Calculus," 1917, pages 406-408. 
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where & is the tabular interval. The area of the next section will be 


3 (uni + Aus TU) and so on, so that the area for a series of sec- 
; i h 
tions will be 3 (и, + 4и, + 21 T 44,44 -2U,4g * + + Au, HU). 


PROBLEMS. 
1 
(25) Divide ji zdz up into ten sections, where 2=х— х=, Tabu- 


late z against z, and integrate by Simpson's rule. Integrate formally, 
and compare the results. 
(26) wis a function of t: 


t 4 
110 (210) 


Tabulate the values of u corresponding to ¢=0, 1,2... 9, 10. 
10 ` 
Evaluate Í и di by Simpson’s rule and by formal integration, and 


compare the results. 
1 
(27) Evaluate Í, E 000 log 10 (1000 1000 dt, either by formal 


integration, or by tabulating from a denary log table and then using 
Simpson's rule. 

Simpson’s rule is a method of numerical integration, as distin- 
guished from formal integration. It is one of the simplest of a 
system of rules that may be obtained, involving the values of the 
function and its differences. Let us now derive a method of numerical 
integration which employs finite differences of various orders. 

We will suppose a tabular function of ¢ and will tabulate its first 
second, third, and fourth receding differences as follows, when each 
difference as tabulated is obtained by subtracting the element on 
the line above from the element on the same line of the preceding 


column, i. e., where a,=f,—f_,; bo d- a4; etc. 
{ s(t) First dif- | Second dif-| Third dif- | Fourth dif- Е 
ference. ference. ference. ference. 
—4 p 
—3 J-a: а. 
—2 S a_2 b2 
-1 fA а; by cA | 
0 Jo a bo Co d, 


з This rule should be followed regardless of the order in which the function is tabulated, i. e., regardless 
of whether t increases or decreases down (or from left to right across) the page. Thus, if the tabular func- 
tion of ¢ algebraically increases as one goes down (or from left to right across) the page, the first difference 
will be positive; if decreasing, it will be negative. Similarly, if the first difference algebraically increases. 
the seccnd difference will be positive; if decreasing, negative, etc. P 
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These differences are called ''receding," for the reason that the 
differences tabulated on any line would have receded if they had 
been tabulated opposite the space between the two elements of 
which they are the difference; sometimes, for clearness, they are во 
tabulated; but it is generally more convenient to tabulate them as 
above. | 

(28) Successively evaluate , f-z H, and f, in terms of f,, a, 
by, с, and de. 

(29) From this deduce a formula for f_n. 

(30) Substitute ¢ for —n. The result is the usual interpolation 
formula in terms of receding differences, whereby f (f) can be calcu- 
lated with values for any fractional value of f. The resulting f (t) will 
lie on à smooth curve within the interpolation interval and, if the 
true f (t) be smooth, will closely approximate it. 

(31) What is the interpolation formula for ¢=4; for t= — 41 

(32) In the tabulation of problem 27, what is the logarithm of 
10.55? Use either interpolation formula from the preceding problem. 

(33) Integrate the formula of problem 30, between the limits —1 
and 0. 

(34) Integrate it between 0 and 1. 

(35) Integrate it between —1 and 1. 

The solution of the last three problems gives us, respectively, the 
formula for integrating by finite differences and two formulas for . 
integrating ahead. These are the fundamental formulas of the 
method of computing trajectories by numerical integration, which 
method will be discussed in Chapter VIII. 

(36) Form the differences of the tabulation in problem 26 and 
integrate each interval ahead, using the formula of problem 34 in 
the earlier stages, and that of problem 35 in the later stages. 

(37) Integrate this same tabulation across, by using the formula 
of problem 33. Compare the two sets of results.? 

(38) Sum up the results in problem 37 and compare with the two 
answers already obtained for problem 26. 

3 In numerical integration the values of the first few integrals will be inexact, due to the absence of first, 
second, third, etc., differences in the integrand. То supply this lack extrapolate back for the missing dif- 
ferences from such first, second, third, etc., differences as are obtainable from the tabulation. A better 
method would be to rewrite the first four or five of the tabulated functions in reverse order, and difference 
this new series as usual. What was the first interval being now the last, the set of differences to be used 
in integrating over this last interval by the formula of problem 33 are now obtainable. The figures now on 
the last line will be found to be the figures that were at the top of the columns of the origínal tabulation, 
except that the signs of the differences of odd order are now changed. By the method here suggested the 
student can construct for himself formulas for interpolating and integrating at the beginning of a table 
with these advancing differences. When this is done, the rewriting of the function in reverse order may be 


discontinued. Опе of the methods of procedure outlined in this note should be followed in this and in all 
succeeding problems of numerical integration. See Supplement H. 
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(39) Substitute 1 for ¢ in the interpolation formula of problem 30. 
The result is an extrapolation formula. 
(40) Using the values .: % do, fefe and the value of f, 


obtained by the preceding problem; evaluate | f dt, | fdt, and 
-1 0 


f i f dt by Simpson’s rule, obtaining the values of E and f} by 
—1 


the formulas of problem 31. Compare the resulting formulas with 
those obtained by problems 33, 34, and 35. 

(41) Derive formula 6 used in problem 18 of Chapter II. 

(42) The last tabular values of t, у, a, b, and c for a time interval 
of one second are | 


— —ä—ͥ— ————————À ͤ ä́6ö6— | ee 


— ——— 0 - - ng me —— — 


Find the value of ¢ corresponding to y=0, using the formula of 
problem 41. Interpolate for y as a check, using the formula of 
problem 30. Carry ¢ to four decimal places. 

(43) The following tabulated values of y” and z/' from an actual 
trajectory computation are given: 8-inch rifle, railway mount, 
model of 1917. Muzzle velocity (V), 594.36 m/s; ballistic coefficient 
(O, 4.0; angle of departure (6), 10°. The meaning of these terms 
need not be considered at this stage. 

Values of x“ and y” at two-second intervals: 


r” t y” 
—368 | 0 | —16.3 
—29.4 | 2 | —18.9 
—23.6 | 4 | —12.2 
~18.9 | 6 | —10.9 
—15.1 | 8 | —99 
—1.9 | 10 | — 9.3 
— 9.4 | 12 | —8&9 
— 7.5 14 | — 8.6 
— 6.3 | 16 | — 8.4 
— 5.4 | 18 | — 8.2 


Obtain у’, x’, y, and z for each value of t by numerical integration, 
checking the final results by Simpson’s rule. Initially z=0, х'= V 
cos ; y=0, y’= V sin G. 
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(44) From this data, calculate values for t, z, х’, and у’ corre- 
sponding to y=0. Ав a check, get the value of y corresponding to 
this value of t. 

The formulas derived in problems 28 to 35 inclusive, and in prob- 
lem 38, are to be used only with receding differences, and should be 
plainly so marked in the student's notebook.* 


QUESTIONS ON CHAPTER IV. 


1. Define integration by finite differences. 

2. State ‘“‘Simpson’s rule. 

3. How are first, second, etc., differences formed in the case of a tabular function? 

4. What is the distinction between ‘‘numerical integration’’ and ‘‘formal inte- 
gration’’? 

5. From theoretical considerations rather than from the numerical results which 
happen to have been obtained in any of the numerical problems, what is your opinion 
of the relative precision of the six methods of numerical integration given in this 
chapter? State detailed reasoning. 

6. State the three formulas for integration by receding differences. 

7. When should each be used? 

8. What is meant by receding differences? 

9. Give two methods of integrating at the beginning of a table. 


* For methods of deriving advancing difference formulas, see Supplement Н. 
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CHAPTER V. 
ELEMENTS OF THE TRAJECTORY. 


From the point of view of the present state of ballistics, the follow- 
ing are the more important elements and features of a trajectory: 

Trajectory.—This term, as here used in connection with computae 
tions, embraces only the standard trajectory of the projectile, moving 
in accordance with assumptions which may be laid down as follows: 

(1) The earth is motionless. (The average eífect of the earth's 
rotation on gravity is included in the assumed value of g.) 

(2) The gun and target are in the same altitude above sea level. 

(3) The preassigned standard muzzle velocity for that type of gun 
is actually obtained. 

(4) There is no wind. 

(5) The atmospheric density varies regulariy with the altitude 
according to the exponential law assumed and is standard at the 
muzzle (15° C.=59° Е. 720 mm. of mercury, 78 per cent saturation; 
1.2034 kg/m’). 

(6) The action of gravity is uniform in intensity, is directed 
towards the earth’s center, and is independent of the geographical 
location of the gun, g=9.80 m/s’. 

(7) In the computation of standard trajectories, the velocity- 
resistance factor, G (v), may be regarded as dependent only on the 
velocity; the standard variation in elasticity of the air being ac- 
counted for by the value of the ballistic coefficient adopted. 

(8) The ballistic coefficient is a constant on the trajectory and is 
as determined from experimental firing. This includes the assump- 
tion that all precessional and nutational effects of the projectile 
may be ignored in computation, and hence the trajectory as computed 
lies in а vertical plane. 

(9) The vertical jump of the gun is as чайга from observation. 
_ (10) The atmospheric temperature (and hence elasticity) varies 

regularly with the altitude according to the law assumed, and is 
standard at the muzzle (15° C. 259? F.),! 
| 1A question arises in connection with assumptions (7), (8), and (10) as to the use of the temperature 
structure. The standard temperature structure adopted by the Ordnance Department is supposed to 
approximate average conditions aloft. Firing conditions on different occasions can be accurately compared 
only when reduced to the same standard temperature structure, and for convenience the above-mentioned 
standard temperature structure is always used as the basis of comparison. Variations in range due to irreg- 
ularities in temperature sloft are computed by reference to this standard temperature structure. But the 
G function used in the computation of trajectories for standard conditions is assumed to be dependent only 
on the velocity. This amounts to assuming, for the trajectory computation, that the temperature is the 
same for all altitudes. The С function was based on a number of actual firings which were not corrected 
for variations in temperature. These firings were probably made under average surface conditions, that is, 
for temperature approximately 15° C.—59° Е. The computed trajectories involve a ballistic coefficient, 
C, which is so taken as to make the range check with firings corrected to standard temperature structure 
aloft. The explicit introduction of variable temperature in the original computations, while possible 


for individual trajectory computations, is not regarded as important, and is not feasible in the case of the 
ballistic tables. 
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(11) The drift (including lateral jump) is as assumed from observa- 
tion. 

Coordinates (x and у). — The coordinates of any point on the 
trajectory, measured in meters. The abcissas (х) are measured 
along the surface of the earth and are positive in the direction of fire. 
The ordinates (y) are measured vertically from the surface of the 
earth and are positive upward. The origin is the muzzle of the gun. 
In the development of formulas and in computation, this coordinate 
system is treated as cartesian, the error being negligible. Note 
the difference between this conception and the tangent-plane con- 
ception of the Ingalls-Siacci ballistics. The present conception, 
being based on a curved earth, obviates the necessity for correcting 
for curvature of the earth. 

Surface of earth. A spherical surface passing through the muzzle 
of the gun and concentric with the earth. 

Muzzle velocity (V).—A fictitious initial tangential velocity of the 
projectile at the beginning of its flight in meters per second. The 
blast may continue to accelerate the projectile for some distance 
beyond the muzzle; so that the ‘‘muzzle velocity” is not the actual 
velocity at the muzzle, but is rather a fictitious velocity which, if it 
occurred at the muzzle and if there were no blast, would cause the 
projectile to travel on the same trajectory as that on which it 


. actually travels. 


Velocity (v).—Sometimes called remaining. velocity." The tan- 
gential velocity of the projectile at any point of its flight, in meters 
per second. 'The z and y components of the velocity are repre- 
sented by z' and y' respectively. 

Acceleration (x or y’’).—The rate of increase of x’ and /, respec- 
tively (in meters per second per second), at any point on the trajec- 
tory. 

Time (t).—The time in seconds elapsed in the flight of the pro- 
jectile from the muzzle to any point on the trajectory. tis the inde- 
pendent variable of the trajectory. (When considered as the time 
when a disturbance takes place, the symbol is fa.) 

Quadrant angle of departure ($).— The angle measured from the 
horizontal to the tangent to the trajectory at the muzzle; sometimes 
called the angle of projection. 

Inclination (0).—The angle measured from the horizontal to the 
tangent to the trajectory at any point on the trajectory. 

Point of fall (under range table assumptions) .—The point where the 
projectile in its downward flight reaches the same altitude above sea 
level as the muzzle. 

(Geographical) range (.X).—The distance in meters from the muzzle 
to the point of fall, measured on the surface of the earth. 

Quadrant angle of Tal (). The. negative “ the inclination at the 
point of fall. 
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Retardation (R).—The retardation due to the resistance of air of 
standard density and elasticity. H is a function of velocity and alti- 
tude. 

R=vE 
СП 
COME 

E 1s called the resistance function; G = a tabular function of v, 

for convenience tabulated to argument ——; C is the ballistic coeffi- 


a 
cient; H is determined by the following exponential law: 


Н = e 
h = .0001036 


Quadrant elevation.—The angle between the horizontal and the 
axis of the bore just before the gun is fired. 

Vertical jump.— The algebraic difference obtained by subtracting 
the quadrant elevation from the quadrant angle of departure. 

Summit.—The highest point of the trajectory. 

Maximum ordinate (%). The y coordinate of the summit. 

Time of flight (T).—The time (t) from the muzzle to the point of 
fall. 

Angle of site.—The angle whose tangent is the ratio of the difference 
in altitude of gun and target divided by the range. The angle is 
positive if the gun is higher than the target. Owing to the fact that, 
in the present ballistics usage, zero altitude is considered as being 
the same elevation above sea level as the gun, rather than as lying 
in the plane horizontal at the gun, the angle measured directly with 
& transit or level is only approzimately equal to the angle of site, 
although satisfactory for short ranges. 

Point of splash.—The point where the projectile, in range firing, 
enters the water. 

Center of impact.—The average position of several points of splash. 

Ballistic coefficient (C).—A purely empirical number, used in the 
formula for the resistance function (E). С is a mean value (con- 
stant over a given trajectory) of the reciprocal of the relative retarda- 
tion. The relative retardation is the ratio of the retardation experi- 
enced by the actual projectile to that which would have been experi- 


enced by a certain fictitious standard projectile," moving at the 


same velocity and altitude. С for any given projectile and muzzle 
velocity varies only as a function of the angle of departure. C may 
be represented as: 
а= 16. 
id? 
where w is the weight of the projectile in pounds, d its diameter in 
inches, and 7 the “coefficient of form," so-called because it is largely 
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dependent upon the form of the shell. This is ап empirical number, 
its value being determined by giving the C in the equation above the 
value necessary to make the range computed with that C equal to 
the actually observed range, when the latter is reduced to standard 
conditions. 

Ascending branch.—That part of the trajectory in which the pro- 
jectile rises. . 

Descending branch.—That part of the trajectory in which the 
projectile descends. 

Plane of projection.—The vertical plane including the line of pro- 
jection, i. e., the tangent to the trajectory at the muzzle. 


PROBLEM. 
(45) Prove that— . 
х’ =v cos 0 
y' = 0 sin 6 


QUESTIONS ON CHAPTER V. 


1. State the eleven standard trajectory assumptions. 

2. Which of these assumptions are always in error? 

3. What corrections in practice have to be applied to compensate for variations 
rom each of the ten assumptions? 

4. Define E, and each of its elements, and explain how each is obtained. 

5. Draw a trajectory lying in the plane of fire and label the following: Gun, surface 
of earth, angle of departure, elevation, angle of jump, maximum ordinate, range, point 
of fall, angle of fall, ascending branch, descending branch, summit. Such as have 
symbols may be labeled by the appropriate symbol. 

6. Draw another trajectory. Mark a dot on it to represent the projectile in flight 
at the end of t Seconds. At that point, draw an arrow to represent the velocity. Draw 
arrows to represent its components, x’ and . Label х, y, and 0. 

7. Define each of the elements of the two preceding questions. 

8. Why does the assumption of a curved earth (see trajectory assumption No. 1, 
‘*Coordinates” and Surface of the earth”) make it impossible to measure the angle 
of site precisely with.a transit or level? 

9. Why is the assumption of a curved earth (see ''Coordinates") more convenient 
than the assumption of a flat earth? 


CHAPTER VI. 
HISTORY OF EXTERIOR BALLISTICS. 


The science of ballistics consists of three parts, namely: Interior 
ballistics, dealing with the behavior of a projectile in the gun; exterior 
ballistics, dealing with its behavior during flight; and ballistics of 
penetration, dealing with its behavior while entering the target. 
This book deals only with exterior ballistics. 

The fundamental differential equations of motion of a projectile 
in flight have been known since the time of Newton. These are based 
on the two components of acceleration, which may be represented as 
follows: 

Horizontal acceleration = — Ё cos 6 
Vertical acceleration = — E sin 0—0 


in which Ё is the retardation due to the resistance of the atmosphere, 
0 the inclination (to the horizontal) of a tangent to the trajectory, 
and g the acceleration of gravity. 

These equations look very simple, and would be if E were а constant 
or were one of certain simple expressions, in terms of time, velocity, 
z,y,or0. But H appears to follow no mechanical law which can be 
algebraically expressed. All that we know about Ё has been derived 
empirically, and even the simplest available approximate expression 
for А renders formal integration out of the question. 

Two great obstacles to the development of exterior ballistics have 
always been (1) lack of knowledge about Ё and (2) the difficulty of 
solving the differential equations. On the basis of the methods pur- 
sured in attempting to overcome these obstacles, the progress of 
ballistics may be divided into three periods. 

The first period, namely all the years prior to about 1865, may be 
called ‘‘the algebraic period." During this period, attempts were 
made to represent 2 by some simple algebraic expression, and thus 
to solve the equations by artifices such as those found in the stand- 
ard methods of calculus. At first it was generally assumed that the 
motion of a projectile through the air satisfied the conditions laid down 
by Newton, under which А would vary as v?. The experiments of 
Robins and Hutton, and the two sets of firings at Metz, furnished the 
first real data on R at velocities sufficient for practical ballistics. As 
a result, № was variously supposed to vary as the square or as the 
cube of v, or as some combination of the two; and, based on these 
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suppositions, there were ‘devised many extremely ingenious solutions 

of the equations, all of which solutions were found in course of time to 

be insufficiently approximate, and most of which were very laborious. 
The second period, lasting from about 1865 to the beginning of the 

present war, may be called ‘‘the period of approximations."  Practi- 

cally all of the physics and mathematics of this period were based 

upon Mayevski’s formula: 

Ат" 

С 


in which È is the retardation, v the velocity, and C the ballistic 
coefficient. It was assumed that velocities could be divided into a 
few intervals (such as 0 to 790 f/s., 790 to 990 f/s., etc.), each interval 
with its own law of resistance, and each with its own constant value 
for A and for ». The experiments of this period, notably those of 
Bashforth and Mayevski, and the Meppen and Gávre firings, were 
accordingly directed to finding how these intervals might be taken, 
and what the proper constants would be for each interval. The 
solutions of the differential equations of motion (notably those of 
Bashforth, Mayevski, Zaboudski, Hojel, Siacci, Didion, Braccialini, 
and Ingalls), attempted during this period, not only made use of the 
Mayevski formula, but also, in developing the various equations, 
used frequent approximations, most of them based upon the assump- 
tion of a practically flat trajectory. Thus the equations were finally 
wrenched into a solvable form. Our familiar Ingalls tables were pro- 
duced in this manner. 

. The third period, extending since the beginning of tho World War, 
may be called ‘‘the period of numerical integration." Since the time 
of Euler it had been known that, by the use of numerical integration, 
the differential equations could be integrated in their original form 
without resort to approximations. But the approximate methods 
were always thought to be simpler and sufficiently precise, until 
certain phases of long-range and high-angle fire introduced during 
the World War compelled a resort to numerical integration. Thus 
ballistics has reverted to the first principles. 

The chief differences between the Ingalls-Siacci ballistics and 
present ballistic methods may be summarized as follows: 

The former used a datum plane tangent to the surface of the earth 
at the gun. The latter uses the actual surface of the earth as datum. 

The former assumed air of a uniform density, this density being 
chosen at such an average for each arc as to give the correct range 
difference for each arc, and hence the correct range for the point of 
fall, but not giving the other terminal elements precisely, and pro- 
ducing undependable results at other points on the trajectory, es- 
pecially in high-angle or long-range fire. The latter assumes the air 


R= 
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of a standard structure approximating closely to observed data: less 
dense the further one gets above the surface of the earth, in accordance 
with an explicit mathematical law.“ 

The former represented the trajectory by equations which had 
been simplified by a step which assumed (for the purposes of this 
step) that the trajectory is a continuous straight line. This assump- 
tion produced grave errors in high-angle fire. The latter solves the 
equations in their strict form as originally deduced from the laws of 
physics. 

The former approximated the observed atmospheric resistance 
function by certain rather roughly continuous formulas. The latter 
uses the resistance as a tabulated function, as derived from experi- 
ment, and smooth throughout.? 

The former found it very difficult to obtain explicit data without 
assuming? the rigidity of the trajectory, which of course is inadmis- 
sible if the target is considerably above or below the level of the gun. 
The latter is adapted to furnishing just as exact information concern- 
ing any point on the trajectory as concerning the point of fall, and 
hence need not assume the theory of rigidity. 

The former possessed no practical method of correcting the range 
for variable wind aloft, or computing the weighting factors there- 
for, and hence often gave results so incorrect as even in exceptional 
cases to have the wrong sign. The latter can compute, with any 
desired degree of precision, the effect of any atmospheric change 
occurring at any point on the trajectory, on the basis of the usual 
physical assumptions. 

In addition, modern ballistics frankly treats the ballistic coefficient 
as being purely empirical; and also introduces two new corrections, 
which become important in long-range fire, namely corrections for 
rotation of the earth and for changes in the elasticity of the air. 
'The derivation of the present day formulas is, on the whole, as simple 
as the derivation of the formulas of the Ingalls-Siacci system, and 
their use in the construction of ballistic tables and range tables re- 
quires only the most elementary mathematics. 

In view of the foregoing comparison, it is obvious why the later 
method is better adapted than the earlier to the high-angle fire, the 
long-range fire, and the fire at a target considerably above or below 
the gun (including antiaircraft fire), which predominate in modern 
warfare. 


1 Cf. '* Physical Bases" (Ordnance Textbook 972), р. 5. 

3 Cf. ‘‘ Physical Bases" (Ordnance Textbook 972), p. 4. 

з This assumption is: To strike a target at a different level from that of the gun, a trajectory may, without 
changing its form, be rotated vertically, about an axis through the gun. See page 68. 
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The numerical integration which characterizes modern ballistics 
is often called the ‘‘short-arc method." Short arc methods have 
existed in the past. That of Siacci (discarded by him) consisted in 
carrying his approxiniations only during a change of say 5? in 0, 
then starting a new set of approximations based on the value of 6 
&t that point, and carrying this new set for the next change of 5? in 
0, etc. But the words short arcs" and “successive approxima- 
tions" should not lead the student to assume that Siacci's method 
of making a series of approximations over successive short arcs bore 
any close resemblance to modern methods of numerical integration. 

Numerical integration has long been used to compute the orbits 
of heavenly bodies, but was applied to the computation of trajecto- 
ries in this country for the first time in 1917. The necessity for these 
more exact methods was realized in this country as the result of 
reports from England and France in which analogous developments 
had arisen. 

The development of the present ballistic methods in this country 
have been eovered in the introduction. 

Let it not be thought, however, that present methods have rele- 

gated the familiar Ingalls tables to the discard. The Ingalls- 
Siacci methods were devised for guns which were rarely if ever fired 
at an elevation of over 15°; therefore their failure to apply to higher 
elevations is not to their discredit. For fire up to 8°, the French 
still use the old methods, even as a basis for their new tables. Until 
the new American ballistic tables are completed, there is no reason 
why the Ingalls tables should not be used for all low-angle short- 
range computations. 
And even in the computation of higher-angle, longer-range trajec- 
tories by present methods, the Ingalls tables have an important 
place, for by their means the observed-range of range firings can be 
corrected to the range-under-standard-conditions, and the first ap- 
proximate value for C can be obtained from this standard range, 
the muzzle velocity and the angle of departure. 

Then, too, the muzzle velocity is obtained from the observed in- 
strumental velocity by means of Ingalls’ tables, which have many 
other practical uses in ballistic experimentation. 

The object of the new methods is thus seen to be not to supplant 
the Ingalls tables, as being something inaccurate and obsolete; but 
rather to treat these tables as sufficient in the field for which they 
were designed, but as needing to be supplemented in the larger field 
of modern artillery fire, for which they were never intended. 
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QUESTIONS ON CHAPTER VI. 


1. Give the dates and characteristics of the three periods of ballistic development. 

2. What have been the two great obstacles to the development of exterior ballistics? 

3. Define the three branches of ballistics. 

4. Compare the present assumption relative to atmospheric density with that of 
the preceding period. 

5. Compare the treatment of the equations of motion. 

6. Compare the assumptions relative to rigidity. 

7. Compare the assumptions relative to atmospheric resistance. 

8. Compare the treatment of the ballistic coefficient. 

9. Compare the coordinate systems. 

10. What new corrections have been introduced by present methods? 

11. Into what three parts does the study of exterior ballistics naturally divide? 


CHAPTER VII. 
THE MOTION OF A PROJECTILE. 


For the purposes of the computation of trajectories and differential 
corrections, the projectile in flight is treated as a particle. The 
motion may thus be confined to the plane of fire, and the various 
effects of the oblique presentation of the projectile to the air (i. e., 
drift and the range effects of yaw) may be temporarily disregarded. 
Drift is subsequently treated as an empirical deflection correction; 
and the range effect of yaw will probably be treated as a differential 
range correction when the gyroscopic action of the projectile has 

been studied somewhat further than at present. 

Ihe atmospheric retardation of a projectile in flight in still air 
depends upon three things, viz., the velocity of the projectile, its 
physical characteristics, and the density of the air. Then the at- 
mospheric acceleration (i. e., the negative of the retardation) can be 
expressed as follows: | 

26. | а = — “CO 
where F is a tabular empirical function of v. This is the F of the 
Siacci ballistics, and is not to be confused with the Е (formerly 
represented by F) of present ballistic methods. H represents the 
actual density in terms of standard density, and C is an empirical 
constant (different for each projectile), employed to make the 
acceleration correspond to the physical characteristics of the pro- 
jectile. 

By throwing all the effects of a change in velocity into F, then H 
and C can be made independent of velocity. By regarding the 
atmospheric density at the gun as constant, H becomes a function 
of y (i. e., the altitude of the projectile above the level of the gun). 
The exponential function (H = 10-*™5y) has been found to be a 
close enough approximation to the physical facts for all practical 
purposes. С is conceived of as an empirical function of the charac- 
teristics of the projectile, and is constant over any given trajectory. 

For standard density, Н becomes unity. For the so-called stand- 
ard projectile, C is unity.“ Thus F'is the acceleration of a standard 
projectile, traveling through standard air, at velocity v. 

1 This sentence is to be taken as the definition of standard projectile." A standard projectile may be 
of any form, weight, etc., provided only that its C is unity. 
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For convenience (as will later appear), F is replaced by 9G, G 
thus being the ratio of retardation to velocity of a standard projectile 
traveling through standard air at velocity v. For further con- 


venience, the argument of the G tables is 160 


Now it has been found by observation that the true G is practically 
dependent upon velocity alone, regardless of С, Н, and 0. The 
very slight effects of changes in С, Н, and 0 on the true С have 
therefore been taken out of G, have been treated as constant through- 
out any given trajectory, and have been merged in the C of that 
trajectory, leaving the G function a function of v alone. 

Equation 26 becomes: 


97. E FH | GH , 


ee 


E being adopted merely as a convenient expression for of. 


The curve obtained by plotting G against v, shows that this 
function changes most rapidly around the velocity of sound 
(v —330 m/s) and apparently has an inflection at that point. This 
suggests a relation between this function and the velocity of sound, 
which leads to expressing G as vB, where B is a function of the ratio 
between v and the velocity of sound. This relation is not accidental. 
It has a physical basis which can be derived theoretically. 

This gives G the dimensions of velocity. Н has the dimensions 
of density. Accordingly C has the dimensions of sectional density, 
i. e., weight divided by length squared. This suggests representing 


C as 5 where ‘w is the weight of the projectile and d its diameter. 


But to make theory correspond to observation, it is necessary te 
insert an empirical factor of ignorance (i) in the denominator, hence 


wW 
C= i 

Let us now resolve the a of equation 27 into horizontal and vertical 
components, as follows: 

28. E = — ky co — Ex’ 

y"- —Evsin0- — Еу 

But there is an additional impressed acceleration, namely, ‘gravity 

(—g). Ав this acts only vertically, equations 28 become: 


( 2" = — Ky’ 
y” = —Ey'—g 
29. * 2 у 
tan 0 = +; 


з Any other dimensions might be chosen for G, Н, and С, provided that they combine to give correct 
dimensions to Е. See supplement С. 
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These are the equations of motion of a projectile at any point of 
its flight, referred to cartesian axes horizontal and vertical, respec- 
tively, at that point. 

But the earth is not flat. Accordingly, cartesian axes, which are 
horizontal and vertical at one point on the earth’s surface, will not be 
so at any other. Furthermore, gravity decreases with altitude. The 
question therefore arises, how to adapt equations 29 to the actual 
state of the earth. 

The most convenient ways would be as follows: 

(a) The tangent method.” Consider the х axis as horizontal at 
the gun. The equations become (see supplement E): 


| a „ 20 | 
z= — Ex Rt 


ve ГА 2 • 
30. y” = — Ey 9% . e 
Party" 
ГА 
tan G=5 


Where Ё is the radius of the earth, and g, is 9.80 meters per second- 
squared. All ranges, altitudes, and slopes are relative to the axes 
rather than to the earth, and the range at least must be corrected so 
as to relate to the earth, even in the case of short-range fire. But 
for purposes of computation. equations 29 are a sufficient approxi- 
mation to equations 30, except in the case of long-range fire. 

(b) ‘‘The curved method.” Measure z in a circle passing through 
the gun and concentric with the earth. Measure y vertically from 
this circle. The equations become (see supplement E): 


"a E 2 2 U. 
s= E. +. 


2 n 
y” = — E/ — g+ Tm: В 


2 2 ут" 
R 


31. 
v= io aa +. 


L (171... 
tan 0-5, (1 Rt | ) 


All ranges, altitudes, and slopes relate to the earth. For purposes of 
computation, equations 29 are a sufficiently. close approximation 
to equation 31, except in the case of extremely long-range fire. 

(c) Trajectories could also be computed by a third method which 
treats the z axis as tangent to the trajectory at the summit. This 
is sometimes called ‘‘the secant method," for it makes the z line 
through the gun secant to the earth. 
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All range tables now computed are based on trajectories computed 
by "the uncorrected curved method, i. e., z is measured along the 
surface of the earth and y is measured vertically from the surface 
of the earth, but equations 29 are used as an approximation to the 
more precise equations 31. The reason for adopting the curved con- 
vention is expressed as follows in the Aberdeen instructions for 
range firing: 

“The level surface may be taken to be either the curved surface of 
the earth or the tangent plane to the earth at the position of the gun. 
For the purpose of the gunner, the former would be the more con- 
venient in case the levels of gun and target were taken from a 
contour map, and the latter in case the levels were determined by 
sighting from the gun. Аз the difference would be insignificant 
except at long ranges, where presumably the first method would be 
employed, it is the custom at Aberdcen to determine the range for 
the curved surface." 


QUESTIONS ON CHAPTER VII. 


. Why is the projectile treated as a particle? 
. Upon what three things does atmospheric retardation depend? 
. What is the standard projectile"? 
. What are the dimensions of E, H, G, and C? 
. What are the relative advantages of the tangent method" and the curved 
method "? 
6. How are the т and y of a modern trajectory computation considered to be 
measured? 
7. What are the equations of motion used in computing such а trajectory? 


QU d» C9 N m 


CHAPTER УШ. 


COMPUTATION OF TRAJECTORIES. 
{Rectangular method.] 


In modern methods, a standard trajectory is computed from 
-given values of ¢, V, and С, by means of numerical integration (see 
Chap. IV) and successive approximations (see Chap. П). 
. The following tables are used: 
Я а v? 
Table of the G Function. 
Log, H= —0.000045y. 


Some computers considerit quicker to compute log , H by subtracting 
20 from z, pointing off four places, and then subtracting from zero, 


than it is to use the last mentioned table; accordingly the H table 
may be omitted, if desired. 
The following blank forms are used: 


“Trajectory sheet," Form 5042. 
* Computing sheet," Form 5041. 


4, 1“, 7 ½ y, and y” computed as follows: | 

For each of these variables there are, on the trajectory sheet, four 
blank columns, the left-hand column being for the variable itself and 
the other three columns being for the first, second, and third differ- 
ences, respectively. 

The initial data are: The muzzle velocity (V) in meters per second, 
the ballistic coefficient (C), and the angle of projection (9). On the 
first line of the trajectory sheet enter the initial values of z, 2’, t, y, 
and /, as follows: 


х —0 
4 = cos ф 
t =0 
y =0 
y'= V sin o 


Then turn to the small computing sheet and compute the first 
column, for £—0. Ор this sheet, log x” and log y^ should be crossed 
out and log v’ (which is a misprint) changed to log y’. Colog C and 
g are constants throughout the computation. Colog С is obtained 

1 Alternative methods of much merit are given in Supplementa A and G. Ofthe three methods, that 


of Supplement G is the one at present favored by the computers of the Technical Staff. 
| 47 
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from the initial data; gis 9.80. On both sheets F should be changed 
to E, to conform to the approved notation. The logarithms are all 
to base, 10. 

The computing sheet is used to obtain z” and “ from the funda- 
mental equations of motion, derived in the preceding chapter: 


"= Ex 
y'-—Ey'—g 
32. „ GH 
4 € 


The process is as follows: Taking the value p 7 from the trajectory 
sheet, look up at the same time log z' and 100 g in the table of logs 


and squares. Similarly look up E y' and y x 


12 
Add 100 and 160 to get 160 With the latter, enter the G table, 


and take out log G. 

With y, as tabulated on the trajectory sheet, enter the H table 
and take out log H; or compute log H from y without using the table. 

Add log G, log H, and colog C to get log E. | 

Add log z’ and log E to get log Ex’. Add log E and log у’ to get 
log Ey’. 

Look up Ex’ and Ey’ from their logs, in the table of logs and 
squares. Be careful to give Ex’ the same sign as x’ and E/ the same 
sign ав y’. 

Add Ey’ and g algebraically, to get Еу’ +g. 

Change the sign of Ex’ and of E/ +0, and enter them on the tra- 
jectory sheet on the same line as the data on which they were based. 
The reason for this change of sign is that z’’ equals minus Er’ and 
not plus Ex’. Similarly for 4”. 

We are now ready to begin a new line on the trajectory sheet. Our 
second ¢ should usually be 1 second. Enter this value in the Time 
column. Two formulas for integration ahead are available for 
breaking into a new line. This first formula is: 


+i ; 1 5 1 
33. f fdt=i(fit zat py + eredi) 


where t is the time of the last complete line and + is the time interval. 
For example, if f, represents the tabulated value of 2“ at time t, then 
the integral represents the increment of x’ during the interval from 
t to t. 

In proceeding from 1 = 0 to t=}, we have no values of a, b, c, etc., 
corresponding to z’’, and so have only the rough approximation: 


„ 1 
* dt 4 т 
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As this integral represents the increment of z’ from time 0 to time 
4, its value should be entered in the a column corresponding to 2’ 
for time 1. Add algebraically to the value of x’ for time 0, and enter 
the sum as the tentative value of z' for time 1. 

Similarly, integrate ahead for a tentative value of y’ for 1. 

We have now broken into our new line, and are in a position to 
get a tentative value for y for time 1, by means of the standard 
integration formula 


34. JI -( M2 0-12 b.—24 6.40 а.) 

Now return to the small sheet, and compute а new columh, as 
before, using the tentative values of 2’, у’, and y for ¢=}. Enter 
the resulting values of x“ and / on the trajectory sheet on the line 
for t=}, and integrate (this time by formula 34) for z', , and y. 
Substitute these improved values for the tentative values previously : 
set down, and repeat the computation on the small sheet. 

Continue this proces of successive approximations until you 
obtain values of z', 2“, y, у’, and /, which check throughout. Then 
integrate for x by formula 34. The smoothness of the differences of 
z is a valuable check on the accuracy of 2’. 

Then proceed, in the same way, to get the values of 2’, z'', y, y’, 
„ and z, for t=4, t=}, and t=1. Each succeeding line will, of 
course, furnish more of a, b, c, etc., for use in formulas 33 and 34. 

After completing the line for t=1, it may be well, though not 
essential, to extrapolate back for values for a, b, and c for t=0, and 
repeat the computation up tot=1 again. This results in smoothing 
out the curve and giving greater precision to the initial steps. By 
using this smoothing-out process, it is often possible to start with a 
one-second interval? and obtain as precise results as a quarter-second 
interval would give without the smoothing-out process. 

After completing the line for t=14, skip a few lines, copy down 
the values for æ, x’, 2“, y, у’, and y” for times 0, 3, 1, and 14, and 
make up new columns of a and b corresponding to each of these, but 
now based on a half-second interval instead of on a quarter-second 
interval as before. 

Whenever third differences are available, the following formula 
for integrating ahead will be found much simpler than formula 33: 


t+i | 
35. i e een 


It is to be noted that the first difference (at) does not enter into 
this formula, and that the increment obtained is to be added to the 
value for the preceding line. For example, in integrating z” at time 


— 


3 See Supplement G for a method of starting at one-second interval. 
24647—21——4 
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20 to get z’ for time 21, add the increment obtained by formula 35 
to the value of z' for time 19. The increment, if obtained by formula 
33, would be added to the value of z' for time 20. When first using 
formula 35, it would be well to use formula 33 too, and compare 


results, as a check. 
Go on from this point at half-second intervals, remembering that 


4 now equals 3. Continue the process until y becomes negative. 
This means that the projectile has pierced the initial plane and has 
passed below the point of fall. 

It is well to carry z” and y” to hundredths, and z’, х, y’, and y to 
tenths of meters throughout the computation. 

Often during the computation, and certainly at its close, the total 
increment of z, z', , and , from time zero, should be checked by 
integrating z', z”, /, and y” by Simpson’s rule, preferably using 
some type of calculating machine. 

Usually, in computing trajectories, a one-second interval is adopted 
shortly after the start, the change of interval being accomplished 
in а manner similar to that prescribed for changing to a half-second 
interval. т 

Intervals greater than one second should not be used by the stu- 
dent, for the reason that data at one-second intervals are frequently 
necessary for the subsequent computation of the differential cor- 
rections.’ 

To get the terminal values of t, z, z', and у’ (i. e., the values: cor- 
responding to y—0), it was formerly the practice to use formula 6, 
and the interpolation formula of problem 30. The following some- 
what analogous formulas‘ are simpler: 


3 The present practice of the Technical Staff (whether the method of this Chapter or that of Supplement 
G is employed) is to double the interval every time the third differences in the doubled interval indicate 
that fourth differences may be neglected. Accordingly, only very short trajectories are computed through- 
out by a one-second interval, the usual intervals at the end being two or four seconds. The differential 
corrections are seldom required to any intermediate time. When they are so required, it is less work for 
the skilled computer to find them by interpolation than by retaining a smaller interval in the trajectory. 


4 Taylor's expansion of y. A; may be written strictly: 
= At /, At / ,. 


=2 = РР At " AL iv 

Дру A. y. AL fy ANI [v er (yi tt eee J 

In the approximate formulas of the text, all terms in y”, and higher derivatives have been dropped. To 
get a closer approximation, make use of: 


A. 


Similarly: 


Ay" 2 tat 
y", © vs ) 


yht TN Paans y’ 2 27 * 


This gives us: 
a- (GOL V^ 04, 
i 


= (V v Ay" 
Ay - ( E Jat 
and similarly for Ar“ and Ar, the other formulas all remaining the same. 
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36; — У 


Note that the At of this method equals 2А? of formula 6. The 
other deltas are the same in both methods. 

Take the values of z, z’, z'', t, y, у’, and y” from the line in which 
y has the smallest absolute values, and keep a strict account of alge- 
braic signs. Solve for At by successive approximations,’ taking as 
а first approximation: 


When At has been found, the terminal values are obtained as 
follows: 


| Az! =x" „Лі 
Az’ 
= , 3 
37. Az (a tt 2 ум 
Т={+ М 


y's = yc Ay' 


Skip a few lines, and tabulate the values corresponding to t= T. 
Compute the angle of fall by the formula: 


tan w= — Ут ` 
T. 
PROBLEMS. 


(46) Compute the standard trajectory which would result from the 
following initial conditions: 


V =792.47 meters per second. 
С 4.0 
e- 2^ 


5 The work on a slide rule is as follows: 

(а) Set the runner to y. on scale D. 

(b) Bring . on scale С up to the runner, 
(c) Opposite 1 on scale D, read the approximate Ai on scale C, and jot it down for a check. 
(d) Opposite (a of /:) on scale C, read the approximate Ay” on scale D. 


(e) Opposite 14 on scale C, read the approximate Ay’ on scale D. 


Bring +2 on scale С up to the runner. 


(g) Opposite 1 on scale D, read a new approximate А? on scale C, compare with one found above. 

(h) Repeat the process from (d) to (д), until no further shift of the slide is necessary. Algebraic signs 
must always be taken into consideration. Also, be careful to see that the various figures are cor- 
rectly pointed off, so that, for example, a At that should be .043 is not used as .43. 

Then At, Ay”, Ay’, y., (аз a check), Ar“, Az’, and A т can be severally read on scale D without change of 


А : ; , „ Ay” Ay’ , „ Ат” 
slide. They will be respectively opposite 1, (a of “.) /i, у” t=? Vitr? (a of 17 /, *. 2 ,and 


ГА 
X'+ + , 
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Start with quarter-second intervals. Do not use the smoothing-out 
process described on page 49. Change to half-seconds after t=1}, 
and to seconds after computing f = 3. 

(47) Use the same initial conditions and half-second intervals. 
After computing the line for t=2, extrapolate back for a, b, and c 
of z” and y”, for t=0, and recompute to £—2. Change to one- 
second intervals at t=3. Compare the results with problem 46. 


QUESTIONS ON CHAPTER VIII. 


. What two mathematical processes are used in computing a standard trajectory? 
. What tables are used? 

. Briefly, just what is a ‘‘standard trajectory"? 

What initial data are used? 

. What equations of motion are used? 

What are the formulas for integration ahead? 

. When is each used? 

. What is the standard integration formula? 

. How are the computations checked from time to time? 


CHAPTER IX. 
DERIVATION OF AUXILIARY VARIABLES. 


A given ballistic coefficient, muzzle velocity, and angle of departure 
determine a standard trajectory, all the elements of which can be 
computed by the methods laid down in the preceding chapter. In 
that chapter we saw that any given trajectory has, for each value of 
і, a definite value for z, y, 2’, у’, z'', and y". The object of the 
present chapter is to derive three auxiliary variables, , у, and р, 
which shall each have a definite value for each point on each given 
trajectory. In the next chapter we shall see how the formulas for 
the various range corrections can be expressed in terms of these 
auxiliary variables. 

In Chapter III we saw how the effect on z at time T, due to a dis- 
turbance at time £4, could be expressed as in equation 17: 


òX =L or TM Nor! TP by’ 


This is a general expression relating to the motion of any particle 
moving under any definite law and subject to a disturbance at one 
instant of time. Let us now specialize this expression, as follows: 
Let the particle be a projectile moving in accordance with the laws 
of motion evolved in Chapter VII. Let T represent the time of the 
point of fall of the standard undisturbed trajectory. Then X is the 
standard range. Let X, hu, у, and p respectively be L, M, N, and P, 
specialized by these conditions. А is the h of the exponential expres- 
sion for Н (i. e., Н —e-h»), 

^ is thus the change in range due to a unit change in х at time 
ta; hu is the change in range due to a unit change in y at time fa, 
etc. Thus these variable conversion factors are seen to have an 
important physical significance. It is obvious that since z does not 
enter into the differential equations of motion of a projectile, a 
change in z at any point can have no effect other than to shift the 
value of z at all subsequent points by the same amount. Therefore 
^ equals unity, and the expression for a range change may be written: 


38. 6X = zb hu dyty dz’ +p dy’ 


The effect of acceleration changes will not be considered in this 
chapter. 
Strictly speaking, the new range (i. e., X & X) will not be the range 
to the point of fall of the disturbed trajectory (i. e., the point where 
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y=0), but rather will be the range to the point on the disturbed 
trajectory where t= Т. The time of the point of fall on the disturbed 
trajectory will be TA & T, where ôT is the change in Т due to the dis- 
turbance. The z effect, at time Т, of the disturbance at time ta 
will differ from the z effect at time 7+8T by some proportionally 
very small part of the very small quantity 5X, and hence this very 
small part of 8X may be disregarded, by theorem 9. 

In equation 38, 6X is independent of t4, by theorem 18. By 
theorem 15: 


d du 
a, 00 -g) ve 


Therefore, differentiating equation 38 with respect to ta, we get: 
ГА du ГА dv ad ГА г 
39. 0-35 tha by thy by +; be" +v ôr T 2 бу +p бу 


Now, by theorem 14, since z' and y'' are each a function of z', y', 
and y only: 
Ox’ Ог’! dx“ 
„ S—; vM 
io ôx 2 a! + Oy’ Аш 


Let us designate 115 by u', and similarly for »' and р’. 


Substituting the values from equations 40 in equation 39 and col- 
lecting the terms, we get: 


41. 0= -(mw +r +p Borers + 5° a 


ОХ 


Now since бу, 5x’, and ôy’ were chosen independently at time ta 
the algebraic law of vanishing coefficients applies, and each paren- 
thesis above equals zero, whence: 


A bx Òr” 0 
hu =r ay P y 
42. | y= 1-9 p Oz 
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Let us now perform the partiat differentiations indicated in equa- 
tions 42. Since 


g” = — Ех’, and 
y" = — Еу —9, 
then, 
дт’! DE 
or 2 5р7 
d“ OE 
57 oy 
dz” __ DE 
43. , ду 

ду"! DE 
Oz’ Е. y дт” 
oy" „DE 
oy — "Oy 
oy"! OF 


Let us now perform the partial differentiations of E, indicated in 
equations 43. Since 


GH 
Р" | Е = CG and 


v? = 2" +y”, and 


Her, 
then, 
55 085-64 gun Ba (LEE 2) 
Or C or G dv Oz v dv 
oE G oH E T 
OE H òG E dd ov d log G 


бу” 0 dy" 0 du oy ( odo 
Substituting the values from equations 43 and 44 in equations 42, 
we get: 
p= — E(z'v- y'p) 
"= 14 By— Ww (29% 
Да rem 
d log 6 


PS — „75/7 
= — hut Ёр— u'y v dv 
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These equations hold for any possible combination of ба, dy, ёх’, 
and dy’ occurring at any point on the trajectory. 

In these equations, the auxiliary variables are, by their original 
definition, functions of the time of disturbance, £4. The elements 
E, z', y', G, and V are functions of t in the original trajectory com- 
putation; but it is evident that their values, for insertion in equa- 
tions 45, must be the values which they have at the point of disturb- 
ance, i. e., at the point ¢=t,. This consideration becomes important 
in the solution of these equations in the next chapter. 

Equations 45 can be solved for any given trajectory by making use 
of the terminal values of д, у, and p, which are as follows: 


cot o 

Ит = 7р 
46. = 
Pr = 0 


These values аге true because, at the point of fall a unit increase 
in y produces an increase in X equal to cot w, but a unit change in 
z' or у” produces no change in X. 

Auxiliary variables may in a similar manner be derived with re- 
spect to time of flight (T) and maximum ordinate (y,). The three 
sets are then distinguished by the subscripts 1, 2, and 3. , 


bX =), be + hus5y + vdr + рду’ 
ST = Nr + hy + v.52’ + p,0y' 
буа Nor + hu,dy + voz + p,dy’ 


Hereafter in this book (except in supplement D) the omission of 
subscripts signifies that the auxiliary variables relate to range 


changes. 
PROBLEMS. 


(48) What are the formulas for X', u',, , and “, 

(49) What initial, terminal, and constant values are there for ^, 
N, ete. ? 

(50) What are the formulas for x“, m's, »', and p’,? 

(51) What initial, terminal, and constant values are there for \,, 
A^, etc. 1 (Note that here the terminus is the summit.) 


The computation of differential corrections has been simplified by 
the discovery of the following so-called new first integral“: 


4. z' 4- huy! Tr / = 0 
This may be derived as follows. Consider equation 38: 


bz + hy, by +v ôr +p ôy =X 
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Restrict the four independent variables (i. e., óz, dy, dz’, and dy’) 
by the condition that they shall represent a change from point to 
point on the same trajectory, in the same infinitesimal time interval, 
dt. This converts the operator 5 to the operator d, and makes the 
range change equal to zero. Thus: 


dz Nd‘ / + ve + рау’ = 0 
Divide through by dt. Then: 
48. z' - huy! Tv + py’ = 


The similar first integral“ for p equals unity; for 11 it equals 


Zero. 

The above derivation gives the following physical meaning to the 
“first integrals," namely, that an infinitesimal change in the four 
elements (z, y, x’, and /) in time dt, which change amounts to 
shifting from one point to another of the same trajectory, has no 
effect on range or maximum ordinate, and makes a proportional 
change in the time of flight. 

Now substitute — Ex’ for z” and — Ey’ —g for y” in equations 47. 
The first integral" thus becomes: 


z' + huyy' — Е (z'v,- NI) — go, = 0 


or: 
49. | 1 Tu + и, = 9р; 
whence: 
д! + huy’ ds ui 
JN Еа 
Differentiating: 


40 TR T RU“ Tui 
o NN 


Substituting these values in the last one of equations 45, namely, in: 


7 7 ИТЦ d log G 
p = — hu, + Ер,—ц 17 (А 5 | 


we get: 
77 / 77 JA d 
z^ hug" huy + a! = ghm + E( + hay’ +) – 90 1 de J 


whence comes the following equation: 


u^ 2y'u E- hy’ — oy'(“ а Ө] 4172 Ег’, 
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which is the equation used in computing the auxiliary variables. 
This equation produces the value of yu,’ and y, for each value of fa. 
The values of p, and », can then be found from the formulas: 


] 
p, gh Nu m’), 
y’ Hi 
di = ~ Рі Ез’ 


The former of these two is derived from equation 47. The latter 
is obtained from the first of equations 45. 
Thus we have, as the three equations of the system: 


w= 2hEy'm +| Е — hy! — gy’ : e] и. +2 Ех 


1 ГА L 
50. n7, (z' + hy' u, + u^,) 


1 
y, = — р, — үу 
1 1^ Е 


The method of solving these equations is given in Chapter XIV. 


* 


PROBLEM. 


(52) Derive the corresponding equations for u,“, p, and »,; and 
P q 2) 


vs 
Us „ Ps, and уз. 
QUESTIONS ON CHAPTER IX. 


1. Why does the ''first integral" 6 equal unity; and the first integral Y, 


equal zero? 

2. What is the physical interpretation of M, hin, ду, and p,? 

3. What is the physical interpretation of the fact that ^, equals unity? 

4. In this book what does the omission of a subscript from the auxiliary variables 
signify? | 

5. The auxiliary variables are functions of what, on any given trajectory? 

6. Of what are т, у, 27, y’, etc., functions, on any given trajectory? 

7. Isit correct to state that the values of x’ and y^ used in the computation of auxil- 
lary variables are functions of і? Explain. 

8. Upon what subjects of mathematics, discussed earlier in this book, is this chapter 
based? ` 

9. Does X+ôX equal the range of the disturbed trajectory? Explain. 

10. Why are the terminal values of u, », and p as given? 


CHAPTER X. 
RANGE CORRECTION FORMULAS. 
IN GENERAL. 


We have seen in the preceding chapter that at each point (fa) on 
any standard trajectory there is a unique value for each of the 
auxiliary variables: Xr, ш, vi, and p,. 

^, represents the change in range caused by a unit variation in z, 
and is a constant. A,=1. 

hu, represents the change in range caused by a unit variation in у. 

v, represents the change in range caused by a unit variation in z’. 

p, represents the change in range caused by a unit variation in y’. 

As only range changes are considered in this chapter, the subscripts 
will now be omitted, to simplify typesetting. 

- The total range change caused by disturbances may be expressed 
by consolidating all the range effects of equations 38 and 21 into one 
range effect, 6X. 


т T T 
51. 5X or Tu у-у óz' p wis [| die+h | a diy [viia 
te Ute te 


T T | т 
Јар I v uta. | p бу' ta 
te to te 


where the to of any of the four bracketed terms is the time of the 
instantaneous disturbance considered in that term, and the fo of any 
of the six integral terms is the time of commencement of the dis- 
turbance considered in that term, and the f. of no two terms need 
necessarily be the same. | 

Any specialized range change will be represented as A.X , with some 
subseript, so as to conform to range-table usage. 

In the formulas of this chapter the auxiliary variables џи, v, p, and 
their derivatives are functions of the time of disturbance, ta. The 
proper value of x’, у’, E, G, etc. (which are normally functions of t), 
to insert in these formulas are the values for f =t4. 

The formulas given in this book are all (with the exception of 
some in Chapter XI) formulas for the increment in range due to non- 
standard conditions. "Therefore, to correct the observed range to 
standard range, or the map range to range-setter range, the quantity 
obtained from the proper formula must be algebraically subtracted. 
This is in conformity with an agreement between the technical staff 
of the Ordnance Department and the chiefs of Coast Artillery and 
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Field Artillery that range tables shall tabulate the range effects of, 
rather than the range correct ions for, all nonstandard conditions other 
than height of site. 


EFFECT OF NONSTANDARD MUZZLE VELOCITY. 


At the muzzle, 
x = cos G 
y'= V sin 6$ 


A 1 m/s increase in V increases т” by cos &, and increases у’ by 
зіп ф. Therefore, for 59V —1 m/s: 


óz' = соз G 
5y’=sin $ 
Therefore, from equation 51: 
52. AX y = (v cos $+ p sin W140 
PROBLEM. 


(53) Derive the alternative form: 


from the last equation of equations 49. 


EFFECT OF CHANGE IN 5. 


By problem 22, ô cos & equals — sin $ 5¢, and ё sin $ equals cos ф бф. 
It is customary to base the correction on а 1-mil change in G. 


1 
ёф=1 mil 1019 radians. | 
А I mil increase in ¢ decreases z' by V 1019 and increases 


y' by V 1010. Therefore, for 4ф = 1 mil.: 
„ py Sing 
дг = — 1019 
‚_ ү, COS 
бу =} 1019 


Therefore, from equation 51: 


53. АХ» 1019 (p cos $ - v sin ф)+,—0 
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PROBLEM. 


(54) Derive the alternative form: 


_ V р u“ sin $ 
АХ, 71019 cos Fr i io 


CHANGES IN ACCELERATION, IN GENERAL. 


From formula 51, the range effect of variations in acceleration is 
seen to be: 


T ( T 
54. X [vastat f p Sy" dta 
to 


to 


We are interested only in such acceleration changes as are caused 
by variations in Ё, inasmuch as variations in z' or y’ (the other factor 
of x” and /“, respectively) are cared for by the terms in ӧз’ and dy’ in 
formula 51. Thus for the purposes of the present derivation we may 
assume that z' and y' do not vary. Accordingly: 


К. | ôC % bH 
ЭБЗЕ 2 E- 2 N HY Gta 


we ATN 80 8 
by 3E pb--ybE- N Stet aT) 


Substituting from equation 55 in equation 51: 
ôC | 9G | 6H 
bX = l -(-% ta x) E(z'v -- y'p) dta 


Substituting ц’ for — E(z'» t y/o), from equations 45, we get: 


à MY, 
56. ix - (-C) eda 


EFFECT OF NONSTANDARD BALLISTIC COEFFICIENT. 


55. 


For а 10 per cent increase in C alone, formula 56 becomes: 
T Я 
* — (—0.1) u’dt, 
te 


57. AAo = 0.1[ut, — Ar! 


EFFECT OF NON STANDARD WEIGHT OF PROJECTILE. 


An increase in the weight of the projectile over the standard weight 
will increase the ballistic coefficient (C) proportionally, and will de- 
crease the muzzle velocity (V) in accordance with the formula of 
interior ballistics, 

SEL 


— 7 — 
— 


р 
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where p is the standard weight of projectile, and n is one of the vari- 
able coefficients of interior ballistics (approximately —0.3), not to be 
confused with the Gâvre n of the temperature formula (equations 59 
to 64). 

A 1 per cent increase in the weight of the projectile will change the 
velocity by 0.01 nV, which in turn will change the range by that, 
multiplied into the range effect of a 1 m/s increase in V (formula 52). 
A 1 per cent increase in the weight of the projectile will also increase 
C by 1 per cent, which in turn will increase the range by one-tenth 
of the range effect of a 10 per cent increase in C (formula 57). Hence: 


58. AX ,=0.1 AX,+0.01 Vn АХ, 


For values of n, see lesson sheet Р of the Ordnance School of Ap- 
plication: ‘‘Tables for Interior Ballistics," Table A; or use the value 
— 0.3. 

It might seem that n should equal — 0.5 instead of — 0.3; for if 
the muzzle energy produced by a given charge of powder were con- 
stant, regardless of the weight of the projectile, then: 


5 (p+sp)(V+8V)=5 pV? 


whence: 
25V_ _ dp 
V p 
6V 5 
y= 0.5 


But a change in the weight of the projectile causes it to stay a 
different length of time in the bore and alters the friction between 
the gun and projectile, so that the muzzle energy is not constant. 
Empirically, — 0.3 has been found to give more nearly correct results 
than — 0.5. 

In range table firings, when the actual velocity is known, and what 
is desired is the standard range, one range correction is made (by 
formula 52) to correct from actual velocity to standard velocity, and 
one to correct from actual weight to standard weight, using only the 
first term of formula 58. The last term of formula 58 is not used, 
for formula 52 covers the whole matter of velocity. 

If the velocities are taken on special velocity rounds, rather than 
on the range rounds, the velocity of the range rounds is obtained by 
adding 5 to the velocity of the velocity rounds, ё У being found by 
the formula: 

BV _ bp 


V р 
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where 6p is the difference obtained by subtracting the average pro- 
jectile weight of the velocity rounds from that of the range rounds. 

In service firings (or in any case where the range obtained by one 
weight of projectile is known, and what is desired 1з the range to be 
expected of another weight, all other things being unchanged) the 
whole of formula 58 should be used. 


EFFECT OF NONSTANDARD TEMPERATURE. 


A 1 per cent change in absolute temperature alone will now be 
considered, merely from the viewpoint of its effect on the elasticity 
of the air. 

Let s represenf the velocity of sound, s, the standard velocity of 
sound, 7 the absolute temperature, and 7, the standard absolute 
temperature for the altitude in question. A normal temperature 
structure has been adopted by the technical staff of the Ordnance 
Department, based upon the normal temperature which occurs at 
each altitude up to ten thousand meters, when the temperature at 
the ground is standard. This temperature structure is incorporated in 
a table, giving | 

10° 

2то 
with y аз an argument. At the ground 7, equals 518.4? in Fahrenheit 
units, corresponding to 59° Fahrenheit.’ 


Represent G as the product of the velocity by a function of : 


Call this function B (3). Then: 


Let s=s,+6s. Then: 


v 
qawa (s T5) 


Expanding by Taylor's theorem: 


G UB (2 2)-e £2 2 4, higher powers or — 


dB ôs 
dv 80 


59. 6G =G—G,= —v 


It now remains to evaluate 23 апа 08 
| * 80 


1 See “ Physical Bases” (Ordnance Textbook 972), рр. 12, 13. Also note, p. 34, ante. 
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, 


60. G=vB= £ 
4) 


(№. B.—This F is not to be confused with the E of modern ballistic methods, which 
n many of the earlier papers was written F.) 


F. Ат 
G = А 0% 
61. B= And 
Differentiating equation 61: 
dB d (Aan G 
62. um AT DCN 2) z(n—2)À,v*7 = (n—2) t 


Differentiating equation 60: 


dB v dd oe v dG -1)-$ v 21066 -1) 
dv vdv G dv 


Therefore the Gávre n may be defined by the — 


v d log G G 
dv 


n = 


The object of the above differentiations was to get n into the 
form of a logarithmic derivative of G. 


Now to evaluate 5 . In the vicinity of s, s varies as Vr. There- 
0 
fore: 


63. ens 


Substituting from equations 62 and 63 in equation 59: 


0G v? dB 8 
а G 2, 900 -35. 


A бт of 1 per cent would make: 


òG 
eu —2)200 


Substituting in equation 56: 


T 
64. AX, = —0.005 f (n—2)u'dta. 
to 


! An and n are here taken, not as constants over certain intervals of v (cf. Chap. VI), but rather as 
smooth slow-varying functions of v. But they vary so slowly, as compared with v, that they may be 
treated as constants in differentiating. 


hmi. 
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А 6r of 1? Fahrenheit would make: 


5G 
qo 2927. 


Substituting this value in equation 56: 
T" 
65. AX, = -f n 2) s dia: 
te “70 


Instead of using the tabular value of r, in computing AX,, the 
present practice is to use the ground value (518.4?) ; and then treat 
AX, as though the value aloft had been used. This produces a slight 
error, which is probably less than the errors in our knowledge of the 
exact effects of temperature changes. 

Formula 64 is simpler to compute and (as used) is more precise 
than 65, and a percentage change is simpler to use with a ballistic 
temperature; hence formula 64 will generally be preferred. 


EFFECT OF NONSTANDARD DENSITY. 


The weight of 1 cubic meter of air under standard conditions 
(59? Fahrenheit, 750 mm of mercury pressure, and 78 per cent 
saturation) is 1,203.4 gm. The values 59? and 750 mm were chosen 
by the Technical Staff for the reason that these values conform as 
nearly as feasible to the general practice of American, British, French, 
and Italian ballisticians during the World War. Seventy-eight per 
cent was chosen to conform with the present practice of the Coast 
Artillery. Ingalls's Table ITI also uses 59°, but uses 752 mm instead 
of 750 mm, thus giving a density of 1,206 gm per cubic meter. Enter- 
ing this table with 59? and 750 mm (i.e., 29.53'^), one gets 1.002, 
which for practical purposes is indistinguishable from unity. 

The exponential law (Н —e-**) was derived on the assumption 
that E varies as the density. And, of course, density varies as weight 
‘per any unit. Therefore H, E, density and the weight of a cubic 
meter at altitude y are proportional. Whence: 


. ôH _ô(weight of a cubic meter) 
UH. 1203.4 H 


At the ground, H isunity. For an increase of 100 gm in the weight 
of a cubic meter of air at the ground: 


6H 100 


H = 1203400831. 


24647—21—5 
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Substituting in equation 56: 
T 
66. A X« 0.0831 f u dta = 0.0831 (ит — hu). 
te 


Changes in atmospheric density aloft are converted into equivalent 
changes in atmospheric density at the ground, by the formula: 


Change aloft 


67. — H m equivalent ground change. 
For an increase of 1 per cent in atmospheric density. 
6H 
F 0.01. 


Substituting in equation 56: 
68. АХъ= 0. of pdt, = O. 01 (ит — иы). 


It should be noted that this is the same formula as that for АХ с 
(formula 57), except that it is based on 1 per cent instead of 10 per 
cent and is opposite in sign. 

Both the formula for 100 gm increase and the formule for 1 per cent 
increase will be used in the computations of Chapter XVI. 


EFFECT OF A REAR WIND.? 


To determine the effect of a wind moving in the same direction that 
the projectile is moving it is convenient to compare the motion of the 
projectile relative to the air and its motion relative to the earth. In 
other words, it is convenient to regard the motion of the projectile as 
though the projectile were moving along a trajectory relative to the 
air and as though this trajectory were itself moving as an entirety at 
the velocity of the wind. 

If a 1 m/s wind is blowing in the direction of the т axis from time 
to to time 7, the x component of the velocity of the trajectory relative 
;to the air at time f, and thereafter will be 1 m/s less than the т com- 
ponent of the velocity of the trajectory relative to the earth at that 
instant, which amounts to considering the 2“ of the projectile as re- 
ceiving at time to an increment (0z^) of — 1; but the trajectory relative 
to the air is itself moving as an entirety at a rate of 1 m/s from time t, 
totime 7. Thus the effects of à 1 m/s wind are two: (a) А minus 1 m/s 
dx’ at time to; (b) a plus horizontal movement of the entire system 
at arate of 1 m/s for T to seconds. From formula 51, we take the 
terms: ' 


T 
ix- | dbx + viz’ 
te 


* Compare ‘‘ Physical Bases? (Ordnance Textbook 972), pp. 11-12. 
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Now as, by (b) above, dz is changing at the same rate as ta, dôr = dt, 
By (a) above, ёх. = —1. Therefore: 


69. АХ dv -- 7-4-7, 
For & 1 m/s wind blowing throughout the trajectory: 
70. АХ, = Т—», 

d EFFECT OF A VERTICAL WIND. 


Similar reasoning applies to the effects of a 1 m/s vertical wind. 
(1) A ôy’ of minus 1 m/s at time (,; (2) an upward movement of the 
entire system at a rate of 1 m/s for T- t seconds. There are no 
variations in z or z’. 

The effect of the decrease in y’,, is: 


bX =p 5y' = — pr 


The effect of the upward movement of the system is: 
T 
X= f hu diy. 
t o 
But, as ôy is changing at the same rate as ta, déy=dt,; hence: 
T 
X =I. ue. 
Ste 


Thus the tofal range effect of a 1 m/s vertical wind, blowing from 
time t to time T, is: 


T 
71. aK. [ udta — ры. 
te 


For a 1 m/s wind blowing throughout the trajectory: 


D 
72. де 29 f EEN 


CURVATURE OF THE EARTA. 


In the Ingalls-Siacci system of ballistics, the practice was to com- 
pute a trajectory by means of rectangular-coordinate equations of 
motion based upon the assumption of a flat earth, and then apply the 
trajectory to the curved earth by referring it to a rectangular grid, 
whose z axis was tangent to the earth at thegun. Of course, in such а 
system, the farther one gòt from the gun, the farther the surface of 
the earth dropped away from the z axis; thus curvature of the earth 
increased the range by Л Хұ= 0.07848 R? cot о, where R was the 
radius of the earth. | 
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But the present convention is that the trajectory, computed exactly 
as before, is now applied to the curved earth by referring it to the 
orthogonal curvilinear grid discussed in Chapter VII. 

Under the tangent system, curvature was taken out of the range- 
firing range to get the tangent range, all range-table computations 
were based on the tangent range, and then the curvature correction 
had to be added in again. 

Under the curved system, the curvature correction does not have 
to be either taken out or put in. Inasmuch as the system is carved, 
curvature of the earth may be disregarded. 


SITE AND COMPLEMENTARY SITE. 


In the Ingalls-Siacci system of ballistics, the practice was to 
treat a difference in altitude between gun and target, by assuming 
that a trajectory whose chord inclined from the gun to the target 
would have the same shape as though its chord were horizontal. 
This assumption was called the ‘‘theory of rigidity of the trajectory." 
It was defined by Col. Ingalls as assuming ''that the relations existing 
between the elements of the trajectory and the chord representing 
the range are sensibly the same, whether the latter be horizontal or 
inclined to the horizontal." 

It is obvious that, all other things being equal, & projectile will 
go farther if the ground slopes down from the gun, and less far if the 
ground slopes up from the gun, than if the ground were level. 
Assuming rigidity, and an angle of site (i. e., the inclination of a line 
drawn from the gun to the target; positive if up) of e degrees, the 
correction to target range would exactly equal the error caused by 
a change of e mils in 4. 

AX. = [АХ] dpe 


But the theory of rigidity is not sufficiently precise for some 
phases of modern fire (see Gunnery for Heavy Artillery, pp. 66-69; 
Gunnery for Field Service, pp. 23-25). The error in the theory of 
rigidity was formerly corrected for by expressing the effect of this 
error in terms of the change in elevation necessary to overcome it. 
This correction was known as the complementary angle of site," 
the words complementary angle" obviously not being used in their 
geometrical sense. The corrections for curvature, site, and comple- 
mentary site were often combined into one correction. 

The present practice is to derive a single range effect from the 
coordinates of a computed trajectory. The following formulas are 
used: 


. 
tan e = 
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The procedure is to plot є against z—.X for each of the last few 
computed points of the trajectory, draw a smooth curve, and from 
it build up a site-correction table. 

Whether the correction based on these formulas be called ‘‘the 
correction for site" or ''the correction for site and complementary 
site, it covers all of the effects of a difference in altitude of gun 
and target. | 


PROBLEMS. 
(55) Prove that— 
PE fe 
V 86 


See the foregoing derivation of formula 59 for effect of nonstandard 
temperature. Suggestion: by performing the indicated division, 


expand (5) into two terms, one of which shall be 5 By 


Taylor's Theorem for a single variable, expand B G + as) into a 
0 


series of B (2) and its derivatives. 
0 


(56) A shell weighing 15.96 pounds attains a range of 5,000 meters. 
A 10 per cent increase in ballistic coefficient would increase that range 
by 164 meters. А 1 m/s increase in muzzle velocity would increase 
the range by 2.6 meters. What range would be expected of by a 
projectile weighing 15.5 pounds, the Ronee charge and angle of 
departure remaining the same? 

(57) A shell weighing 16.1 pounds is fired at the same elevation 
and in the same gun as that of the preceding problem. The muzzle 
velocity is measured and found to be 20 m/s above standard. What 
range should be expected ? 


QUESTIONS ON CHAPTER X. 


1. What is the difference between AX and 6X? 

2. Give a complete list of the range effects discussed in this chapter, and state as to 
each the unit change on which itis based. 

3. How does the correction for weight of projectile in range-table firings differ from 
the correction used in service firings, and why? 

4, Between what limits should the formulas of this chapter be integrated? 

5. What are the proper values of 2’, y’, E, G, etc., to insert in the formulas of this 
chapter? 

6. How is curvature of the earth corrected for? 

7. Is complementary site included in the correction for site? 

8. Given that AX w, —50, map range=10,000, what range would you set for a vertical 
wind of +2 meters per second? 

9. Would it be possible for AX, to equal —25? Give reasons. 


CHAPTER XI. 
ANGLE OF DEPARTURE CORRECTION FORMULAS, 
EFFECT OF JUMP. 


Vertical jump equals angle of departure minus angle of elevation. 
It is therefore positive if up. The gun is laid by the angle of eleva- 
tion, but it is the angle of departure which controls the shape of the 
trajectory. 

Jump obviously affects the range through the formula for AX. 
The object here will be, however, not to evolve & formula for the 
range effect of jump, but rather a formula for measuring the jump 


N 
JUMP 
J 
AE : 
d cost 
SCREEN 
FIG 6 
itself; inasmuch as it is ¢, rather than X, which is always corrected 


for jump. 

A special blank form is used at Aberdeen Proving Ground for the 
computation of jump. Most of the symbols. used conflict with other 
symbols of ballistics, but are nevertheless here given just as they 
appear on the form. 

S is the horizontal distance, in feet, from the muzzle to the jump 
Screen. 

E is the elevation at which the gun is laid for the jump firing. 

D —S sec Е. 

J is the measured vertical error, in inches, at the jump screen, + 
if up. 

8 is the distance from the trunnion axes to the rauzzle. 


t is the time of flight from muzzle to screen, 1.e., approximately D, 


N is the fall, in inches, due to gravity =} g?=3 (32.16. 12- #)= 


193 2. 
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The pivot of jump has been empirically determmed to be about 
halfway from the trunnions to the muzzle, rather than at the trun- 
nions, jump being due to a combination of rotation about the trunnions 
and other motions, principally recoil. 

7 is the angle of jump. 


dic 1 cos (EJ) 
12 (D+5) 
the reason for the 12 being to reduce (D +5) to inches like (N+). 


But 7 is so small ше сов (E--j) is practically cos E, and 7 (in 
1 


| киды) is practically in L. Sin 1’ is .——. Therefore 


3438 


3438 (N --J) cos Е  573(N4J) | 
12 (D+; 2 = 00р +в) sec E 


EFFECT OF CANT. 


44. 


Although the most marked effect of cant (i. ө., trunnions out of 
level) is the effect on deflection, yet it also has a slight effect on the 
angle of departure. 

See the effect of cant" in Chapter XII. From the development 
there, it is evident that cant has the following effect on angle of 
departure: 


Аф (in mils) = — 1019 cos i (tan $ —tan E) 
Converting into minutes, 
75. Ad = 3438 cos à (tan = tan E) 


The 2 of these formulas is the angle of cant, and has no relation 
to the form factor ?, nor to е? of the jump formula. 

When laying a gun with a quadrant or with a sight shank that 
is capable of being leveled, the indicated elevation is the true eleva- 
tion, and no range correction for cant need be made. Even in the 
case of guns laid by range drum or by nonlevelable scale, the correc- 
tion is apt to be inappreciable unless the trunnions or base ring are 


badly out of level. 
ANGLE OF SITE. 


The point of splash in range firing is of course at a lower level 
than the gun. The exact difference may be ascertained by com- 
paring the height of the trunnions above mean low water and the 
height of the tide. 
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The difference will be so slight that the theory of rigidity" 
(see p. 68, ante) may be assumed without necessitating any com- 
plementary correction. 

This difference, divided by the range (being careful to convert to 
the same units), gives the tangent of the angle of site. (See the 
definition of ‘‘angle of site’’ in Chapter V.) As small angles are 
approximately proportional to their tangents, and as the tangent 
of 1’ is улуу, then the angle of site (in minutes) is equal to the above 
ratio multiplied by 3438. 

This angle must be added to the quadrant angle of departure to 
give the ¢ which would have produced the same X, if firing on the 


level. 
IN GENERAL. 


It is often convenient to know what elevation correction is necessary 
to offset any one of a number of nonstandard conditions affecting 
range. 

Suppose, for instance, that the disturbing cause is a 1 m/s helping 
wind. This is found, by formula 70, to produce a certain AX. 
Find, by formula 53, what Аф will produce AX, numerically equal 
but opposite in sign. This Аф is the desired correction. 

Similarly can be found the elevation correction necessary to offset 
any other disturbing nonstandard condition. 


PROBLEM. 


(58) A gun is laid at an elevation of 2°. From the muzzle to the 
jump screen the horizontal distance is 100 feet. The length of the 
gun from trunnion axis to muzzle is 20 feet. The shell pierces the 
jump screen 30 inches above the point of bore sight. The muzzle 
velocity is 2,000 foot-seconds. What is the jump? | 


QUESTIONS ON CHAPTER XI. 


1. Why is the pivot of jump considered to be halfway between the trunnions and 
the muzzle? 

2. Define vertical jump. 

3. Define cant. 

4. Define angle of site. 

5. How do you find the elevation correction necessary to offset the range effect of 
some disturbing cause? f 

6. In deriving the jump formulas, why is g taken as 32.16, instead of 9.80, as in 
the rest of the book? 


CHAPTER ХП. 


DEFLECTION FORMULAS. 
EFFECT OF CANT. 


Cant occurs when one trunnion is higher than the other. The 
tangent of i, the angle of cant, is obtained by dividing right wheel 
above left" by the ''dis- 
tance between levels." 

Suppose a gun is ele- 
vated E degrees and bore- 
sighted at & point on & 
jump -screen S+s feet 
from its trunnions, and 
then is slowly elevated to 
$ degrees. (See Fig. 7.) SIDE VIEW OF GUN AND SCREEN (TRUNNIONS LEVEL) 
If the trunnions are level, FIG.7 
the line of bore-sight will 
travel vertically up the screen a distance equal to 


(S + в) (tan $—tan E). 
Notice that in the case of cant the pivot is at the trunnions. But this 


А в 
.-RIGHT ABOVE LEFT B «c 


«Тї DISTANCE BETWEEN POINTS 


m 


REAR VIEW OF GUN 
(TRUNNIONS INCLINED) G 
FIG 8 REAR VIEW OF SCREEN 


FIG 9 
vertical line will be tilted 2 degrees to the left, if the trunnions are tilted 
р through i 4 degrees by raising the 


P right trunnion above the left. 
G (See Figs. 8, 9, 10.) The line of 


V8 9+5 -------- ----•  bore-sight will now cut the screen 
GROUND PLAN (S--s) (tan - tan E) sin 7 feet 
FIG 10 to the left of its original position, 


and hence will have a deflection 


whose sine equals | E 
tan - tan -— 
(S + 8) 60) віп ?, i. e., 
sin D= — (tan ¢—tan E) sin 1, 
the minus sign being introduced because, deflection to the right being 
taken as positive, and“ right wheel above left" being taken as posi- 


tive for cant, a positive cant produces a negative deflection. 
73 


74 COURSE IN EXTERIOR BALLISTICS. 


Now, since small angles vary approximately as their sines, and 
since the sine of 1 mil is urg, 
D (in mils) =- $22 =1019 sin D 


46. D = —1019 (tan - tan E) sin 2. 


EFFECT: ОЁ LATERAL JUMP. 


Let I be the lateral jump, in inches, measured on the jump screen, 
plus if to the right. Let: be the deflection angle of side jump. Do 
not confuse this ; with the form factor 7, nor with the i of the cant 
formula. The other symbols are the same as those used in the dis- 
cussion of jump in Chapter XI. Notice that here, as in the case of 
vertical jump, the pivot is taken as half way between the trunnions 
and the muzzle. 

The horizontal distance, in feet, from the pivot to the screen is 
845 sec E. But 5 is so small compared with S, and sec E is so 


nearly unity, that the expresssion may be written S +5 without 
material error. 


07 I 
OE ир e 
12(8+5) 
in mins.) 518 t 34381 5731 
“ smil’ ( 5) 28 +8 
12 (8+5 
2 
7 : 170 7 
71. i (in mil)-5 gF; 


. Lateral jump is usually not computed, but is included in the 
„drift“ as tabulated in range tables. 


EFFECT OF CROSS WIND. 


A 1 m/s cross wind is treated much the same as a 1 m/s range wind 
or a 1 m/s vertical wind (see Chap. X? А cross wind is positive, if 
blowing from left to right. 

А +1 m/s cross wind, blowing from time to to time Т, causes the 
trajectory relative to the air to bend abruptly to the left at time to, 


by an angle whose tangent is 2, Therefore the deflection in meters 
at time T due to this cause alone will be to r- ti, as 1 to * 1. e., 
rai) to the left. | 


te 
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But the trajectory relative to the air will itself be moving to the 
right at arate of 1 m/s for (T- %) seconds, or a distance of (T — to) 
meters. Therefore the net deflection, in meters, is: 


78. D. F Ed 


the values of z and x’ being those at time tẹ. 
In mils, assuming the wind to blow throughout the trajectory, 


. 19. D, = 1,019 (X 25 
їо 


DRIFT. 


Drift is due to the gyroscopic action of the projectile. 

Drift is at present treated as wholly empirical. In range firing 
the total angular deviation from the line of fire is measured (posi- 
tive if to the right, negative if to the left). From this is subtracted 
algebraically the effects of cant and cross wind. The remainder is 
the drift (plus lateral jump, of course, which, for range-table pur- 


poses, is usually included in the tabulated value of drift). 


PROBLEM. 


(59) In the firing of problem 58, the right trunnion was 1 inch 
above the left, the distance between the points where the levels were 
taken being 3 feet. The shell pierced the jump-screen 9 inches to 
the left of the point of bore-sight. The gun is then elevated to 10° 
and attains a range of 6,421 meters. There is а 9.5 m/s cross wind 
blowing from left to right. The time of flight 1з 17.1 seconds. The 
shell lands 40 meters right. What is the drift in mils? For the 
purposes of this problem do not include the lateral jump in the 
drift. 


QUESTIONS ON CHAPTER XII. 


1. If deflection is given in meters, how should it be converted to mils? 

2. Why, in the case of cant, is the pivot considered to be at the trunnions? 
3. What is the cause of drift? 

4. How is drift determined? 

5. What does the range-table value of drift include, besides drift proper? 


CHAPTER XIII. 
ROTATION OF THE EARTH. 
IN GENERAL. 


The rotation of the earth has two effects on a projectile in flight: 
(a) The higher the projectile goes, the more must its velocity bo 
altered, in order to maintain the same linear velocity relative to the 


FIG il 


earth; (b) centrifugal force offsets to some extent the attraction of 
gravity. | 

One of the assumptions on which the computation of & standard. 
trajectory is based is (see Chap. V): 2. The earth is motionless. 
(The average effect of the rotation of the earth on gravity is included. 
in the assumed value of g.)”’ 

The object of the present chapter will be to ascertain the difference 
between equations of motion based on the assumption of a motion- 
less earth and those based on the assumption of a rotating earth, 
to separate out from this difference such effects as the standard 

1 For a nonmathematical treatment of this subject from another viewpoint, see Physical bases“ 


(Ordnance Textbook 972), pp. 9-11. 
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trajectory includes in g,, and to base a rotation correction on the 
remainder. 

Inasmuch as most discussions of the rotation of the earth adopt 
the convention of rectangular axes (the z axis being tangent to the 
earth at the gun), and as the corrections thus derived can be shown 
to be equally applicable to the convention of a curved grid (with z 
lines concentric and y lines radial), this chapter will develop the cor- 
rections on the rectangular basis. 

Let us consider a projectile fired from a point whose latitude is J, 
in a direction whose azimuth is a, measured clockwise from the south. 

Let us consider two sets of rectangular axes, coincident at the 
instant the gun is fired. The origin is on the equator at the same 
longitude and altitude above sea level as that of the gun; the y 
axis is vertical upward; the z axis is horizontal to the east; the 2 
axis is horizontal to the south. Inasmuch as the orbital motion of 
the earth has a negligible effect upon the trajectory, we shall con- 
ceive of the axis of rotation of the earth as fixed in space. 

Now consider that one set of axes (designated by m) moves with the 
earth; and that the other set (designated by f), although coincident 
with the moving set when the projectile leaves the gun, thereafter 
remains fixed in space, and does not rotate with the earth. 

Let ¢ be the time which has elapsed since the gun was fired. Let 
R be the radius of the earth. Do not confuse this E with the symbol 
for retardation. Let Q be the angular velocity of rotation of the 
earth.? | 

At time ta, the angle between the two sets of axes is Ota; so, by the 
familiar formulas of analytic geometry: 


m = Zt COS Qta — (ус R) sin Qta 
80. 5 sin Qta + (yr + E) cos Qta 


Zm = Zt 


The subscript m means referred to the moving set of axes. The 
subscript f means referred to the fixed set of axes. 
Differentiating equations 80, we get: 


81. 3] m 'i sin Qta + V'rt cos Qt, ＋ Tm 


fom = cos Qa tr sin Qta -m - RQ 
2 mm — Eft 


The meaning of the double subscript is as follows: For instance, 
z' means the time derivative of z;, which derivative is left referred 
to the fixed system of axes. If we were to transform 2777, by formulas 
analogous to equations 80, so as to refer to the moving system, we 
should get Z' em: 

3 Since there are 86,164 mean solar seconds B 8 sidereal day, 


Q= 50,164 00007292 radians per second. 
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Differentiating equation 81, we get: 


х! ЕТ 2271er cos Qt4— y" et sin Qta — 29' mm + mD 
V''mmm = T tf sin Qa T ttt COS Ota + 2 mm2+ (ym + R) O 


77 77 
2 mmm fff 


82. 


The meaning of the triple subscript is as follows: For instance, 
z'. means the time derivative of 2777, which derivative is left re- 
ferred to the fixed system of axes. If we now transform 2’ fr by 
formulas analogous to equations 80, so as to refer to the moving sys- 
tem of axes, we shall get х''г„. These formulas are: 


1 ttm = tt C08 Ma-) rer Sin Qta 
V''ttm = ege Sin Qa Ty ere cos Qta 
77 ah 

2 ftm =2 tft 


83. 


Substituting from equations 83 in equation 82, we get: 


2 md = tm 2) mm + tm? 
84. 4j ^ nimii =Y ttm + 2z! e + Ym? + RO? 


77 Т; 
2 mmm=2Z ffm 


The symbols with the subscripts mmm are the components of rela- 
tive acceleration; those with the subscripts ffm are the components 
of absolute acceleration; both being referred to the moving set of 
axes, which set is that used in computing trajectories. 

Equations 84 can now be simplified by dropping the negligible 
terms Im% and y,,0?. Even with a gun firing 100 miles, neither of 
these terms could exceed 0.00001 on апу part of the trajectory. 
Thus: 

85. Y rel = be + 2 rei Ro? 


1 rel 4s 27 rel 
z’ “Fel =2’ ‘abe 


Now, regardless of whether we consider the earth as moving or as 
motionless, what we are interested in is the motion of the projectile 
relative to the earth. If we assume that the earth is motionless, all 
terms involving Q in equation 85 would drop out, and we should have: 

1 rel & abe 

86. rel == Y abe 

2“ rel 2 abe 

If we consider the earth motionless, except for its centripetal 
acceleration, the result is to add 20° to the right side of the ‘second 
equation of 86. Thus: 

87. rel nd MP + RO 


1 rel = Z' abs 
2 rel 22 “ate 
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Equations 87 represent relative accelerations under the standard 
assumptions of trajectory computation. Equations 84 represent 
relative accelerations as they actually exist. Therefore the effect 
of the rotation of the earth may be found by subtracting equations 
87 from equations 85. 


óz'' = —2y'Q 
88. y” = + 22'0 
ó2/' = 0) 


Let us now shift the axes until they become the axes used in tra- 
jectory computation. First rotate the system northward about the 
center of the earth through an angle of І degrees. The origin will 
then coincide with the gun. Equations 88 become: 


89. by"! = +20(z’ cos l) 


wr- — 20(y' cos 1+2’ sin 1) 
62'' = 4-20 (z' sin l) 


Let us now rotate the т and 2 axes clockwise around the y axis, 
through an angle of (a—270) degrees. The z axis will now point 
in the direction of the line of fire; z^ will equal zero. Equations 89 
therefore become: 


90. by’’ = —20z' cos L sin a 


b= +2Qy’ cos l sin a 
oe“ = + 2Q(y’ cos l cos a+2’ sin D) 


Now equations 90 were derived on the assumption of a rectangular 
system of coordinates, whose z axis is tangent to the earth at the 
gun (the ‘‘tangent method’’) instead of on the assumption of a sys- 
tem whose abscissas are measured along a circle concentric with the 
earth and whose ordinates are measured radially from that circle 
(the ‘‘curved method"). 

It now remains to be shown that equations 90 apply equally to 
the curved method. For this purpose equations 90 must be derived 
So as to relate to instantaneous axes respectively vertical and hori- 
zontal at the point where the projectile is at that instant,’ instead 
of axes respectively vertical and horizontal at the gun. 

In other words, any arbitrary point on the trajectory will be 
taken, and equations 90 will be derived with respect to motion at 
that point. Then, as that point was any point on the trajectory, 
equations 90 will apply to all points on the trajectory. That point 
will be called the instantaneous point. 

Now revert to page 77, and consider (in place of the two sets of 
axes there mentioned) two sets coincident when the projectile reaches 


3 That is, the same axes as those of equations 111 in Supplement E. 
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the instantaneous point. The origin is on the Equator at the same 
longitude and same altitude above sea level as the instantaneous 
point. Except as otherwise stated, the development is the same as 
the preceding. | 

ta is the time which has elapsed since the projectile was at the 
instantaneous point. E is the distance from the center of the earth 
to the instantaneous point. 

Equations 80 and 84 are developed for a point on the trajectory 
very near to the instantaneous point. This near point is then allowed 
to approach the instantaneous point as a limit, with the result that 
the terms z,0? and /“ become zero and equations 84 become 
equations 85. 

Equations 86 and 90 evolve as before,* with no change in their 
derivation, except that AO now represents the centripetal accelera- 
tion at the instantaneous point instead of at the surface of the earth. 
But as gravity is assumed to be constant, regardless of altitude 
(which assumption can be shown to cause ap inappreciable error even 
with the most powerful guns), RQ? can still be considered to be in- 
cluded in g. 

Equations 90 are thus seen to represent the increments of acceler- 
ation due to rotation, and to be correct for either the tangent or 
the curved method. In the tangent method they refer to a rec- 
tangular coordinate system whose y axis is vertical at the gun. In 
the curved method they refer to the directions which are horizontal 
and vertical at the particular point on the trajectory. 


RANGE EFFECT. 


The range effect of rotation of the earth may be found by insert- 
ing the values of ôr” and ôy” from equations 90 into equation 54 
of Chapter X, with the following result: 


T 
91. Aa = 20 соз {зїп а | (vy’—pr’)dts 
Е 0 
20 equals 0.000 14584. 


DEFLECTION EFFECT. 


Consider now the deflection effect (D in meters) of an increment 
in lateral velocity (6z’) occurring at time £4. Then: 


D _ 62’ 

X х! 
Whence: 

D- be’. 


4In the curved method, equations 90 relate to instantaneous cartesian axes which are horizontal and 
vertical at the instantaneous location of the projectile. To transform them, зо as to relate to the curvilinear 
grid of the curved method, we might use equations 120 of supplement E. But ifthe increments be regarded 
as analogous to infinitesimals of the first order, the only effect of such a transformation will be to add infini- 
tesimals of a higher order, containing Ё, which terms can hence be disregarded. Hence equations 90 may 
be considered as remaining unchanged. 
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By differentiation, followed by integration, as in the derivation of 
equation 21, we get the following expression for the deflection effect 
of an increment in lateral acceleration, occurring throughout the tra- 
jectory: 

92. р- [`+ 62'' dta. 

0 


* 


Whence, substituting the value of 52'' from equations 90: 


rA 4l" 
93. Da=2Q cos I cos a [ y'dt4 +290 sinl E (X-) dta. 
0 0 


x! 


For convenience, the following rotation coefficients“ have been 
adopted: 


T 
A = 0.00014584 f (x’p—y’v) dta. 
0 


T 
94. B =0.00014584 f (X —2) dta. 
0 
T — 
0 —0.00014584 f = - 2 y'dta. 
0 


Note that B is not the B function of atmospheric resistance, and 
that Сіз not the ballistic coefficient. 

Therefore the range and deflection effects of rotation of the earth 
become, from equations 91 and 93: 


AX g — А cos { cos a. 


95. | Ро ВЫ О cos bee a 


A close approximation to the three rotation coefficients is the value 
which would obtain in a vacuum,‘ namely: 
A=QX T (cot 5—5 tan ¢). 
В=ОХ Т. 
0-5 x T tan e. 


PROBLEMS. 


(60) Prove that tm% < 0.00001, as stated on page 78. 
(61) Derive equations 89 from equations 88. 
(62) Derive equations 90 from equations 89. 


* See Physical Bases (Ordnance Textbook 972), p. 11. 
24647—21—6 
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QUESTIONS ON CHAPTER XIII. 


1. What are the two effecta of rotation of the earth? 

2. For which of these need the range be corrected? 

3. Do the formulas derived in this chapter apply to the “tangent method” or to 
the ‘‘curved method” of trajectory computation? 

4. What acceleration is used in the computation of a standard trajectory? 

5. What acceleration should be used to account for all the effects of rotation of the 
earth? 

6. Explain the ratio given under “ Deflection Effect." 

7. Explain the figure 0.00014584 in equations 94, 

8. Firing in vacuo, what value of e would cause the range effect of rotation"offthe 
earth to be zero? 


CHAPTER XIV. | 
COMPUTATION OF DIFFERENTIAL CORRECTIONS. 


The computation of differential corrections is based upon the fol- 
lowing equations, derived in Chapter IX. 


Hl , / d d log G G , ыў 
u” == 2Һ Еу u+| Е-м'—ду ie: bdo )| +2 Ex 
wmf u“ dta 


te 
96. и = ) иа А 


29 1 (sl TM и”) 
/ 
dn ў -Er 

The method of procedure is to integrate 2“ ahead to get an approxi- 
mate value of р’, integrate the approximate ш’ to get an approximate 
value u, and then get an approximate value of “ by substituting the 
approximate values of и and “ in the first above formula. Then 
integrate u“ to get р’, and p’ to get и; continuing the process of suc- 
cessive approximations for any one line until the values check. 
Thus we employ numerical integration and successive approxima- 
tions in а manner very similar to the method of computing a standard 
trajectory. 

The data is taken from a trajectory computed by the methods of 
Chapter VIII. The computing form is made by taking a sheet of 
paper about the size of the trajectory sheet of the trajectory compu- 
tations, having horizontal blue lines one-fourth inch apart. Rule a 
vertical line 14 inches from the left-hand margin, then one one-half 
inch from the first, then one one-half inch from that, then every 
inch or three-fourths inch across the page. 


Two special tables are used, being entered with 100 as an argument. 


ГА 


These tables are: 


Table of f (v) = 10° (A +9 ©) 


Table of (n—2) 
88 
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On the computing sheet, the times run across the page instead of 
down it, and run from T to 0, instead of from 0 to T. 
The first four columns are filled in as follows: 


u (approx.) 
p’ (approx.) )): 


10° cr vd» )--- 


(2)=0.01 (ur — t) 
Dens. —(2)X8.310.... 
F 
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7—1, — ꝛ 


W,-T-—t—»....... I 
(14) =0.001 B..----.- 
C 


—— — | ann | e —ÀÀ— M —BGRü: 


> оъ 20 = о ә 
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The fourth column above (“number of decimals”) need not be 
copied. It is given above merely to indicate the correct number of 
decimals to which to carry the computations. 

Put the total time of flight ( T) at the head of the fourth column, 
and at the head of the succeeding columns put all preceding whole 
values of t, from Т back to zero. 

Fill in lines 4, 5, 6, 7, 42, 50, 26, 15, 11, and 31 with data from the 
trajectory computations, and lines 60 and 9 with data from the 
special tables. It is optional whether or not to fill in on the trajec- 


4 
tory computation sheets themselves the values of E, X — z, L, etc. 


Line 8 is computed from line 7, using + 0.0001036 as the value of 
h. Line 10 is computed from lines 5 and 9. Line 12, from lines 
10 and 11. 

The numbers in the column labeled “Next” show which line to 
proceed to, after completing the line in which the number occurs. 
In the absence of any number in the “ Next" column, proceed to 
the next line below. These numbers have no relation to the italicized 
numbers in parenthesis in the left-hand column. 

The work, up to line 12, has been line by line across the page. 
But from now on (until line 23 is completed in all columns), we pro- 
ceed to solve each column by numerical integration and successive 
approximations. 

The procedure is to start with values for u and u“, obtained from 
the following formulas. For time 7, the following are precise: 


1/0. 


For the next two columns the following are approximate: 


" р 27 
. N (Tt) V^ 
ra M (T. 


u’ 
4 Af. 3 (T—1). 


The values of » are useful merely for a check. 
Thereafter the approximate values of u and д” (to be entered in 
lines 2 and 3) may be found by integrating “ ahead by formula 35: 


tH a 
[fao e Oed. o ) 
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to get the increment of ; and then integrating u' by formula 34: 
А 1 1 1 1 
Säi (fizun) 


to get the increment of и. и is positive throughout; “ is negative 
throughout. j 

Although we are integrating from time T toward time zero, the 
student should proceed exactly as in the method laid down in Chap- 
ter VIII. Wherever it is necessary to change signs in integrating, 


T 
that fact is indicated by a minus before the f in the left-hand column 
| te 


in lines 19 and 22. The first, second, etc., differences will be formed, 
proceeding from Т toward zero'; but the time interval (i) will be 
treated as positive. For a mathematical explanation of all this, see 
supplement B. 

Note that the time interval between the first and second column 
is different from the uniform time interval of succeeding columns. 
The first difference obtained by differencing the first and second 
columns is used only in integrating for values for the second column. 
There should be no second difference in the third column, no third 
difference in the fourth column, etc. 

The procedure, therefore, is as follows: Insert the precise values 
for time T in lines 2, 3, 19, and 22 of the Т column. To get the first 
precise value of џ, it will be necessary to use straight multipli- 
cation and division, preferably on a calculating machine, as 
logarithms will not give results precise to enough places. Insert the 
approximate values in lines 2 and 3 of the next two time columns. 
Use these values in lines 13 and 14, and thus get three values of n“ 
in line 16. Form the first differences of д”/, and integrate for the 
increment of џи’, using formula 34 above. Similarly integrate for 
the increment of и. The values thus obtained will usually check 
so closely with the approximate values, that u'' will not have to 
be recalculated. 

Start each succeeding column by integrating р’ ahead to get 
the approximate increment of u’. Add this algebraically to the last 
precise value of “, to get the new approximate “.. Integrate this 
approximate џи’ to get the approximate u. Use these values in lines 
13, 14, and 16, and proceed as before. 

After completing the fourth or fifth column of computations, it may 
be well, though not essential, to extrapolate back for values of u“, p’, 
and their first, second, third, etc., differences for a time equal to the 
next whole number larger than 7', i. e., one whole second larger than 
the time recorded at the head of the second column of computations. 


1 That is, a, equais/,—f;. See the footnote on p. 30. 
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Starting from this point we can now perform our integrations, with 
a constant time interval and a full set of a, b, c, etc., for use in the 
integration. This results in smoothing out the initial steps and in 
giving them a greater precision. 

When р’ and и have been computed clear across the page, we are in 
а position to compute the differential corrections. АП of the following 
work is done a line at & time clear across the page. 

First comes the effect of an increase of 100 gm in the weight of a 
cubic meter of air (using lines 24 and 25); from formula 66: 


AX = 0.0831 (ur — ut). 


The value for time zero is the range-table value, the other values being 
used as & basis for weighting-factor curves, as will be explained later 
(see Chap. XV). 

Then comes (in lines 26-34) the computation of the auxiliary 
variables p and ». 

Then come (in lines 35—49) the three rotation coefficients, formula 94: 


А —0.00014584 (^ (%-) dta 
B —0.00014584 f (X-) dta 
C= 0.00014584 f 7 (45%) y’ dta 


which enter into the formulas for the range and deflection effects of 
the rotation of the earth. A, B, and C are respectively obtained by 
multiplying by 1.4584 the values, for time zero, in lines 40, 44, and 48. 

Then come (in lines 50—59) the effect on deflection of a cross wind, 
and the effect on range of a range wind or a vertical wind, 1 m/s 
being taken as the unit wind; formulas 78, 69, and 71: 


X—z 
D, -T -[. 
AX wx - T—t— to 
AX wy = — рі, +h pdt, 
to 


The value for time zero is the range-table value, the other values 
beng used as a basis for weighting-factor curves, as will be explained 
later. (See Chap. XV). 

. Then comes (lines 60—66) the range effect of the change in elasticity 
due to & 1 per cent increase in temperature; formula 64: 


T 
AX,=—0.005 f  (n—2) шд 
to 
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The value for time zero is the range-table value, the other values 
being used as a basis for weighting-factor curves, as will be explained 
later. (See Chap. XV.) 

This completes the calculations of the columns. We now proceed 
to certain formulas which are based merely upon ‘initial or terminal 
data; formulas 53, 52 and 57: 


y | | 
AX ,— 1919 p cos ó—» sin 2 


АХ „= [» cos G +p sin ф] o 
AX,=0.1 (- ur) 


These are, respectively, the effects of a 1-mil increase iv ¢, a 1 m/s 
increase in V, and a 10 per cent increase in C. 
They can be checked by the following (see problems 54 and 53): 
u'sin ф 


[р cos G- sin ф] -U TE Ex’ 


&=0 


A cos 9 
АХ; = | Ex’ to 0 


[4 
EP for time zero being taken from line 33. 
The computation sheet is completed by the following summary of 
effecta: 
1 wm /s change in V. 
1 mil change in ¢. 
- Rear wind, 1 m/s. 
Vertical wind, 1 m/s. 
Cross wind, 1 m/s. 
100 gm/m increase in density. 
10 per cent increase in C. 
1 per cent increase in absolute temperature. 
Rotation A. 
Rotation B. 


Rotation C. 
PROBLEM. 


(63) Compute the differential corrections complete from the fol- 
lowing data: V=579.1 m/s; C=2.3; ф:= 5°; T=8.824 sec. 
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t 0 1 4 
FFF 576. 9 520. 2 401.2 
V 50. 5 36.2 1. 8 
Е.... 0. 1075 0. 0994 0. 0732 

D 
ium 3354 2719 1609 
2. х Ес 0 548 1913 

t 5 6 7 8 8.824 
FCC 374. 7 353. 3 336.2 322. 4 312. 9 
7 —7.8 —16.9 —25.7 — 34. 2 — 41. 1 
РОИ 0. 0636 0. 0540 0. 0453 0. 0388 0. 0344 

v? Зет 
1000 . 1404 1251 1137 1051 996. 1 
FC 2301 2664 3009 3338 3599 


QUESTIONS ON CHAPTER XIV. 


1. Why are the times, tabulated in the first row of the computing form, designated 
by t, instead of ta or t? Why by t in problem 63? 

2. Why the minus sign in the terminal value of д? Compare equations 46 of Chap- 
ter IX. 

3. Why is no first difference formed between the first two columns of computa- 
tions? 

4. How is the first part of the computation smoothed out? 

5. The values from what columns are the range-table values? 

6. For what are the values from the other columns used? 

7. What difference in methods or algebraic signs, from those of Chapter VIII, is 
used to compensate for the fact that in the present chapter we are integrating from T 
to earlier times? 

8. Suppose a tabular function for times 4 seconds, 3 seconds, and 2 seconds, i. e, 
fof, and fa. What is the value of a:? What would this have been in the notation of 
Chapter VIII? 


CHAPTER ХУ. 
WEIGHTING FACTORS. 


In considering how weighting factors are derived from the compu.» 
tations described in the preceding chapter, let us take as an example 
the weighting factors for range wind. 

The computation sheet furnishes, among other data, the effect of & 
1 m/s wind ор the range. This is entered in the range table, so that 
a battery commander, by multiplying this by the number of meters 
per second in the range component of the wind, can calculate the 
wind correction to apply to his map range. 

But this step by the battery commander is based on the assump- 
tion of a uniform wind at all altitudes, which is a condition that sel- 
dom, if ever, exists. Тһе wind constantly changes both its velocity 
and direction from one altitude to another. Velocity usually in- 
creases with altitude. Accordingly, some means must be devised for 
calculating the velocity and azimuth of & purely fictitious uniform 


wind which would have an effect on range equal to the combined 


effect of all the actual winds met by the projectile in its flight. "This 
fictitious wind is used by the battery commander just as though it 
were the actual wind. It is called the “ballistic wind." The Brit- 
ish call it the equivalent uniform wind," which is a very apt name. 
In the field, the meteorological service measures the average direc- 
tion and velocity of the actual wind for successive strata of 250 


meters each. 


Our problem is now to determine, for any given trajectory, what 
weight is to be given to the wind of each stratum in making up the 
ballistic wind. This is done by determining what proportional part 
each stratum plays in making up the total correction for a 1 m/s wind 
on the computation sheet of the preceding chapter. The procedure 


is as follows: 


From the trajectory sheet of the original trajectory computations, 
find the time corresponding to y’ = O, and the maximum ordinate (i. e., 
the value of y corresponding to this time,) by methods analogous to 
those used in getting the terminal values: 

ГА 


. = 5 4 A (At +1) (At+2) | 
97. 31 4ſſ 
на мш Ses a 


1 The Technical Staff practice is to use the method for slide rule given in footnote, Chapter VIII, for 


finding terminal conditions, adapting it to the present conditions. 
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taking the values of y' and its a, b, с, etc., and of y and its a, b, с, etc., 

from the line whose y’ is nearest to zero. In the first equation, use 

the a, b, c, etc., of y' and solve by successive approximations, taking 
— 477 

for a first approximation At= =”. In the second equation, use 


the a, b, c, etc., of y. The result is y,, the maximum ordinate. 
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FIG 12 
Against t tabulate y from the trajectory sheet and T—t,—», from 
line 53 of the differential computation sheet, and calculate " and 


Leue, as follows: 
T—» 


This tabulation gives for each value of to: The altitude of the pro- 
jectile; the range effect of a 1 m/s wind blowing from time ¢, to time 
T; the ratio of the altitude to the maximum ordinate; and the ratio 
of the aforementioned wind effect to that of & 1 m/s wind blowing 
throughout the entire flight of the projectile. 
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T—t—v 
Now plot, on a small sheet of coordinate paper, —R— against 

0 
ET and connect the points by a smooth curve. (See fig. 12.) This 
curve will run from the point (1, 0), representing time zero, to the 
origin, representing time 7; and will be tangent at its summit to the 


line J =1. 
8 ; 
Let AB be drawn parallel to the axis of abscissas, a distance of any 
5 above it. Then EB represents the proportional effect of a 1 m/s 


wind blowing from the time the projectile first reaches altitude y, 
until the point of fall. EA represents the proportional effect of a 1 m/s 
wind blowing from the time ,, 
the projectile reaches altitude : — iA 
y in its descent until the point ГЛ 
of fall. Therefore AB repre- 8 B 
sents the proportional effect „ 
of & 1 m/s wind blowing above E 
the altitude y. | * ENN 

On a second sheet of coordi- 5 Pinaceae 
nate раат plos 48 (mews «ECE 
ured to the left from the line ; CCC 


T 
2 Wokomatesingonedeny NL 


of unit abscissas) against A 


(See fig. 13.) Then FB equals E task ккк het T 
1— АВ, equals the propor- 
tional effect of a 1 m/s wind о . 2 3 4 5 & 7 8$ 3 10 
blowing below altitude y. | | 

The proportional effect of FIG 13 


a wind blowing throughout any stratum is the difference between 
the FB corresponding to the top y of the stratum, and the FB cor- 
responding to the bottom y of the stratum. 

Thus, to get the weighting factors for the trajectory in question, 


mark off on the curve like that of figure 13, points whose А are respec- 


к 0 250 500 750 
tively —» ——» —» — 

Us Ys Ys Ys 
first and second points is the weighting factor for the first stratum; 


the difference in abscissas between the second and third points is the 
weighting factor for the second stratum, etc. As a check, the sum 
of all the weighting factors should be unity. 

Of course, it would be possible to calculate the weighting factors 
directly from the curve like that of figure 12, but it is more convenient 
to have the weighting-factor curves in the form of figure 13, as it is 


» etc. The differencein abscissas between the 
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usual to plot on the same sheet the curves for various angles of 
departure of the same gun. 

It should be noted that, for any given trajectory, the weighting 
factors for making up a ballistic wind to use for range corrections are 
quite different from the weighting factors for making up a ballistic 
wind to use in deflection corrections. Therefore, according to present 
methods, there are for each combination of C, V, and Фф, two ballistic 
winds: (a) The range ballistic wind, whose range component is used 
as a basis for range corrections; and (b) the cross ballistic wind, 


whose lateral component is used as a basis for deflection corrections. . 


The foregoing method for getting weighting factors for range bal- 
listic wind is equally applicable to cross ballistic wind and ballistic 
density, the last requiring slight modifications of a fairly obvious 
sort in referring changes aloft to equivalent changes at the ground. 
The treatment of ballistic elasticity has not yet been decided upon. 

For each combination of С and V, a chart can be made up showing 
(for instance) the weighting factor curves for range ballistic wind for 
four or five values of G. The curves for other values of & can be inter- 
polated, if only we have acurvefor¢=0. But of course no trajectory 
can be computed for ¢=0. Nevertheless, there has been deduced 
the equation of the limiting curve*which the range wind weighting- 
factor curve approaches as & approaches zero. The equation is: 


98. I- - — Eh л =  - Eh 


in which р is the abscissa and k the n: of any point on the 
curve; and л is the Gávre n, to be obtained from the table of (n— 2) 2 


entered with 100 using the standard muzzle velocity. 
The equation of the limiting curve for density weighting factors is: 
3 xcd 

99. 1—p=7(1—k)*+ 701 -k з 

The equation of the limiting curve for cross wind weighting factors is - 

100. 1—р=(1—Ё) 

The limiting curve for elasticity weighting factors is the same as 
that for density; but, owing to the fact that elasticity weighting 
factors become infinite in certain cases, the use of such weighting 
factors is not to be recommended. 

It will be observed that in the case of density, elasticity, and 
cross wind, the limiting curve 1з independent of muzzle velocity and 
ballistic coefficient. In the case of range wind the limiting curve 
depends on the muzzle velocity, but not on the ballistic coefficient. 

The weighting-factor curves for ¢ = 0 can be plotted on the weight- 


ing-factor charts with very little extra work, and will assist materially 
in the interpolation of new curves among those computed. 


XV. WEIGHTING FACTORS. 95 


By comparing a large collection of weighting-factor curves at 
Aberdeen, three mean wind weighting-factor curves have been 
deduced, whose equations are, respectively: 


1—р=1.11 (1—k)4—0.11 (1—k)? 
101. 1—p=0.74 (1—k)4+0.26 (1—0)? 
1—p=0.36 (1—k)4+0.64 (1—k)? 


For any given battery, one of these curves can be chosen, which 
will give a sufficiently approximate ballistic wind for,both range and 
deflection corrections for actual field service. 

Similarly the following single density weighting-factor curve has 
beén deduced: 


102. 1—p=0.48 (1—k)4+0.52 (1—k)'h 


Prior to the World War a single curve, known as the '' time curve" 
or vacuum curve" was used for all weighting-factor purposes. 
These names are due to the fact that this curve weights the various 
Strata of atmosphere in proportion to the time the projectile would 
spend in each, if the trajectory occurred in vacuo. Its equation was 
the same as that of the limiting curve for cross wind (equation 100). 

The present indications are that the second equation under 101 
will be adopted by both branches of artillery? for all wind weighting- 
factor purposes; and equation 102 for all density-weighting factor 
purposes. | 
PROBLEMS. 


(64) The maximum ordinate. of problem 63 was y,=97.9. Con- 
struct the range-wind and the cross-wind weighting factor curves on 
the same sheet of cross-section paper. 

(65) Plot equations 102 and 100 on a single sheet for comparison 
of the old and the modern mean density weighting-factor curves. 

(66) Plot equations 101 and 100 on a single sheet for comparison 
of the old and the modern mean wind weighting-factor curves. 


QUESTIONS ON CHAPTER XV. 


1. Define ''ballistic wind." 

2. What is the British term for ballistic wind?” 

3. What are weighting factors?” 

4. Explain the two ballistic winds used in present-day methods? Are these the 
components of a single ballistic wind? 

5. As no trajectory can be computed for ¢ = 0, how are the weighting factors 
obtained for this angle of departure? 

6. Of what use are the curves for ¢ == 0? 

7. Are weighting factor curves used to determine ballistic elasticity? 

8. Define ‘ballistic density.“ 

9. What is meant by the time curve"? 


3 As a result of experiments at Fort Monroe, the Coast Artillery has adopted the second of equations 
101 for all wind purposes. 


CHAPTER XVI. 
CONSTRUCTION OF A RANGE TABLE. 


A range firing consists of a number of rounds, usually 10 to 20, 
fired at each of several elevations, say, at 5°, 15°, 25°, 35°, and 45°, 
The exact location of each point of splash is plotted by means of 
intersecting azimuths from at least four observation towers. 

Throughout the firing the meteorological conditions are observed 
from time to time by aeroplane, pilot balloon, etc. At the gun a 
detailed record is kept of the time of firing, quadrant elevation, and 
ot all variations from standard, such as weight ot projectile, cant of 
trunnions, etc. Time of flight is taken by stop watch as a check on 
the computations later to be made and as a basis for the rough deter- 
mination of maximum ordinate (y, —4.05 T). 

Separate rounds are usually fired through jump and chronograph 
screens, as a basis for determining jump and muzzle velocity. 

Given the range-firing records for a specified gun and projectile, 
with prescribed values of weight of projectile and muzzle velocity, 
to construct a range table for this gun and projectile, the procedure 
is as follows: 

Divide the total rounds fired into groups, each аар identical 
ranges (except tor slight differences due to dispersion) because of 
having been fired at the same elevation and on the same day. 

Compute, for each group, the mean values of all measured quanti- 
ties. These measured quantities are muzzle velocity, weight of 
projectile, observed range, deflection, right wheel above left, time of 
flight, ete. 

With average values of the atmospheric conditions, and with 
weighting factors from a similar range table, compute a tentative 
ballistic range wind, cross wind, temperature, and density. If no 
similar range table is available, use the mean weighting-factor curves 
of Chapter XV. 

Compute, for each group, the angle of departure, by correcting the 
quadrant elevation for any individual error in the quadrant, for any 
inclination between quadrant seat and axis of the bore (determined 
by applying a clinometer to the gun) for jump and for height of site. 
(See Chap. XI.) 

Unless velocities were taken on the range rounds themselves, esti- 
mate the mean muzzle velocity of each range group, as follows: Let 
ёр be obtained by substracting from the mean projectile weight of 
the group in question, the mean projectile weight of the velocity 

96 
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rounds fired on the same day. Then the estimated velocity of the 
range group will be the algebraic sum of ôV and the mean velocity of 
these velocity rounds, ôF being obtained from the equation: 


See the explanation preceding formula 58. 

Estimate the ballistic coefficient for each group. This may be 
done in a number of ways. The observed ranges should first be 
roughly corrected for nonstandard conditions, using а similar range 
table, or the Ingalls tables, or Alger’s charts (which are a graphic 
renresentation of the Ingalls correction formulas), or French charts, 
or the A. L. V. F. tables, or the G&vre tables of September 15, 1917. 
The last-named two are French tables. An Americanization of the 
A. L. V. F. tables has been published by the Ordnance Department 
(War Department, Document No. 983, Confidential), and an Ameri- 
canization of the Gâvre tables by the Coast Artillery Board (mimeo- 
granhed—title: Artillery Ballistic Tables“). 

Then with ¢, V, and the corrected .X as arguments, enter '' Ingalls’ 
Ballistic Tables" (printed by War Department, 1918), or a set ot 
French charts, or a set of Alger's charts (Journal of the U. S. Artil- 
lery, Dec. 1919, p. 585), and take out C. It should be noted that 
the C obtained by either of these methods is not the C desired. The 
C of the Ingalls tables and the Alger charts is the Siacci C, so-called: 


2275 
103. С, = Bid? | 


The С of trajectory computations in modern ballistic methods is the 
normal О, so-called: 


104. Cn ich. | 

The coefficient of form, (?) of the two C's is also slightly dif- 
ferent. But a tentative value for Cy can be obtained by extracting 
В and f. from C, The use of French tables or charts is, however, 
much better.! 
8 = sec ¢ 


105. Joss x) = —0.00012 P. 


1 I. is the mean value of H for the trajectery. This is assumed to be the same as the Hat an altitude 
a 


two-thirds the maximum ordinate. The Ingalls formula for maximum ordinate in terms of time of 
flight is y.=4.05 ТЗ. Therefore: 


log „(+ ) logie H= —0.000045 (#) ==—0.00003 . —0.0001215 Тз. 
24647—21—7 
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The French C and the normal C are bound together by the follow- 
ing approximate relation: 


ү Cy + logo Cx = 7.0570 — 10 
Or. Cx = 0.001140. 


The Americanized French tables give Cx direct, instead of Cy. 
Plot the approximate values of C thus found against ¢, and draw a 
smooth curve. For each value of ¢ at which firings were made, take 
the smoothed-out value 


106. 


на tb ty of C and the standard V 

12222 and compute a trajectory 

TEET by the methods laid down 

-----n zan in Chapter VIII. | 

TERTE | Compute the differen- 
SEF LL tial corrections by the 
7 methods of Chapter XIV 
254 for each of the computed 
АШИ trajectories. 


The variations from 
standard for which the 
differential corrections 
are to be computed are as 
follows: 
| aR 1 m/s change in 
© "ANGLE OF ELEVATION IN DEGREES muzzle velocity. 

1 per cent change 
FIG. 4 in atmospheric 
density. 
1 per cent change in absolute temperature (as affecting elas- 
ticity). 

1 m/s range wind. 

1 mil change in elevation. 

1 per cent change in weight of projectile.” 

1 m/s cross wind. 

From these results, compute the changes in range for the variation 
of actual conditions from standard (atmosphere, muzzle velocity, 
weight of projectile, and elevation) at the time of firing, and correct 
the observed ranges to standard ranges (i. e., to the ranges which 
the given ¢ and V, and the estimated C would have produced under 
standard conditions). 


LE ILE SE EH i 


es 
3 Only the first term of formula 58 should be used in getting the standard range from the observed range. 


In getting the estimated velocity of the range rounds from the observed velocity rounds, use en T as 


before. (See the explanation which precedes formula 58.) 


XVI. CONSTRUCTION OF A BANGE TABLE. 99 


Compare these corrected ranges with the ranges obtained from the 
trajectory computations. If the differences are considerable, say, 
greater than 5 per cent, correct the values of the ballistic coefficients 
at which the trajectories were computed and recompute the trajec- 
tories and the differential corrections for the new values of C. Then 
repeat the work indicated in the preceding paragraph. 

. The range-elevation curve should now be constructed. Correct 
each quadrant elevation (i. e., observed elevation) by adding the 
angle of site and the corrections for the quadrant and quadrant seat. 
We do not here correct for jump, as what is here wanted is the angle 
of elevation rather than the angle of departure. Plot the corrected 
ranges against the corresponding elevations as thus obtained. (See 


fig. 14.) 
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Certain elements of the trajectory and certain differential соггес- 
tions should now be plotted against the angle of elevation correspond- 
ing to the angle of departure from which they were computed. (See 
fig. 15.) On one sheet of cross-section paper should be plotted— 

. Slope of fall. | 
Angle of fall (о). 
Time of flight (Т). 
Terminal velocity (vr). 
Maximum ordinate (y,). 


a 


з Both the :omp ite 1l ani the observe i times of flight are plotted. The curve is then drawn, giving by 
far the greater weight to the computed values, but giving some weight to the observed. 
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On one sheet should be plotted (See Fig. 16) the range changes due to— 
1 m/s increase in muzzle velocity. 
1 per cent decrease in atmospheric density. 
1 per cent increase in absolute temperature (only as есир 
elasticity). ! 
1 m/s following ud 
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1 mil increase in elevation. 
1 per cent increase in weight of projectile.‘ 

On one sheet should be plotted (see fig. 17) the deflection effect, 
in mils, of— 

1 m/s cross wind. | 
Drift (including lateral jump). 

Compute the probable error for range and deflection. This is 
obtained by multiplying the mean deviation by 0.845. Usually 20 or 
more rounds are fired at each of several elevations for the determina- 
tion of dispersion alone. The dispersion of the range groups is also 
considered if these contain enough rounds. Plot the 50 per cent zones 
against elevations. The 50 per cent zone is twice the probable error. 
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FIG. 18 


Any additional data which may be needed should be treated in the 
same way. MN 

АН curves should be smoothed out. Practically no smoothing 
should be necessary, except in the case of the drift curve and the 
probable error curves. 

From these curves take the necessary data for the construction. 
of the range table. 

Typical curves are shown in , figures 14 to 18, inclusive, which are’ 
based on the computation of a range table for the 75 mm gun, firing 
a Mark IV projectile at 1,900 f/s. 


4 The whole of formula 58 should be used here. 
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PROBLEM. 


(67) Construct the necessary curves from the following data: 


155 mm G. P. F. gun firing 94.7-pound shell at 2,410 f [s velocity. 


Number of rounds considered 
Clinometer elevation 
Range (meters) 
Deflection (mils )) 
Time of flight (seconds.) 
Weight of projectile (pounds) 
Mean deviation in range (meters). . 
Mean deviation in deflection (mils) 
Ballistic range wind (m/s)......... 
Ballistic cross wind (m/s) 


Ballistic density................... 


Ballistic temperature (°F) 
Right wheel above left (meters). 
Distance between wheels (meters). 
Angle of bore sight 
Height of trunnions above mean 
low water (meters). ............ 
Height of tide (meters)............ 
Vertical jump (minutes) 
Weight of projectile for velocity 
rounds (pounds). . 
Velocity of velocity rounds (f/s).. . 


—— 22 е 


— 


| 


Date. 
June 9. June 9. June 9. June 10. June 14. > 
25 13 12 20 “20 
6° 47 | 16° 35 10235 | 25° 5 | 35° 5 
6,477 11, 288 11,453 | 14,142 | 16,564 
5.9 | 10.7 8.1 6.3 7.4 
13.23 | 29.85 30.06 41.89] 55.0 
93.48 93.67| 33.53 93.52| 93.37 
44 | 54 61 107 94. 
0.5) 0.5 0.6 0.7 0.6 
42.3] 43.4) 43.7 —54| —88 
+181 42.5 +406] —29| —9.0 
0.978 | 0.982] 0.981] 0.972] 0.965. 
70° 70° 70° 65° | . mne 
—0. 0033 |—0.0027 |—0.0027 |—0.0003 | —0. 0027, 
0.503 | 0.503| 0.503) 0.503] 0.503: 
6° 4/| 16° 35’ | 16° 35’ | 25° 5 35? 57: 
5. 94 5. 94 5. 94 5.94 . 5.94 
0.54| 0.40 0.40 0.61] 0.43 
2.0 -20| -20| 2.0 . 0 
93.33 | 93.33 93.33] 93.35] 95. 64 
2,301} 2,391) 2, 301] 2, 401 2,410 


Questions ON CHAPTER XVI. 


l. Give a very general outline of 


the steps in the construction of a range table. 


2. How may the ballistic coefficient be estimated? 
3. Why do you suppose y, equals 4.05 Т?, instead of 4.02 T?? 


4. 
corrected for jump? 
5. To convert the observed range 


In constructing the range-elevation curve, why is not the observed elevation 


to the standard range, should the results of solving 


the formulas of Chapter X be added or subtracted? Why? 
6. What curves are likely to need smoothing out, and why? 


SUPPLEMENT A. 
TRAJECTORY COMPUTATION BY THE TANGENT RECIPROCAL METHOD. 
The tangent reciprocal method of computing trajectories is a 
variant of the rectangular method described in Chapter VIII. It 


is based on the following equations, involving three auxiliary vari- 
ables c, с’, and “: 


/ ) 
9 25 = zm (i. e., tan 0 divided by 9) 
1 (i. e., minus the reciprocal of the horizontal com- 
z’ ponent of velocity) 


c! Ea’ 


These three variables are “critically varying", i. e., they change 
at such rates that inaccuracies produce the least possible effect on 
the results. This enables the use of longer time intervals with the 
tangent reciprocal method than vith the original rectangular method. 

For the tangent reciprocal method, change the trajectory sheet as 
follows: Use the — Er’ column for y. Use the velocity, mean height, 
and time columns, and the first half of the f y dt column for y’. 
Use the last half of this column for time. Use the y column for c, 
the y^ column for ос”, and the - Ey’ —g column for c". 

On the small computing sheet, label the rows as follows (the 
numerals show the order in which tbe rows are used): 


107. 
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t 0 1 2 Order 
ee — ess 8 
I). ⁵ ⁵ Жа» 9 
S sects Constant. 
))..öü деда 10 
E.. — ашыры 1 
logo дедек ра 8 11 
CONG 4... eese uo oS ОАК 2 
NOW жие NA ERE E Constant. 
TOP ER (c SES 3 
log A E E EE 4 
Ihr 5 
УЛОО асрдан: 6 
ohr К Т 
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The logarithms are all denary. 
To start the computation, enter on the trajectory sheet following 
initial values: 


х= 0 
z’= V cos ¢ В 
y = 0 
y’= Уво ф 
; 1 
А dini". 
PL 
„ qe 
g 


Enter on the small sheet log H=, for tO; colog C (constant 


472 ^3 
хорзор log 9 (constant 8 Take 100 and 108 from 


the table of squares; add them to get a 1 © iL and with this as an argu- 


ment, get log G from the G Table. 

Add log G,log H, and colog C to get log E. Add log E to log o’ 
to get logo. Enterc''onthetrajectory sheet. Note that initially c 
is plus, and с’ and 6“ are minus. This completes the computation 
for ¢=0. 

To start a new line, integrate c'' ahead to get the increment of c'. 
Integrate a’ to get the increment of c. 

Turn to the small sheet. Set down logg’. From this, get colog a’. 
Set down log c. Add colog a’, log g and log с to get log у’. ! 

Enter the table of logs and squares with oe z' (ваше as colog с” 


and logy’ as — and take out 2’ „ 100 100 and 100 Add the last 


two to get 100 With this enter the G table and Ew out log G. 


Integrate y’ to get the increment of y; and, with y as an а 
get log H from the formula 


log H= (10 —0.000045y) — 10. 


Add log G, log Н and colog С, to get log E. Add log E and log c^. 
to get log o”. 

Integrate o to get the increment of с’, and е” to get the increment 
of c. If these new values check sufficiently close, the work on that 
line is completed. 

Proceed in the same way with each successive line. 


SUPPLEMENT B. 


EXPLANATION OF THE SIGNS IN THE COMPUTATION OF DIFFERENTIAL 
CORRECTIONS. 


Whatever may be their physical significance, all of the integrations 
of Chapter XIV are, mathematically speaking, performed from a 
later to an earlier time. 

In integrating from a later to an earlier time, the first differences 
of the integral must receive an algebraic sign opposite to that of the 
integrand; in other words, a positive integrand produces an alge- 
braically decreasing integral and vice versa. Thus the integral: 


ih ( ) 44 requires that a change of sign be made when numerically 


A | 
integrating; the integral: — [ ( ) dt, does not. But it is bad 
T 


psychology to use а notation in which plus requires a change of sign, 
and minus does not; hence, although integrating backward, we shall 


T to 
use the symb 1 f in place of the symbol f . Then, since 
to T 


«f^ )ds-- f ) dta, and 
- Г ) dta= f ) dta, 


& minus wil mean to change signs, and a plus will mean not to 
change. . 

Chapter XIV employs two sorts of integrals The integrations to 
get к’ and и are of one sort. Since n“ is the derivative of u then, u 
is the indefinite integral of u’. That is 


в ГА 
pom | д. 


But the values of д are fixed by being known at time T. Hence: 


te te T 
m= | dts u. |да | ий. 


in line 22 of the computation sheet. li 

This equation means that и is the indefinite integral of »’; that its 
value is known at time 7; that to find its value for an earlier time, 
one integrates from 7' to this earlier time; and that one 88 sign 
in integrating. Similarly: 


T 
P PEN "dta u — [| n" dta. 
te 


The other integrals of Chapter XIV are of a different sort, namely 
definite integrals from time t, to time Т, representing the effects, at 
| 105 
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time Т, of causes starting at time to. Thus, in line 40 of the compu- 
tation sheet: 


(8) = + J. (7) dta. 


Logically we ought to integrate from each t, to T, proceeding 
from right to left, as follows: 


10. 16.760 16 14 12 11 10 
(5) . 0 0.08! 0.301 0.52] 0.621 0.73 
n erc err О —0. 22 —0.10 —0.11 |......... 
(8)=+ оаа "m 008] 00 — .... OO — 
te 
84. f. (ууйд....... 0.411 0.38] 0 Ll... eese 
te 
(8)=+ f. (0d. 1.23 1.20 0.820 .Ul . 
(8)=+ it (Ida. 1.80 1.77 1.39 ов) о |......... 
(8)== + E (7)dta....... 2.47 2.44 2.06 1. 24 0. 67 0 
Ete. Ete. Etc. Ес Ее.) Ete. | Е. 


But, if we are willing to give up logical arrangement, in order to 
save labor and space, we can obtain the same results as follows, inte- 
grating from left to right: 


fo. 16.760 16 14 2.5 11 10 


TT 0 0.08} 0.30] 0.52] 0.62 0.73 

ИРИСИ — e — 0.22, 0. 22 0.10 0.11 
T 

(8)=-+ J (7)dta....... 0 0.08| 0.41 1.23] 1.80 2.47 


In either way, the plus sign in front of the integral signifies that 
no change of sign is made in integrating. 

The heavy black line between two columns is & convenient device 
to indicate a change in time interval. 

To recapitulate, although all integrations in Chapter XIV are made 


T 
from time T to time to, we shall use the integral symbol + [ ‚о that a 


minus sign before the integral shall serve notice that the increment of 
the integral is to be given a sign opposite to the sign of. the inte- 
grand, and so that a plus sign shall serve notice that no such change of 
signs is to be made. 


l 


SUPPLEMENT C. 
DIMENSIONS OF BALLISTIC SYMBOLS. 


In deriving ballistic formulas, it is frequently convenient to test 
them by the theory of dimensions. For this purpose we shall employ 
the following three dimensions: Length (Г), time (Т), and mass (М). 
Unity will be expressed by 1. 

The following are the dimensions of the principal ballistic symbols: 


Symbol Dimensions 
"d nU m ] 6⁰ʃ˙mm ½ ͤ AU...... 8 
%% KC p E LIT 
1. and E T L/T? 
————————— L|T 
ОЧОККО ОЛ ОООО K ОЕР: MJI 
ООСО ЕЛУ — TEE 1/1, 
С Tp AA E À———— — MG м/т? 
| COE" p M 
JJV ————— A Seeeke 1 
И hii 
10° СЕ оре .... YL ` 
‚ЖЕН КЕ —— —— À——s À———— kꝛ0h;0 НА 1/T 
Pp —————— Á——————ÉÓ  SÜ€ L 
ЖОО —— —— — ⁵ð L/T 
GET ИККЕ ЕЕН ИЛК ЛЕК И ЕС LIT? 
PONG СРР reieun bee a ТЛ N a T 
WONG E ET EOE essences Л КС Г ieee TT 1 
All trigonometric functions. ........................... eee eee eee 1 
All angles. ............... )jGCCCͥõ;ͥͤõͤ ĩðͤ SEU STE 1 
7; TTE 1 
АШ :ехропевіе осеннее КАКИЕ 1 
All of the foregoing dimensions are fixed by physical laws except 
the dimensions of G, H, and С. ME | | | 
- There is а wide latitude possible in choosing the dimensions for 


these three symbols, provided only that the dimensions of a are Jy 
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Thus C may be treated as a dimensionless constant and H a dimen- 
sionless ratio; G will then have then same dimensions as E, which 
convention would have much to commend it. Or v@ could be 
regarded as having the dimensions of its Mayevski equivalent A, vn, 


n-i 


1. e., 1 whence G would have the dimensions m 
The dimensions actually adopted for G, H, and C were arrived 
at as follows: G was regarded as vB 9 = being treated as а dimen- 


sionless ratio (see equations 59 to 64 in Chap. X). H was given the 
dimensions of density, and C the dimensions of sectional denaity. 
The following physical laws were thus expressed: 


The usefulness of the tabulated dimensions is as follows: To test 
the dimensional eorrectness of a formule, substitute for each symbol 
its dimensions, disregarding algebraic signs and non-dimensional 
numerical coefficients. If the dimensions of all the terms are alike, 
it is dimensionally correct. Very often typographical errors and 
mistakes in derivation can thus be readily detected. 


SUPPLEMENT D. 
ANTIAIRCRAFT FIRE. 


At the time of writing this book the methods of computing dif- 
ferential corrections for antiaircraft fire are in such an unsettled 


state that it is thought best not to describe them in detail, but merely 


to refer to them in a general way. 

In Chapter IX we saw that the use of three auxiliary variables, 
u, », and p, and an auxiliary constant, ^, enabled us to express a 
change in X, T, or y, in terms of charges in z, y, z', and y’ occurring 
at any time, the auxiliary variables being functions of the time of 
the disturbance. The expressions for 6X, ö T, and ду, in that chapter 
may be generalized into 


108. s — )eX&x hy Sv dx’ +p 8% 


where the parenthesis represents an effect of any given nature occurring 
at any given point on the trajectory. Thus 6X is the effect on z at 
the point of fall, ôT is the effect on f at the point of fall, dy, is the 
effect on y at the summit. Similarly we might use the same general 
form to express any given effect at any other point on the trajectory. 
The auxiliary variables here used are not to beconfused with thespecial 
cases considered in Chapters IX, ete. 

In the general form given above, бт, dy, dz’, and ôy’ are arbi- 
trary small increments occurring at the time of the cause (which 
we shall call £4). 

z, y, 1, y’, 4 and y" are functions of the trajectory and of the time 
(t) of the point on the trajectory at which their value is taken, whether 
this be the time of cause (f= ka), or the time of effect (f= te), or some 
other time. 

^, u, », and p are functions of the trajectory, of the nature of the 
effect considered, of the time of cause (ta), and of the time of effect 
Ke). . 
5 ( ) occurs at time te and is a function of the trajectory, of the 
nature of the effect considered, and of бт, бу, dx’, dy’, ta, and te. 

Let us now consider some particular trajectory and some particular 
effect, such as change in z. Let us successively consider some fixed 
value of te, the time of effect, so that it successively equals 0, 2, 4, 
| d. 

Then for each value of te: 


109. дт, = or + hy у-у or +p dy’. 


For each value of te, we can get a different set of values of А, и, >, 


and р as functions of ta alone. 
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The entire computation of Chapter XIV is in effect repeated for 
each value of te, running each computation backward from o = t, to 

—0, the result being a tabulation of differential corrections, from 
which can be determined the effect on the æ coordinate of a projectile 
at any point in its flight, due to a disturbance at any preceding point 
or points in its flight. 

The terminal values used to start each ot these computations are: 


S cot 6 
Е mu 
= (0) 
110. | КЖ" 
p 20 
= (), 


This is a very tedious performance. Accordingly there have been 
devised several alternative methods, each involving variables aux- 
iliary to the auxiliary variables, i. e., bearing much the same rela- 
tion to the auxiliary variables as the auxiliary variables do to the 
elements of the trajectory. 

In the latest method, two sets of values for и and и” are computed, 
each based on the assumption that equation 108 holds true and that 
\ —0; one based on the assumption that Ли is zero and и’ unity at 
the gun, the other on the assumption that hu is unity and џи’ zero 
at the gun. These four variables, called u,, и’, и, and u,, are func- 
tions of the time of cause alone. In this method the integration is 
performed forward. 

Variable coefficients K, H., L., L., M, M, N.,, N., and A are 
also computed, these being functions of u, u., u, and u’, and of 
the time of effect. This A is a quantity, and not an operator. 

Attempts are now being made to simplify this system, Until 
this is accomplished, it is not thought advisable to present the matter 
to students in any more detail than is here given. 


SUPPLEMENT E. 
DERIVATION OF TWO EQUATIONS OF CHAPTER VII. 


The precise equations of motion of a projectile in the tangent 
method (equations 30) or the curved method (equations 31) have 
been derived in various blue prints of the Ordnance Department. 
A brief skeleton of the derivation is here given. 

The equations of motion referred to an instantaneous system of 
cartesian axes, whose origin is the instantaneous position of the 
projectile, and which are respectively horizontal and vertical at that 
point, are: 


z” pu — Ez 
* -—Ey'-g 
111. a =27” 7477’ 3 
tan 0-7, 
T 


The tangent method assumes, in place of an infinite number of sets 
of instantaneous &xes, & single set, namely, the instantaneous set 
whose origin is the gun. The only effect of this assumption is to change 
the basis of the H function and the manner in which gravity enters 
into the equations. Also, in any system, gravity (g) ought to be ex- 
pressed in terms of the constant surface gravity (go). 

If y is the angle at the center of the earth subtended by the flight 
of the projectile up to its arrival at the point X F, if p is the distance 
from the center of the earth to the projectile and E the distance 
from the center of the earth to the gun, and if gravity varies in- 
versely as р?, then: 


P Р сов? 2Y | 
112. 589. 1 9. TE Py nM (1K 3 1 - -) 


The components of gravity are: 


9,72 ein V-, (1-25 + TE )s- erus ) 
Jy=9 cos —9. (IK * 90 T ) 


Therefore the equations of mot'on are as given in equations 30 of 
Chapter VII. 


1 Coordinates z and ӯ are those of the instantaneous system. Coordinates X and Fare t obe of tne 
tangent method. Coordinates т and y are those of the curved method. 


113. 
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Also the correct basis of the H function becomes: 


_R+Y xX? 
114. p—E-———--— 22 * 


The curved method assumes a single system of orthogonal curvi- 
linear coordinates, whose = is measured from the gun along a circle 
concentric with the earth, and whose y is measured vertically upward 
from this circle. 

The equations of motion of the instantaneous system are converted 
into polar coordinates y and р as defined above, and thence into the 
orthogonal curvilinear coordinates. 

The first conversion is derived as follows: Consider & point P, on the 
trajectory near the instantaneous point P, and denote by A y the 
angle at the center of the earth subtended by PP,. "Then: 


r = o Sin A 
115. J =р cos Ay —0P 


where 0 is the center of the earth. Ж 

Differentiate twice, with respect to time, and then let Р, approach 
P. Cos A will approach unity, and sin A y will approach zero, and 
we shall have as conversion equations: 


ў =p’ 
E —2p'y/ + ру’! 
ў''=р''— py” 


The relation between the polar coordinates (y, p) and the coordi- 
nates (z, y) of the curved system is: 


116. 


e 


T 
117. VR 


p=yth 


Differentiating twice, with respect to time, we get as conversion 
equations: 


QU 

y =R 

118, cay 
a 


p! =y" 
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Substituting from equations 118 and 119 in equations 117: 


— z^ 

E =r Р 
119 dii 2 ˙ 
° sth : y 21 7 
4X =T (85). R 

z' " . 

y =Y Р 


Substituting from equations 119 in equations 111, correcting g for 
altitude, but not for obliquity, and neglecting y in comparison with 
R, we get the equations of motion as given in equations 31 of Chapter 
УП. 
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SUPPLEMENT F. 
A DERIVATION OF THEOREM 1. 


An inspection of any number of formulas of ordinary differentia- 
tion will show that they are homogeneous with respect to the ordi-- 
nary differential operator d.! Therefore, let us assume, as is indeed 
the case, that there is no possible form of differentiation which. 
violates this principle. 


Let 
ъ=], (A, F, Z, ) 
X =}, (2, Y, 2, ) 
Y =f, (£, Y, 2, .) 
L=f, (£, Y, 25. ) 


etc. 
Then, by the foregoing assumption: 
du =AdX+BdY+... (a) 
dX=Cdx +Ddy +... (b) 


where A, В, О, D, E, F, etc., are undetermined coefficients. 
Substitution from equations b, e, etc., in equation a gives: 


du=(AC+BE+ ...) dt+(AD+BF+ ...)dy+... (d) 


In equation а, u was considered as a function of X, Y, Z, etc. In 
equation d, v is considered as a function of z, у, z, etc. 
Now divide equation a through by dX. Then: 


du dY 
d; 4t Bayt pi 


If, in this expression, Y, Z, etc., be regarded as constant, then 


becomes oe by definition of the partial derivative, and all the other 


derivatives of this equation vanish. Thus A is identified as 
OUxyz . . . 


Ox 


1 This means that, if we treat d as an algebraic quantity, instead of as an operator, each term will contain 
d to the same power as any other term. 
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Similarly: 


Now divide equation d through by dz, and regard y, z, etc., as con- 
stant. The equation thus becomes: 


Our : 
52 =AC+BE+... 


Substitute the identified equivalents of A, C, B, E, etc. 
Then: 


шуа... duxrz ... OXi;y..., Olxyz... ӘҮ... 2. 
oz ox ` Oz t-5y дт * eu 


SUPPLEMENT G. 
NEW METHODS OF TRAJECTORY COMPUTATION. 


Since this book went to press, a new method of trajectory com- 
putation and a new method of starting computations have been 
devised by a member of the Technical Staff. Thereis room now only 
for a bare outline of these methods, without going into the proof of 
the formulas involved. 


TRAJECTORY COMPUTATION. 


The method of trajectory computation may best be described by 
comparing it with that of Chapter VIII. 

The first four or five lines are computed as laid down in that 
chapter, or by the method hereinafter described. Thereafter, the 
present method makes use of anti-differences of y” and x” 
These bak: written A^! and A-. Thus, if we have & tabulation of A a 
A1, y”, a, b, c, etc. (cf. p. 30), A~ is the first difference of A7; y’ 18 
the first difference of A-!, and hence the second difference of A a ete. 

Now, since, for instance, 25 and 30, 93 and 98, 1000 and 1005, 
all have the same difference, and hence either set could equally 
‘well serve as A`! to produce y’’=5, it becomes necessary to deter- 
mine, in some manner, the initial values of A-t! and A*. By insert- 
ing in formula 121 (post) the values of у’, y’’, and the a, b, and c of y’’ 
for any given line, the value of A^ for that line can easily be 
ascertained. Do this for two or three successive lines; and, if neces- 
sary, adjust the values of A“! thus found, so that the already tabu- 
lated values of “ will be exactly the first differences of A". 

Similarly get the values of A, from formula 122 (post) for several 
successive lines, adjusting as before. Similarly get values of the 
A"! and A? of z”. 

From this point on, proceed practically as in Chapter VIII, except 
that we obtain z', z, y', and y, directly by the formulas of this sup- 
plement, instead of Az’, Az, Ay’, and Ay, as in Chapter VIII. 
Integrate ahead by formulas 33 or 35 of that chapter. For the ten- 
tative value of y, to use in getting log H, we may disregard all terms 
of formula 34 except: 


t 
120 Í y! dt=1 (y'.— За) 
t-i 


Compute the tentative values of y’’ and z'' on the small sheet, as 
in Chapter VIII, and enter them on the trajectory sheet. Obtain the 
new A'!, by algebraically adding the y” of its line to the A~! of the 
line before. Obtain the new A^, by algebraically adding the A 
of its line to the A? of the line before. Similarly for 2“. 


A Technical Staff т, accompanying a Pampo Trajectory, gives additional formulas which 
further lighten the labor of integration i in special ca ne 
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Then, instead of getting the increment of y' by formula 34 of 
Chapter VIII, get y’ itself by: 


1 1 
b.— 27 €) 


i : 1 1 
121. 9% (A 32 / „ 12 434 


and, instead of then integrating у’ to get the increment of , get 4 
itself direct from by: 


122 - m=i (a-a 1p y” — 11 


i t^ 7940 


Note that the first difference, a, does not occur in this formula; 
and that 7 із squared. The first two terms in the parentheses can be 
consolidated into Az. The last term in the parentheses . can 
usually be neglected or estimated. 

. 2’ can be found from formula 121 using z'' and its differences, and 
2 (whenever needed) from formula 122, in each formula substituting 
z's for all the y's. 


No columns of differences oft у", y, z',orz need be carried on the 


trajectory sheet, unless desired as a check on the smoothness with 
which these functions are developing. Also, at least one less column 


of differences need be carried for and 2“ than in Chapter VIII. 


Thus, in place of the 25 columns of figures tabulated on the trajec- 
tory sheet. in the method of that chapter, the pe method can 
get along with 17 columns, as follows: 


„Ass, Аг, z' 5 a, b, Ku PPS 


The chief advantages of this new method are: It necessitates fewer 
columns of figures; its integration formulas are more rapidly con- 
vergent than those of Chapter VIII; fourth differences, which are 
apt to be very erratic in practice, are not employed in this method; 
&ny errors in integrating by this method are.noncumulative, except 
of course as affecting z'' and y'' through the small-sheet computations, 
whereas in the method of Chapter VIII errors have this effect and 
in addition cumulate in their own columns. All that has to be paid 
for this gain is a slight artificiality of method, and the a ыш ss 
labor of starting the A and A 2 columns. Е 


31fa ‘and ¢ are running so smoothly that att, e and cvi can be estimated, the ipllowing gan be subs 
ee for forma 121: 


1 IND 
-3[5. + а. — 12 (0 200 ( ten) ]. : О, — b 


This femula has ọbvious advantages, if the time-interval is two seconds. 
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STARTING THE TRAJECTORY. 


The new method of starting trajectories is applicable either to the 
foregoing method or to that of Chapter VIII, except on very flat 
trajectories. 

The values of the various elements for t=0 are computed аз in 
Chapter VIII. Tentative values for log Ex’ and log Ey for t=i are 
then extrapolated, by subtracting from the values for t=0 the 
respective quantities DL Ez' and DL Ey', found by the formulas: 


123 10* DL Ex =AE+ By’ 
: 10 DL Ey' = 10* DL Ex' +0 


If 5 is other than one second, multiply each of these by d before 
subtracting. | 

The factor 10* is included for convenience in connection with four 
place logarithms. The decrements found by formulas 123 represent 
the change in the fourth decimal place of the logarithm in question. 
The logarithms mentioned in this chapter are all denary, although 
natural logarithms were used in the derivation of the formulas. 

Log A and log B are given on a one-page table (“ Auxiliaries for 


Starting the Trajectory’’) with the argument 100 Dis from this table, 
by means of the argument y’. These symbols should not be confused 
with the three rotation coefficients of Chapter XIII, nor with the 


ballistic coefficient C. 
The following is the explanation of the symbols involved: . 


DL Ee! =$ log Ex- 
A=10' An = 10 0040010 (see p. 64) 
B= U(h+-255,)10" 
00 
M= log e = 0.434294 e е e 
Only a rough interpolation for log A and log B is necessary. For 
& y' less than 100, multiply y' by 10, enter the table with the result 
and multiply by 10 the C thus obtained. 
The tentative values of log Ex’ and log Ey’, resulting from the use 
of the foregoing formulas, yield tentative values of z'' and y”, which 


can be integrated in the manner of Chapter VIII, to у’, y, and z', 
with which to compute the column for t=ż on the small sheet as in 


that chapter. 
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Repeating the same process upon the final figures of the column 
for t=, we obtain tentative values of у’, y, and z for f= i. And so 
on, until enough figures are tabulated on the trajectory sheet to supply 
-sufficient differences for integration ahead. The use of this special 
method can then be discontinued, and the trajectory completed, 
either by the method of Chapter VIII or by the method discussed 
at the beginning of this supplement. 

This special method of starting trajectories makes possible a very 
close approximation of z” and y”, and thus obviates the repeated 
recomputation of the first few time intervals of the trajectory.’ 


з It is claimed that this method will enable the computer to begin a trajectory with one second intervals. 
There is some doubt as to the validity of this claim for very high velocities or very low ballistic coefficients. 
"This method is to be used merely to yield approzimate values of z” and , on which to base the usual 
-method of successive approximations. It is mathematieally equivalent to obtaining correct values for 

a’ and y'" at t=te, and assuming that this value remains unaltered during the interval. It therefore 
materially reduces the number of approximations necessary to get a satisfactory value of z' and . 


SUPPLEMENT Н. 
NOTE ON ADVANCING DIFFERENCE FORMULAS. 


The student’s attention is directed to the fact that, in most math- 

ematical books that treat of interpolation and integration, advancing 
differences are used almost exclusively, instead of the receding differ- 
ences used throughout this book. Therefore, the formulas of most 
textbooks on finite differences will not be those of this book. 
The computations of ballistics that involve interpolation or inte- 
gration are mostly of such a nature that the work continually lies 
at or very near the last one of those values of the functions already 
derived and tabulated. (In the formula for “integration ahead," 
for instance, one limit of our integral, and, in fact, the entire integra- 
tion interval, lies beyond the time of our latest established and tabu- 
lated value of the function.) Hence, we are continually so situated 
that the only differences available to us in this region are receding 
differences. We are, therefore, unable to use the advancing differ- 
ence formulas usually given in textbooks. 

The student who wishes to construct advancing difference for- 
mulas, for use in special cases, may use either of the methods given 
in this note. It is, however, advised that all such formulas be plainly 
labeled to indicate that they are to be used only with advancing 
differences. 

Advancing differences are formed thus: 


а= ї, — f, (instead of f,—f_,), 

а — f, —f, (instead of f, — f,), etc. 

b,— a,— a, (instead of a,—a_,), etc. 
That is, f, — f, is now called a, (instead of d); f; — 2f 4- fo, which, being 
the old а, di, was called b., has now become a, - d= bo; etc. 

It is to be noted that, whether receding or advancing differences 
are used, the subtractions are always made with same direction; i. e., 
we always have f, —f,, da d, etc., never f, —f, or а, dz. 

One method of obtaining any advancing difference formula from 
the corresponding receding difference formula (see footnotes, pages 
17 and 31) is as follows: 

In the formula for receding differences, change the algebraic sign 
of every odd-ordered difference, the even-ordered remaining un- 
changed. Also, change the sign of every t (and therefore of every 
At) throughout the equation. 

The last condition requires that any integral shall have the sign of 
dt and also of both the limits changed. Thus the formula: 

+2 
+ 0 = terms in f, a, b, c, etc. (receding), 
120 
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will become 
—1 


—2 —2 
+ f(- dt) = — /й= 4 fat = terms in f, a, b, c, etc. (advancing); 


the coefficients of f, a, b, c, etc., in one equation, being respectively 
equal to the coefficients of f, a, b, c, etc., in the other equation, but in 
the second equation the coefficients of a, c, etc. (the odd-ordered differ- 
ences), are each opposite in sign to the same coefficients in the other 
equation. Referring to the schedule of page 30, it is noted that if 
the first equation above uses fo, do, bo, Co, etc., of the receding differ- 
ence notation, then the second equation will use fo, do, bo, co, ete., of 
the advancing difference notation, which will be respectively equal to 
Jos а, ba, са, etc., of the receding difference notation of the schedule. 

The preceding process is reversible: i. e., exactly the same changes 
are to be made in transforming any formula for advancing differences 
into the corresponding formula for receding differences. 

A full set of advancing difference formulas can be derived ab 
initio, by tabulating f, f, f, etc., and the advancing differences; 
and then performing the work of problems'28 to 35, and 39 and 41; 
noting that the algebraic signs of all subscripts, limits, and values 
of t in general, in the statements of those problems, are to be changed. 

The expression: “corresponding formula," which has been used in 
the preceding discussion, needs some explanation. If a formula 
uses differences receding (or advancing) from f, to interpolate for 
fn, the corresponding formula in differences advancing (or receding) 
from f, will interpolate for f_n. Ifa formula using differences receding 


n 
(or advancing) from f, gives the value ot | fä the “corresponding ” 
m 


formula in differences advancing (or receding) from f,, will give — fit. 


For example, the formula in receding differences, for integration 
ahead from f, transforms into a formula using advancing differences 
for integrating over the interval back from f,. 

It should be noted that even in Chapter XIV, where the function 
to be integrated is tabulated in reverse order of time, yet, since the 
differences used recede with respect to the order of tabulation, only 
receding difference formulas are to be employed. (See footnote, p.30.) 
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rotation of the earth..................... 80 Rotation of. See Rotation. 

Density: ln tg M —— T 35- 
Standard atmospheric................... 34 | Elasticity correction, same as Temperature 
old and modern assumptions distin- correction, q. v. 

Kulisse.... „ агага унынан 36. 

Effect of nonstandard...... — € 05 range-elevation curves......... € . 98,100 
computation o 88 | Etror, Probable......... S елешн 101 
curves Showing............... сагу ев 100 | Euler............ 7 E e. 09 

Mean weighting factor curve............. 95 } Exponentialfunction........................ 26 

weighting factors: .. 94 | Exponential lav «e sees 43, 65 

Departure: Exterior ballistics: 

Angle of............ ——— 30 Histor.] - 38 
Effect of change in................... 60 Obstacles to development of........... . 38 
Computation of this effect 89 Three parts ok T 

Derivative, Partial. See Partial Derivative. F. 

Descending branch....... — »... 8 sees . 43, 64 

Didion............. TX —— — —— ' FCC EN b esi URS CURE ee, 97 

Differences Fall: 

: Advancing.................... 17 note, 120, 121 Angle оЇ............................- e 35 
formulas for interpolation and inte- Curve of angle of................. РУР . 99 
gration obtainable from  reced- Curve.of slope of......................... 98. 
10g... ness 17 note, 31 note, 120, 121 Рош&фо{...........--+.---.-... S 80 
ditto, derived ab initio...... — 121 | Finite Difference 28. 
not serviceable in ballistics.......... 120 Integration by. See Differences. 
ile aes . 28, 120 Interpolation by. See Differences. 
See Interpolation, Numerical Inte- Fire: 
gration. Anti-áireraft.............. ——ä— V . 40, 109 

Recedin gg. . 30, 120, 121 at target above or below gun............. 40 

formulas for interpolation and inte- High-angle.............. e 39, 40 
gration derivable from advanc- Long range 5 39,40 
RE ——ë 81 note, 121 Low-angle . eeaeee sess. A 

ditto, derived in problems . 30,31 | Firings: 

required in ballistics.............-.. 17, 120 Gávre. ......... TTT SB 8 

reason given 120 iC A 39 

Sign of both kinds of differences the game. 120 1/· A ³·»wuw-w 38 

See also Antidifferences; Numerical In- J7ÜÜĩ ⅛ð K ĩ 8 62, 96 

egrition. Flight, Time of. See Time. 

Differential corrections: Form, Coefficient of. ................ — 38 
Computation of ‚· 9 ꝙ.ä＋³ ... qꝑß . 83 | Forms: 
explanation Of agnus 105 for computation Кн КЕК О A RS 4T 

Differential, Total. ....... ———— M Un for jump corrections. .......... eus WERE 70 

Differential variations e 18 | Formulas: : 
Classification of disturbances............ 19 Differential Variations. See Differen- 
Disturbances proportional to time 26 tial. 

Instantaneous disturbances............. . 20 Integration. See Differences. . 

Varying disturbances.................. — 25 Interpolation. See Differences, Inter- 

Differentiation, formulas apply to variations. 21 polation. 

Differentiation, Partial. See Partial. Franklin, Philip............. 25. 3 — 8 

Dimensional analysis.................. -..... 108 | French Ballistic Coefficient...........-.----- 98 

Dimensions of ballistie quantities € MÀ 107 | French Ballistic Tables............ 3 9T 

Disturbances. See under Differentia] varia- French Charts —' 97 

tions. Functions: 

Drift: Рейпеа................ ———— 10 
includes lateral jump........... ОТРУТ" 74, 75 осы . EEEE waves 44,63 
Measurement oùů . 75 Exponental. ........... 3 3 36, 43, 65 
атта UU JJ ͤ — .... 36, 44 

E. "EIRMOD NEN: к КЕ . 36, 

Е(е!еуаНоп).........................„... 70 Dimensions of........... „ 44, 107 

E (resistance function ꝛꝛꝛõꝛꝛ 36,44 / 32 ͤ see 34 
Dimens. ons of........ 3 ГУТ 107 obliquity......... 8 wines sues . 11 
distinguished from F................... . 43,04 variation with altitude.................. 111 
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Page. Page, 

GAvre: | Metz, Firings at! 38 
firings................... dotado ed Nu a 39 | Mitchell, Prof. H. H......................... 8 
N — ЕВР 9 64 and note. | Motion of projectile treated as a particle... .. 43 
/ РОРИ НЫ ООЛО ГЕНИ .... 97 | Moulton, Prof. F. R 7 

@Gronwall, Dr. T. H.......... EUM 8 | Muzzle Уеіосібу............................. 35 

Effect of nonstandard.................... 60 
H. Computation 0ї..................... 89 

Boro Wiwiesss Er eS FENESTEREREESEER EQUIS cove 36,43 | M. V. See Muzzle Velocity. 

V eee 36, 53 N. 

Dimensions of......... rere „ашалы „„ AUF О, seas ТЫНЫМ ³ / КОНК 24 
Hoar, Capt. Roger Sherman................. San 64 and note. 
Hojel..... % ĩͤ V AS 39 | n (of interior ballistics). ) 62 
Al,. pM 38 |. Newton’s law of air гезізіапсе............... 38 

I Nomenclature; Standard, established........ 8 
` Numerical Integration. (See also Differ- 

i (angle of cant)............ "———— ҮҮТ 71,78 епсез.) 

i (angle of side jump).................... e. 73 adopted from astronomy... ............. 41 

1. efficient of form). ...................- 36, 44,97 called **Short-Arc Method’’............. 41 

Ignoranee, Factor off. 44 Computation of trajectories ys 47, 103, 116, 118 

Impact, Center of.................. Мр кА 36 derivation of ſormu lass 31 

Enclination................ ————— ÁEE , 48, 49, 117 

Increment, Totaallll . 2 p 39 

Ingalls-Siacci methods: Simpson’s Rule. See шша. 

' distinguished from modern methods...:. 39 

Rigidity of trajectory under... .. 68 | Onorator: о. 

Ingalls solution эзе» e 39 des x 12.20 

Ingalls Tables: ыенен а кош А : 
still useful....................-- eee 41 „ . Cee a ete gd 20 

М ted by new t ables... E ИЕС 
.  usein computing range tables........... 97 Co formulas ш 5 а и казо 23 x 
 Entegral, First. нано шш ануе e , 
ntegrating factor (8)............ META сот р Ordinate, Maximnm. i. 4e шыннан, x 

integration: Approximate formula for................ 96 
Formal e „„ 5 9 
"edd Differences, Numerical Osgood, Differential ond Integral Calculus .. II note, 

Interior balllstiecss 2.00000000000000 .. 62,97 A nose: 

Interpolation: P. 

Formulas for.......... — — Wess 17 51 | P...... 2222 8 24 

Derivation of ................ 81, 50,120,121 | p..................-..-. 535 ͤ Ä 

Slide- rule methods. ................... 51 note. Partial derivative: 
Introduction КИРЕРГЕ 7 ambiguity of the symbol................ 10 
. J subscript notation....................... il 
р | Рага1@ШегепНаЦоп...................... - 10 

PII ое... 7 general theorem 11 

Jump: / ⁰˙· AA vote . 114 

Lateral.................. 9 2 •—— — 73 special formula kk 12 
Included in drift........ — 74 | Physical Bases of Ballistic Tables Compute 
Vertical. — —(— . 36 | ion... ... notes on pages 40, 63, 68, 76, 81 
Measurement ої. .................... 70 | Pivot: 
Effect of, on er C 70 Of cant. coner 73 
Pivot oll... — 71, 74 of jump eee eee „ 71,74 
L. анан of projection 37 
par nt: 

L / NE ER NM FERME 24 КТИТОР an 35 

: M. Of Splash u.c Lieu ð адаар 36 

M ee oe vui ТОГОЛ da a EM Dua Cus 24 | Primes represent time-derivatives.. ......... 18 

Mathematical requirements for a reader of Probable Error........... —— СҮЗ "— 101 

n ã ..... 8 9 | Projectile: 

Maximum Ordinate......................... 36 Motion ОС ioo see sna br e EN EOM ERAS US 43 
Approximate formula for................ 96 Variations in weight of... ............ w.. 61 
Gurs... EA 99 Standard, defined........... N 43 note. 

Mayevski’s formula for atmospheric retarda- Projection: 

rr... y sed eeu Meu 39, 64 Plane of....... 74644424 . — 37 

Meppen, Firings at.......................... 39 АП @ k...... Os Es Es 85 
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Q Page. , Standard: Page. - 
Quadrant: r 8 34 
angle of depart ur 35 ргоўееШе...........-..-.-.-.-...---.----.- 43 note. 
ae 35 temperature. 34 
elevation 30 {га]ес{огу............--------...------.- 34 
Ll PP 96 | Successive approximations. . ..............-- 14 
advantages of the method............... 16. 
R. distinguished from mere succession of 
P p CE ⁵ ⁵ĩ— 88 36 approximationn sss. 41 
Range e eee cece eecees 35 | Symbology, Standard, established.......... 8 
Range-elevation res... 98 | Symbols, defined: | 
Range: ЖОЛЧО vowing ANSA VERE ss ONE UE se .. 81,88, 89° 
effects rather than corrections tabulated . Р. e — ͤ FFC 64 note 
effects in general........................ 59 B (function of velocity).................. 44, 63 
formulas for effects of: B (rotation coefficient)............... 81, 88, 80 
change in acceleration. . ............. 61 C (ballistic coefficient)... .. Saber cust aea 86, 97 
change in angle of departure. 60 C (rotation coefficient)...... — 81, 88, 89 
change in muzzle velocity 60 )))) ³·³AA ³ 8 98 
complementary ste 68 [MH 97 
curvature of earth 67 6 . S We Ea nel ГЕ 97 
nonstandard ballistic coefficient. .... 61 D. oM d SEI RE INE TAI 18: 
nonstandard density................ 65 ҮР МЛГЕ ⁰⁰ С С 15 
nonstandard temperature (elasticity). 63 B/!» . 8 81 
nonstandard weight of projectile. . ... 61 d (diameter) РЕА JJ ⁵—8 36- 
rear ул/їп@..................-......... 66 d (operator)................ 8 20 
rotation of the earth................. 80 E (elevation) .......................-.... 70 
I! УСЛУГ 68 E (resistance function) .................. 36, 44 
vertical wind........................ 07 p MP -——————— —— € 48,04 
WING eo ee seo КЕ эре EE 67 “ЖОРУККА ГС SAU ENTM 8 9T 
Range-firinnnn gg.... 96 ice Sedat ce АКЕ UU UR 36, 44 
records, how used.. ..................... 96 | rop E —— РОНА 34,80, 111 
Range-table, Construction of... ............. 96 ?.—— À € e 36, 43 
Receding differences. — 31 J ХОРА „ 36, 58 
See also Differences, Numerical Integra- i (angle of cant). — — — 11, 7 
tion. . (angle of side jump ..................- 73 
Rectangular method of trajectory computa- i (coefficient of form)... ................. 36, 44 
ssd“ 8 47 р ER ——M e nuiosa 71 
Retardation... .......... cc . 36 | MER NR 24 
Approximate expressions for.. .......... 38 М........... JFF! аа Ку ОЛАН — 24 
can not be algebraically expressed. 38 N RE M —— СОГОТ 24 
Mayevski's formula for.. ................ 39 n (@йх+тгө)..............--.--..---- 64 and note 
“Rigidity of the trajectory,’’ Theory of...... 40, 68 n (ofinterior ballistics) .................. 62 
Ritt; Dr. J. Ec 8 3 FFC — 24 
e ß Le ECKE UI SUE 38 v E m 62 
Rotation coefficients. 81 J7J½V%VCCCͥõͤĩͤĩ;k ͤKx.ß. ences Due T duh 96 
Computation ooo 88 8 ааваасаа нажах 63- 
Rotation of the eart gn... 76 о PENES аге 36- 
coefficients. ...................... eee eee 81 ü S 19, 35. 
computation ot 88 F a E AAV СЕ Ue 19, 85 
deflection effect.......................... 80 | AEE TR Em 35- 
Effect of— ONE Inn x И ү чө 35- 
in curved method................... 80 J DU rci UL LER qd E 36- 
in tangent method.................. 79 Б. ТОРЕ nape NCR POR RATE Sener 24 
Range effect of. ......................... 80 An 35 
т and y, z and y, and X and Y dis- 
B. tinguished........................-. 111 note 
Short-are methods 4l and eee I 35 
Siacci solutilo nnn 39, 41 z" and ũꝶ· .. 35 
See also [ngalls. ГООГО ТО аем 36 
Signs, Algebraic, explained.................. 105 Т EEE DA E esce ЛАЛУУ СС N 
Simpson’s Rule... dc. 29 AX oibo РЕР - 61 
Checking dd 50 AX COP РРР 66 
Proof of S277 ае а өөө Ө Vie өс eV ai" ее» P 29 AX •vlhllllllllll 92-972» ое . 67 
Site, Angle оЇ....................:........... 36 C)] —H— р нан ez 
согге°Ноп...............-................ 68 AX SL ² PH 00 
Splash, Point of............................. 36 / ]˙Ü ] AAA НОРИИ 67 
Square root, Method for extracting. ......... 15 JJ/%%% yt 8 TO 67 
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Symbols—Continued. Page 
%% T 68 
AX 9 9ꝙ＋ꝗꝓ9 VVV terres 60 
Гар. СОСЕН РЕТИ 8 aga ME . 81 
JJ.“ ³ . 20 
C ————— —Ü 68 
7 ĩↄðſ ³ĩðA 8 35 
À, M, ”, rr 53 
/ AE y 8 63 
S 35333 ͤͤ ĩͤ b О 85 
r EE ТК y a EE 35 

Symbols, Dimensions of.................... 44, 107 

T. 

; c ENS 36 
Nep m 19, 35 
ТОКТОЛСО 19, 35 
Tables. See specific subject matter. 
Tangent-reciprocal method .................. 103 
Temperature, Standard...................... 34 

Effect of nonstandard ................... 63 

Computation of...................... 88 
Curves showing...................... 100 

Time of flght........................... e eese 36 
Curves Of Pr PE 99 
observed by stop-watch................. 96 

E,, 00000000000 00000s0000o . 18 
Computation of. See Computation. 

Rigidity of................... „ 40, 68 

Sn, 34 

V. 

| pc rM 35 
N ЕАР РТР РЗК C 85 
Variables: 

Auxiliary. See Auxiliary. 

Critically varying........ esse ..s.......- 103 
Variations, Differential. See Differential. 
Velocity, Muzzle. ................... баан 35 

Effect of nonstandard ........... 25006206. QU 

Computation of...................... 89 
Curves showing....... РГР — 100 


| Velocity, Muzzle—Continued. Page. 
Remaining.................. 8 35 
Terminal, Curves ll... 99 
W. 
CCC ò ͤ P EHE AR AV SR AA RBS seen 36 
Weight of projectile, Effect of nonstandard.. 61 
Curves showing this effect. 100 
Weighting factors: 
Cross wind............... eee ce cer eee 94 
Curves for zero $....................... 94 
A/, e erue ee une ua. ARE EE eS. 94 
mean eur ð 95 
Range wind e ... 91 
Wind: 
outer оосар — ..-. 91, 94 
Computation of effects G ... 88 
Curves showing effects ooo. eig des 100 
Effects of cross wind.................. e. 74 
Effects of rear wind................... . 66 
Effect of vertical wind........ — 67 
Equivalent uniform 91 
mean weighting faetor curves 95 
old апа new methods of treating wind 
N. ³ð ⁰ КУЧКА 40 
weighting factors РРР ЕР —-—— 91 
X. 
J) 4 24, 111 note 
r E A E OE TTE E EET 35, 111 note 
Eis C e e note 
J ³¹A—⁴A EERS "— 60 
Y. 
КРРТ N RE —(—— Habe $$ 
7 e 22 . 35, 111 note 
JJ ⁵ ↄ §S3õ· x 
55 e e e e es 
J ⅛ðÄVu Г e — 35 
y cee — 8 e . 35 
2. 
Zaboudski................... — (M а 39 
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